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SUMMARY 


The primary objective of this study was to establish an accurate method to 
predict the cyclic-life of high chamber pressure regeneratively cooled rocket 
engine main combustion chambers constructed of a channel configuration. This 
objective was to be accomplished by a refined analytical model that accommo- 
dated the experimentally observed cyclic-creep phenomenon as an additional 
life limiting factor. This cyclic-creep phenomenon is characterized by mid- 
channel wall thinning and hot wall deformation which is cyclic dependent but 
independent of engine mainstage duration. The basis for this study was the 
channel deformation and wall thinning observed during hot-fire cyclic-testing 
of the Space Shuttle Main Engine-Main Combustion Chamber (SSME-MCC) and other 
combustion chambers tested at NASA-LeRC and NASA-MSFC; as typically noted in 
Fig. 1. 

This study consisted of (1) determining the cyclic life of the baseline SSME- 
MCC based on measured calorimeter heat transfer data, (2) determining the life 
sensitivity of local hot spots caused by injector effects, and (3) assessing 
life enhanced redesigns. The overall program methodology applied for this 
study is presented in Fig. 2. These analyses utilized an upgraded analytical 
model which accounted for cyclic creep. Cyclic channel deformation and wall 
thinning geometry were updated for each analysis to account for the accrued 
distortions. The cyclic life results of these analyses were then compared to 
the typically used low-cyclic-fatigue (LCF) analyses that do not address the 
cyclic creep phenomenon. 

The analytical model was characterized with the basic SSME-MCC channel geome- 
try. First-order sensitivity studies were performed between a nominal full- 
power-level heat transfer condition (100% hg) and a condition of 135% hg that 
may be related to an extreme local hot spot condition over a width of five or 
more channels. The analytical model was further sensitized with artificial 
pressure gradients and artificial temperature gradients across the channel 
wall. 

The results of this study show that cyclic creep is the dominating factor in 
defining the reusable life of a combustion chamber. It is typically shown 
that the usable cyclic-creep life is approximately 30 to 50% of the commonly 
used low-cycle-fatigue life. This reduction in usable life is basically 
attributed to (1) the accrued channel wall thinning as a result of cyclic 
creep and (2) the ultimate strain or tensile instability criterion. The range 
of biaxial tensile instability was determined to be between 5 and 8% wall 
thinning, which is approximately equivalent to the average effective strain. 

The analytical model depicts the mid-channel hot wall being deformed slightly 
inward toward the combustor centerline as the channel wall is thinned due to 
cyclic creep. The degree of channel wall thinning relative to hot wall defor- 
mation is characterized by the channel geometry and surface temperature. The 
lower wall temperature condition reflects less wall thinning for a given chan- 
nel wall deformation or number of duty cycles. These analytical results are 
supported by experimental combustion chamber data which have indicated varying 
degrees of wall thinning for similar hot wall deformation. The analyses also 




Figure 1. Throat Section of SSME-MCC Showing Thinned 
Hot-Gas Wall 
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STUDY PROGRAM INPUTS 



Figure 2 . Overall Program Methodology 
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show that continued wall thinning persists after hot wall deformation is cur- 
tailed. This appears to occur in the projected range where tensile instabil- 
ity is initiated. This characteristic has also been experimentally evidenced 
on hot-fired combustion chambers. 

Four potential life-enhanced designs were analytically assessed at the extreme 
135% hg thermal condition. These designs stemmed from the cyclic-creep sensi- 
tivity analyses of the baseline SSME-MCC channel geometry. These designs 
spanned a wide spectrum of channel geometry and manufacturing feasibility. A 
maximum usable life benefit of a factor of three was shown for the maximum 
risk design which incorporated two life enhancement features. The next best 
design is considered currently feasible with moderate risk and reflected twice 
the usable life. The other two designs did not provide any usable life bene- 
fits and consisted of a moderate risk design and a minimum risk design with a 
minimum SSHE-HCC design change. 

The primary technology gained from this study is the establishment of a cyclic 
creep life prediction methodology along with a working analytical model capa- 
ble of evaluating innumerable channel wall design concepts in a cost effective 
manner. 
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INTRODUCTION 


The reusable life of the Space Shuttle Main Engine (SSME) and future engines 
are greatly influenced by the cyclic life of several major components that are 
exposed to a high temperature combustion gas environment. The main combustion 
chamber (MCC) liner is typically exposed to the most severe environment. The 
SSME and future engines operate in excessof 3000 psia chamber pressure at 
approximately 6500 R combustion gas temperature. These conditions produce a 
heat flux of approximately 100 Btu/in.^-sec in the life limited throat 
region. The thrust chamber reusable life is attained by regenerative cooling 
a copper base MCC liner configured of rectangular cooling channels. During an 
operational duty cycle (mission) the MCC liner is subjected to a large tran- 
sient thermal gradient which imparts a high thermal cyclic strain to the MCC 
liner hot gas wall. 

Typically, low cycle fatigue life is a function of the cyclic strain range, 
the material properties, and the operating temperature. The reusable life is 
normally defined by the number of strain cycles that can be accrued before 
severe material degradation occurs. Reusable life is normally signified by 
the initiation or propagation of surface cracks. Hot-fire testing of channel 
wall combustors at Rocketdyne and NASA-LeRC has shown significant mid-channel 
wall thinning or deformation during accrued cyclic testing. This phenomenon 
is termed cyclic-creep and appears to be significantly accelerated at elevated 
surface temperatures. The sensitivity of this phenomenon to surface tempera- 
ture is evidenced by the non-uniformity of channel wall deformation around the 
circumference of the combustors. History has shown the non-uniformity is 
first evidenced by local discoloration (blanching) of the copper base alloy 
hot gas wall signifying higher surface temperature hot spots. Hot-fire test 
history also shows that minor changes in operational level significantly 
degrade the cycle life in these hot-spot regions. 

The purpose of this study was to provide an analytical model that accounted 
for cyclic creep channel deformation and establish the methodology necessary 
to upgrade the SSME-MCC design and future engine designs. The analytical 
model was to provide the capability of accounting for the coolant channel wall 
deformation cyclic history. 

The approach taken to accomplish the study objectives consisted of three tech- 
nical tasks as follows: 

Task It SSME Main Combustion Chamber Life Analysis consisted of performing 
an SSME MCC life analysis at standard operating conditions that considered 
cyclic-creep damage as an additional life-limiting factor. Analyses included 
thermal, structural, and failure mode analysis. Two-dimensional thermal and 
structural analyses were performed in the critical throat region. The struc- 
tural analysis was to include elastic-plastic behavior of structures subjected 
to constant and variable amplitude cyclic mechanical and thermal loads, non- 
linear material properties variation as a function of temperature, and non- 
linear geometric updating. The combined analysis was to identify the failure 
mode and predict the combustion chamber life. 
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Localized Temperature Ext r emes on SSME Main Combustion 
g ^ber Lite was investigated to define the impact of local hot-gas-wall tem- 

Ls simii injector anomalies on reusable life. The analysis 

of thrsSME wltt ll i" Task I. The duty cycle was identical to that 

the SSME with the exception of the varying local thermal conditions. Para- 

define the cyclic life sensitivity to local 
heat flux and material temperature variations. 

^ ^four^'vial?!f Vf consisted of identifying a minimum 

rAB.fifc ^^^^®-®"*'ancing main combustion chamber designs based on the 

^ L®"'* selected designs analyzed utilized min- 

CLn. / existing SSME main combustion chamber. One of the four 

from feL S .® ®®’" *=°"®®Pt (maximum life), and deviated extensively 

oroLcArt current design. The remaining design was a combination of the other 
proposed designs. The use of hot gas wall coatings was excluded from the 
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DISCUSSION 


THERMAL ANALYSIS 

The reusable life of the Space Shuttle Main Engine (SSMB) is greatly influ- 
enced by the cyclic life of several major components that are exposed to a 
high temperature combustion gas environment. The SSME main combustion chamber 
(MCC) liner is exposed to the most severe environment. It operates in excess 
of 3000 psia chamber pressure and 6500 R combustion gas temperature. These 
conditions produce a nominal heat flux of 100 Btu/in.^-sec in the throat 
region. 

The SSME-MCC reusable life is attained by regenerative hydrogen cooling the 
copper base (NARloy-Z) MCC liner configured of rectangular cooling channels. 
During an operational duty cycle (mission) the MCC liner is subjected to a 
transient thermal gradient approaching 860 K (1500 R> . The peak SSME MCC heat 
flux region, located 3.1 cm (1.22 inch) upstream of the geometric throat, was 
selected for this study. A two-dimensional thermal model was used to define 
the life cycle transient temperature history as a function of time. The MCC 
average axial thermal conditions were determined from hot-fire calorimeter 
chamber measured heat transfer rates. These type of measurements cannot dis- 
tinguish circumferential heat transfer variations commonly associated with 
injection/combustion abnormalities. SSME-MCC regenerative cooled cyclic life 
testing has shown that failures occur in local hot spot regions. These hot 
spot regions are identified by hot gas wall surface discoloration (blanching) 
and typically exist over a width of 5 to 10 channels. 

This section presents the thermal analyses for: 

1. The baseline SSME-MCC sensitivity study between 100% FPL hg level and 
1A0% FPL hg level. 

2. The baseline SSME-MCC channel wall thinning sensitivity study between 
the nominal channel wall thickness of 0.071 cm (0.028 inch) and 
deformed mid-channel wall thickness of 0.010 cm (0.004 inch). 

3. Four life-enhanced channel redesigns, evaluated at 135% FPL hg level. 
SSME-MCC BASELINE DESIGN 

The baseline design for this study was the SSME-MCC. The SSME Main Combustion 
Chamber is shown in Fig. 3. The selected baseline operating level was Full 
Power Level (FPL) with the heat transfer rates noted as 100% hg. This power 
level is 9% above the SSME Rated Power Level (RPL) . The FPL chamber pressure 
is 3250 psia and the MCC mixture ratio is 6.0 (oxygen/hydrogen). The SSME 
duty cycle, including the start and cutoff transient conditions, were used for 
all analyses. The SSME-MCC chamber pressure and mixture ratio duty cycle are 
depicted in Fig. 4 and 5. 
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The minimum life region is located slightly upstream of the geometric throat. 
The peak heat flux and maximum thermal gradient across the MCC NARloy-Z liner 
occurs at -3.1 cm (-1,22 inch). The two-dimensional thermal conditions at 
this location were used for all analyses in this study. 

Gas Side Heat Transfer 

The SSME-MCC heat transfer rates were based on subscale 40K thrust calorimeter 
chamber testing. The subscale test article is shown in Fig. 6. It employed a 
staged combustion system and utilized main injector elements identical to the 
SSME. The combustor length (injector to throat) was identical to the SSME- 
MCC. This water-cooled calorimeter chamber provided local axial heat flux 
measurements over a range of chamber pressure up to 1750 psia. Chamber pres- 
sure scaling criteria were attained from the 40K calorimeter MCC and subse- 
quent testing on two 40K hydrogen cooled chambers up to 3000 psia chamber 
pressure. These hydrogen cooled chambers were configured of the SSME-MCC 
channel geometry. Comparison of local heat flux measurements and total heat 
load data provided local heat transfer scaling to 3000 psia chamber pressure, 
which adhered to the classical scaling relationship of 

hg cc (Pc)^*^ 

Since the purpose of the two 40K hydrogen cooled chambers was a life demon- 
stration, these combustion chambers utilized the identical SSME-MCC throat 
region channel geometry and hot gas wall contour. The local gas-side heat 
transfer conditions of the full-size SSME-MCC were obtained by correcting for 
local hot gas wall mass velocity (hg (pV)^*^) where the SSME-MCC axisym- 
metric geometry differend from the 40K-MCC. This correction is effectively 
attained by integrating the measured 40K calor imeter-MCC heat transfer data 
with the Rocketdyne analytical boundary layer computer program. This allows 
evaluation of any chamber and nozzle-geometry on an empirical basis. The 
injector end heat transfer conditions are based on experimental data. This 
empirical approach was subsequently verified by measured SSME-MCC heat load 
measurements . 

Coolant Heat Transfer 


The hydrogen convective coolant coefficient was based on the Rocketdyne hydro- 
gen coolant correlation as noted below: 


he = C„ G C (j5, a 
w 


where 


Tg = bulk temperature 
T^ = wall temperature 
Cp = heat capacity 
G = mass flux 
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0c = curvature enhancement 

C - 

H 0.92 + ^/FTS l6(c*) - 8.48J 

where 

g(c*) = 4.7 fop g* >70 

g(c*) = 4.5 + 0.57 c*0-75 for c* < 7.0 

A coolant curvature enhancement (0c> of 1.4 was used in the life-limited 
throat region at X = -3.1 cm (-1.22 inch). This enhancement was expreimental- 
ly evaluated by laboratory tests for the SSME-MCC geometry. 

The thermal computer-program logic diagram used in this study is shown in 
Fig. 7. The temperature distribution of the coolant channel is determined by 
a sophisticated 277 node dimensional grid model as shown in Fig. 8. The model 
uses a finite difference technique. The nodal temperature distribution is 
stored on tape and subsequently used with the structural/stress analysis com- 
puter code. The same nodal network was used for the thermal and structural 
analyses. The coolant channel thermal distribution is shown in Fig. 9 at 
3.1 cm (1.22 inch) upstream of the geometric throat. The maximum hot-gas wall 
temperature is 866 K (1100 F) at steady-state FPL (100% hg) conditions. The 
coolant bulk temperature and combustion gas temperature profiles are shown in 
Fig. 10. Noteworthy transient and steady-state thermal conditions of indi- 
cated analytical model nodes (Fig. 11) are shown in Fig. 12 over midchannel 
and Fig. 13 over midland. 
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APSA/STAMP 

• Use finite element numerical techniques to 
generate mesh grid for MCC channel wall at 
both deformed and undeformed conditions. 

• Define each nodal volume, heat conduction 
area and heat transfer resistance. 


ROCKETDYNE BOUNDARY LAYER COMPUTER PROGRAM 

• Analyze the combustion chamber boundary layer 
characteristics at subsonic, sonic, and super- 
sonic flow. 

• Calculate the hot gas theoretical convective 
film coefflc lent , ad iaba t i c wa 11 temperature 
and pressure along the MCC liner wall. 


Incorporate 40K calorimeter 
thrust chamber te^t data. 


• Provide 2-D channel wall temperature 
distribution at transient and steady- 
state condition for simulating the 
"start-up”, "Full power," and "Shut- 
down" engine conditions. 

• Store all MCC channel wall temperature 
distributions at each time slice on the 
DEAP stress tape for stress analysis 
purposes. 


REGENERATIVE COOLING DESIGN/ANALYSIS COMPUTER PROGRAM 

• Analyze the regenerative coolant flow behavior 
along the cooling jacket. 

• Predict the thermodynamic and transport properties 
of the coolant fluid at each engine operating 
condition. 

• Calculate the coolant film coefficient, coolant 
pressure drop, and bulk temperature rise. 


Figure 7. Thermal Computer Programs Flow Chart 
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Figure 8. SSME MCC Cooling Channel Geometry and Nodal Structure 
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Figure 9. Coolant Channel Thermal Distribution at 100% FPL hg 
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Figure 10. Combustion Gas Temperature and Hydrogen Coolant 
Temperature Duty Cycle Profiles for SSME at 
X=-3.1 cm (X=-1.22 Inch) 
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Figure 11. SSME MCC Coolant Channel Geometry 
at X=-3.1 cm (X=-1.22 Inch) 
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Figure 12. SSME MCC Liner Midchannel Temperature Response 
at X=-3.1 cm (X=-l,22 Inch) at FPL Operation 
(100% FPL hg) 
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Figure 13. SSME MCC Liner Midland Temperature Response 
at X=-3.1 cm (X=-1.22 Inch) at FPL 
Operation (100% FPL hg) 
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HEAT FLUX (HOT SPOT) SENSITIVITY 

Hot-fire testing of channel wall combustors at Rocketdyne and NASA-LeRC has 
shown significant midchannel wall thinning or deformation during accrued cyc- 
lic testing of regeneratively cooled combustors. This phenomenon is termed 
cyclic-creep and appears to be significantly enhanced at elevated surface tem- 
peratures. The sensitivity of this phenomenon to surface temperature is evi- 
denced by the nonuniformity of channel wall deformation around the circumfer- 
ence of the combustors. History has shown the nonuniformity is attributed to 
injector effects which may not necessarily be detectable by laboratory flow 
checks or physical measurements. This nonuniformity is first evidenced by 
local discoloration (blanching) of the copper base alloy hot-gas wall, signi- 
fying higher surface temperature hot spots. MCC regenerative cooled life hot- 
fire testing has shown that failurs occur in these circumferentially located 
hot-spot regions. The SSME-MCC axial thermal conditions were determined from 
measured 40K subscale hot-fire calorimeter chamber heat transfer rates. How- 
ever, these type of measurements cannot distinguish circumferential heat 
transfer variations commonly associated with injection/combustion 
abnormalities . 

Two hot spot thermal sensitivity studies were conducted at the life limited 
axial location -3.1 cm (-1.22 inch) to provide insight to the subsequent life 
sensitivity evaluation. These consisted of: 

1. Heat Transfer (hg) Sensitivity 

Thermal (hg) assessment of nondeformed channels between an 
hg level of 100% FPL and 140% FPL was conducted. The 
influence of the gas side heat transfer (hg) level on the 
MCC hot-gas wall temperature over midchannel is shown in 
Fig. 14. The coolant channel thermal distribution is shown 
in Fig. 15 through 18 for hg levels of 110%, 120%, 130%, 
and 140%. The hot-gas wall temperature sensitivity to hg 
level is about 65 K (120 F) per 10% variation in the heat 
transfer rate. 

2. Midchannel Wall Cyclic Creep Distention 

Thermal assessment of midchannel wall thinning sensitivity 
was conducted at the 120% hg level. These results are sum- 
marized in Fig. 19 which depicts midland and midchannel 
hot-gas wall temperature and mid-channel coolant wall tem- 
perature as a function of midchannel wall thickness. The 
channel wall geometry utilized for this preliminary assess- 
ment and the resultant thermal distributions are shown in 
Fig. 20 through 23 for four cyclic-creep distention levels. 
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Figure 15. Coolant Channel Thermal Distribution (110% FPL hg) 
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Figure 17. Coolant Channel Thermal Distribution <133% FPL hg) 
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MINIMUM WALL THICKNESS, INCHES 



Figure 19. Influence of Midland and Midchannel Hot-Gas Wall and 
K) Coolant Wall Temperature as a Function of Midchannel 
^ Wall Thinning for 1207. FPL hg 
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Figure 20. 


Cyclic-Creep Distention Isotherms at 
Minimum Wall Thickness = 0.056 cm (0.022 Inch) for 120% FPL hg 
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LAND WIDTH 
WALL THICKNESS 


CW - 0.040 
LW - 0.045 
WT - 0.028 


0 * 1100 
P - 1150 


CHANNO. WIDTH INCH 


Second Channel Geometry Cyclic-Creep Distention Isotherms at 
Minimum Wall Thickness = 0.041 cm (0.016 Inch) for 120% FPL hg 




Figure 22. 


Cyclic-Creep Distention Isotherms at 
Minimum Wall Thickness = 0.025 cm (0.010 Inch) for 120% FPL hg 
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LIFE ENHANCED DESIGNS 


life-enhanced designs were assessed at the 135% hg level. This 

ihi« S T} effected hot spot condition 

his elevated thermal level may be related to the SSME-MCC liner thermal con- 

ditions near th^e time of incipient cyclic life channel failure and probably 

n f^rthrthermirsen roughening. These designs stem- 

til analyses and the cyclic-creep channel disten- 

Iririevel baseline SSME-MCC as related to hg level and 


The selected designs, noted below, span a wide spectrum of channel geometry 
and manufacturing feasibility. The first two MCC liner designs can alomo- 
SSME-MCC without any changes to the inner mold-line (same hot 

1/Ni IloHrVh-^2'^ mold-line (same structural jacket with the same 

Sclmnl? 7 thickness). These two designs require minimum to moderate 

technology advancements. The third and fourth designs require changing either 
the inner or outer MCC liner mold-lines. Either change would be cltir 
molll*"^ combustion chamber tooling plus a new structural jacket (outer 
fill hoT accepting a smaller throat diameter with a slightly modi- 

eas wall contour. The third design requires moderate technology 

advancement while the fourth design requires entensive technology advancement^^ 


1 . 

2 . 

3. 

4. 


Contoured Midchannel Coolant Wall 
Increased Number of Channels 
Keel-Rib (Midchannel Fin) 

Combined Slotted Hot Gas Wall 
channels 


Midland with increased number of 


These four life-enhanced designs are pictorially depicted in Fig. 24 and com- 

rrcm (°1 22%nch?^uL% located 

J.x cm (1.22 inch) upstream of the geometric throat. Pertinent dimensions 

ZT\ll and thermal operating conditions are noted for 

the'sSME-MC?wtv''°"?"^'°"® 'represent a constant coolant AP equivalen^t to 

straints al^^ d^Tsen^rerf considerations, guidelines, and con- 
straints are discussed separately for each design as follow: 


GUIDELINES 


• Minimum Channel Width 

• Minimum Land Width 

• Minimum Wall Thickness 

o Maximum Channel Height/Width 
o Coolant Liner Material 

• Hot Gas Wall Contour 
o Coolant Pressure Drop 


0.076 cm 0.030 in. 

0.076 cm 0.030 in. 

0.063 cm 0.025 in. 

5:1 

NARloy-Z 
Same as SSME 
Same as SSME 
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• THERMAL CONDITIONS ARE FOR 135% 
OF FPL hg EXCEPT AS NOTED FOR THE 
100% FPL hg SSME DESIGN CONDITION 



WALL THICKNESS 
CHANNEL WIDTH 
LAND WIDTH 
CHANNEL HEIGHT 
NO. OF CHANNELS 
Twg (MID CHANNEL 1081 K 
Twg (MID LAND) 

Twc (CHANNEL BASE) 670 K 


0.071/0.079 (0.028/0.031 IN.) 
0.102 CM (0.04 IN.) 

0.102 CM (0.04 IN.) 

0.257 CM (0.101 IN.) 

390 

(1486 F) 

1066 K (1458 F) 

(745 F) 



WALL THICKNESS 
CHANNEL WIDTH 
LAND WIDTH 
CHANNEL HEIGHT 
NO. OF CHANNELS 
Twg (MID CHANNEL 
Twg 'MID LAND) 

Twc (CHANNEL BASE) 


0.071 CM (0.028 IN.) 
0.102 CM (0.04 IN.) 
0.102 CM (0.04 IN.) 
0.312 CM (0.123 IN.) 
390 

963 K (1273 F) 

984 K (1311 F) 

520 K (475 F) 


INCREASED NUMBER OF CHANNELS 



WALL THICKNESS 

0.064 

CM 

(0.025 IN.) 

CHANNELWIDTH 

0.076 

CM 

(0.03 IN.) 

LAND WIDTH 

0.076 

CM 

(0.03 IN.) 

CHANNEL HEIGHT 

0.257 

CM 

(0.101 IN.) 

NO. OF CHANNELS 

540 



Twg (MID CHANNEL) 977 K 

(1298 

F) 


Twg (MID LAND) 971 K 

(1288 

F) 


Twc (CHANNEL BASE) 616 K 

(649 F) 



WALL THICKNESS 0.071 CM (0.028 IN.) 

CHANNELWIDTH 0.102 CM (0.04 IN.) 

LAND WIDTH 0.102 CM (0.04 IN.) 

CHANNEL HEIGHT 0.257 CM (0.101 IN.) 

NO. OF CHANNELS 390 

100% FPL 135% FPL 

Twg (MID CHANNEL) 848 K (1066 F) 1069 K (1464 F) 

Twg (MID LAND) 847 K (1065 F) 1066 K (1458 F) 

Twc (CHANNEL BASE) 528 K (490 F) 669 K (74“ F) 


(4) SLOTTED 


HOT GAS MID-LAND WITH INCHEASED NUMBER OF CHANNELS 



WALL THICKNESS 
CHANNELWIDTH 
LAND WIDTH 
CHANNEL HEIGHT 
NO. OF CHANNELS 
Twg (MID CHANNEL) 
Twg (MID LAND) 

Twc (CHANNEL BASE) 
MID-LAND SLOT =0.01 


0.071 CM (0.028 IN.) 
0.076 CM (0.03 IN.) 
0.102 CM (0.04 IN.) 

0.3 CM (0.118 IN.) 

460 

1006 K (1351 F) 

1006 K (1350 F) 

643 K (698 F) 

WIDTH X 0.040 DEPTH 


Figure 24. Comparison of Baseline SSME 
MCC Design and Life Enhanced Designs at 
3.1 cm (1.22 Inches) Upstream of the 
Geometric Throat 
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CONTOURED MID CHANNEL WALL DESIGN 


The contoured channel wall represents the minimum change life enhanced 
redesign. The same coolant channel geometry is utilized with the exception of 
increasing the wall thickness (<10%) at midchannel in an attempt to counter- 
act the structural life phenomenon of midchannel thinning. The only design 
change to the current SSME-MCC is contouring the milling cutter to provide a 
slightly convex coolant channel inner wall. This change can be accommodated 
without an engine system rebalance and without any changes to the MCC coolant 
liner or structural jacket mold lines. 

The thermal model and impact on the two-dimensional hot gas wall temperature 
of the MCC coolant channel as a function of wall thickness is shown in 
Fig. 25. The midchannel wall thickness was varied from a nominal value of 
0.071 to 0.079 cm (0.028 to 0.031 inch). A hot gas film coefficient of 135% 
of FPL was used in this study to be consistent with other life enhanced design 
analyses. The results show that the maximum channel wall temperature 
increases 17 K (22 F) when compared to the nominal channel configuration for 
the maximum wall thickness considered. 


3A 


4Twg, 



/ITwg, F 





INCREASED NUMBER OF COOLANT CHANNELS DESIGN 


This life enhanced redesign utilizes a typical approach to increase life by 
reducing hot gas wall temperature while maintaining the same coolant delta 
pressure (AP) and therefore the same pump discharge pressure. This is 
accomplished by reducing the coolant channel width, wall thickness, and land 
width which results in increasing the number of channels from 390 to 540. 
This approach provides a larger surface area exposed to the coolant for the 
same hot gas wall surface area. 

This redesign is the logical life enhanced approach that has been used over 
the past several years for future high pressure MCC design studies. The 
design constraints used for this redesign are shown in Table 1 and compared to 
the current SSME-MCC design parameters. The redesign constraints are typical 
of those used for previous NASA and Air Force studies. 

TABLE 1. COMPARISON OF THE DESIGN PARAMETERS CONSTRAINTS OF THE 
LIFE ENHANCED INCREASED NUMBER OF CHANNELS REDESIGN AND 
THE SSME-MCC DESIGN IN THE LIFE-LIMITED THROAT REGION 


PARAMETERS 

SSME MCC 
DESIGN 

INCREASED NUMBER OF 
CHANNELS REDESIGN 

NUMBER OF CHANNELS 

390 

540 

MINIMUM CHANNEL WIDTH 

0.1016 cm 

0.0762 cm 


(0.040 in.) 

(0.030 in.) 

MINIMUM LAND WIDTH 

0.10922 cm 

0.762 cm 


(0.043 in.) 

(0.030 in.) 

MAXIMUM CHANNEL HEIGHT/ 

5:1 

5:1 

WIDTH RATIO 



MINIMUM HOT GAS WALL 

0.07112 cm 

0.0635 cm 

THICKNESS 

(0.028 in.) 

(0.025 in.) 


This life enhanced design is considered to be currently feasible from a manu- 
facturing viewpoint, but does require some fabrication verification. Areas 
requiring fabrication verification are: 

1. Reduction of the minimum land width from the current 0.10922 cm 
(0.043 inch) on the SSME-MCC to 0.762 cm (0.030 inch). Fabrication 
verification is required to demonstrate adequate exposed surface area 
of the land for bonding the electrodeposited copper strike (to pre- 
vent hydrogen embrittlement) and the electrodeposited nickel closeout 
structure. 
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2 . 


Reduction of the hot gas wall thickness from 0.0.07112 cm (0.028 
inch) to 0.0635 cm (0.025 inch). The bending stress of this redesign 
is more favorable than the current SSME-MCC as reflected by the ratio 
of wall thickness to channel width. However, fabrication verifica- 
tion is required to assure that coining or bulging of the hot gas 
channel wall does not occur from the pressure of the milling cutter. 

The life enhanced 540 channel design axial hot gas wall temperature profile is 
compared to the SSME-MCC design at 135% FPL hg in Fig. 26. The life limited 
throat region temperature distribution is shown in Fig. 26 for the 540 channel 
design and Fig. 27 for the SSME-MCC design at 135% FPL hg. This 540 channel 
design realizes a reduction of 92 K (165 F) in hot gas wall temperature. 

The life enhanced redesign can accommodate the current SSME-MCC or any uprated 
version structural jacket and hot gas wall mold lines without any changes. 
The 540 channel design is optimized by increasing the channel height in the 
combustion zone which allows a more normalized hot gas wall temperature axial 
profile (see Fig. 28) and provides a AP tradeoff to allow the channel height 
reduction at the aerodynamic throat required to meet the SSME-MCC hot gas wall 
mold line. 
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Figure 26. 540 Channel Design Thermal Distribution for 135% FPL hg 
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Figure 27. SSME MCC Baseline 390 Channel Therinal Disbribufcion 
for 135% FPL hg 
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TEMPERATURE 



TEMPERATURE 







KEEL RIB (MID-CHANNEL FIN) DESIGN 


This life enhanced design employs finned channels from the MCC throat plane to 
3 inches downstream of the injector (11 inches upstream of the throat) as 
depicted in Fig. 29 to realize the maximum available reduction in the throat 
region hot-gas wall temperature. 

The finned channel configuration is a straight uniform rectangular shape 
placed at mid-channel. The channel has a bi-width design typical of the SSME- 
MCC as shown in Fig. 30. The channel width is 0.062 inch in the combustion 
zone and narrows down to 0.040 inch in the throat region. The finned channel 
height and fin height are determined by the design constraints shown in 
Table 2. The logic diagram for determining the channel height is shown in 
Fig. 31 for a given fin height to channel ratio. For a constant coolant 
pressure drop, the ratio of wetted perimeter to the flow area is determined by 



where WP, A, and p are the wetted perimeter, flow area, and coolant density, 
respectively. Subscripts 1 and 2 denote finned and nonfinned channel config- 
uration respectively. Since the total heat input for the finned channel is 
usually 1 to 5% more, the density ratio effect is approximated by 0.978. 

The cooling effectiveness of the fin was parametrically defined. Figure 32 
shows the mid-channel and mid-land temperature as a function of fin height to 
channel height ratio. As noted in Fig. 32, the mid-channel is cooled more 
efficiently due to the fin effect. The maximum thermal efficiency is obtained 
for a fin height to channel height ratio of 30% for the SSHE channel geome- 
try. This is the keel-rib design used for life enhancement structural evalua- 
tion as shown in Fig. 33. Typical 2-D temperatures are shown in Fig. 34 for 
the fin height configuration selected in Fig. 34. The life enhanced redesign 
provides a 83 K (150 F) decrease in mid-channel hot-gas wall temperature as 
compared to the current SSME-MCC design. 

COMBINED LIFE ENHANCEMENT DESIGN 

(HOT-GAS WALL MID-LAND SLOTS WITH INCREASED NUMBER OF CHANNELS) 

This life enhanced design consists of cutting axial slots in the liner hot gas 
surface over mid-land. The slots reduce the total effective strain in the 
liner by minimizing the compressive tangential strain. The effectiveness of 
this life enhanced design depends on the impact the slots have on the liner 
temperature and total heat load. A selected slot geometry was analyzed in 
detail and a parametric analysis was made to determine the effects of signifi- 
cant variables. 
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Figure 29. SSME MCC Contour Plot Noting Region of 
Finned Channels (Keel-Rib Design) 
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Figure 30. SSME MCC Coolant Channel Width Profile 
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TABLE 2. SSME MCC FINNED CHANNEL (KEEL-RIB) 
DESIGN CONSTRAINTS 


PARAMETERS 

CONSTRAINTS 

NUMBER OF CHANNELS 

390 (SAME AS SSME) 

CHANNEL WIDTH, BI-WIDTH 

0.157/0.102/0.157 CM 
0.062/0.04/0.062 IN. 
(SAME AS SSME) 

LAND WIDTH, MINIMUM 

0.109 CM 
0.043 IN. 

(SAME AS SSME) 

HOT-GAS WALL THICKNESS 
(INCH) 

0.084/0.066/0.084 CM 
0.035/0.028/0.035 IN. 
(SAME AS SSME) 

MAXIMUM CHANNEL 
HEIGHT/WIDTH 

5:1 





NOMENCLATURE; 

X = MULTIPLIER 
FH = FIN HEIGHT 
C.W = CHANNEL WIDTH 
C.H = CHANNEL HEIGHT 
A = FLOW AREA 
WP = CHANNEL WETTED 
PERIMETER 

SUBSCRIPT; 

1 = NON-FINNED 

2 = FINNED 


Figure 31. Procedure for Determining Finned Channel Parameters 
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SSME MCC Finned Cooling Channel at X = -3.1 cm (-1.2 inch) 
Mainsbage Temperature Profile for a Constant Coolant 
Pressure Drop at 135% FPL hg 
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Figure 33. Keel Rib Geometry and Nodal Structure 
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Figure 34, Keel Rib Thermal Distribution at 135% FPL hg 
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The heat transfer within the hot-gas wall slot is a function of 

1. The radial mass transfer between the combusiton gas and the slot gas 

2. The axial flow of the gas within the slot 

3. The heat flow to the slot walls 

The mass transferred to the slot is from the boundary layer that has developed 
on the chamber liner wall. The boundary layer gas is at a lower temperature 
than the mainstream combustion gas. The transfer of 'this gas into the slot 
represents a complex flow condition which may result in a very thin boundary 
layer at the slot corners, causing high local heat transfer rates. This anal- 
ysis considers only axial gas flow and does not attempt to account for these 
boundary layer corner effects. 

Analytical Approach 

The slot configuration at any axial position was shown schematically in 
Fig. 24. A one-dimensional compressible gas flow analysis is used to deter- 
mine the mass flowrate within the slot. The slot mass velocity is adjusted 
through mass exchange between the slot and undisturbed hot-gas wall static 
pressure. Once the slot flowrate profile is determined, the heat transfer 
between the slot gas and the liner can be determined and, hence, the liner 
temperature can be calculated. The equations used for arriving at the solu- 
tion at each axial position are described below. 

1. Slot Gas Mach Number 

The slot gas Mach number is calculated from the following compressible 
gas flow relationship. 

where 

w = slot gas flow rate 

A = slot cross-sectional area 

Y = ratio of slot gas specific heats 
R = gas constant of slot gas 

P = slot gas static pressure 

Tq = slot gas stagnation temperature 
M = slot gas Mach number 

2. Friction Factor 


The friction factor for the slot gas is calculated from the Darcey- 
Weisbach equation, given below. (This relationship is typical of that 
used for tube and channel flow, ) 



1.74 


1 ,2c 18.7 ■ 
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where 


f = Darcey-Weisbach friction factor 

c = slot surface roughness 

Dh = slot hydraulic diameter 

Re = Reynolds number of slot gas 




3. Heat Transfer 


where 


fSllJws^ determined as 
lahi^n^ slot gas convective heat transfer coefficient is calcu- 
lated using a empirical equation for flow within passages. 


hg = 0.023 


Re®-® Pr®-“ 


hg = convective heat transfer coefficient 
k = thermal conductivity of slot gas 
Pr = Prandtl number of slot gas 

Then the heat transfer per unit mass of slot gas is calculated from 
Q = hg A (T„g - To)/w 

where 

Q = heat transfer per unit mass of slot gas 
A = slot wall heat transfer area 
^wg ” slot wall surface temperature 

4. Enthalpy Charge of Slot Gas 


mass transfer ® Sas depends on whether there is 

tr!Lp^ /u mainstream to the slot (Aw >0) or mass 

ransfer from the slot to the mainstream (Aw >0). The change in 
Slot gas enthalpy per unit mass of slot gas can be written as " 


AH 


C 

P 


tT - T ) , 

o slot 


Aw 

w 
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where 


AH = enthalpy change of slot gas per unit mass of slot gas 
Cp = specific heat of slot gas 
T = Tq of mainstream gas for Aw > 0 

T = Tq of slot gas for Aw > 0 

5. Equation for Change in Stagnation Temperature 

Once the heat transfer (AQ) and enthalpy change (AH) of the slot 
gas are known» an equation for the charge: in stagnation temperature 
can be written 


ATq = (AQ + AH)/Cp 

By substituting for AH and dividing by CpTQ, this equation becomes 



T 

o 


6 . 


C T 
P o 


T ~ T 


o ^ 

w 


Equation for Change in 


Flowrate 


Using the table of influence coefficients for constant specific heat and 
molecular weight, the equation for the change in flow rate as a function of 
pressure gradient, area change, stagnation temperature change, friction fac- 
tor, and Mach number can be written as 


Aw 

w 


f ^1 

-YM^(1 + M^) 

AT 

g 

-yM^ 1 + (y-I)M^ 

Ax yM^ AA 

1 PI ^ 

2 

T ^ 

2 

.» + 

D. , , 2 . 


L 1 - M -1 

o 

2 (l-M ) J 

h (l-M ) A 


2yM^ + V" 


1 - M 


The equations for the change in stagnation temperature and flowrate were 
solved simultaneously for ATq and Aw. By marching axially along the 
slot in small increments (Ax), the flowrate and stagnation temperature pro- 
files for the slot gas flow were determined for a particular slot geometry. 

The stagnation temperature of the free stream gas is assumed to be 6500 R 
throughout the length of the combustor. The liner wall temperature to which 
the slot gas transfers heat is assumed to be 1000 P. For the major portion of 
the liner (injector to throat), this is a reasonable average temperature of 
the surface of the slot periphery (as verified by subsequent two-dimensional 
results) . 
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The free stream static pressure profile corresponds to the local wall opera- 
ting condition. This should be the most severe operating condition relative 
to heat transfer in the slot. The surfaces of the slot are assumed to be at a 
finish of 40 raicroinches rros for calculation of the friction factor. 

At the beginning of the slot, the stagnation temperature of the gas in the 
slot is set equal to the mainstream gas. An initial value is assumed for the 
flow rate in the slot. This is not at all critical since the initial flow 
rate value was found to affect the flow solution for only a short distance at 
which time the solution becomes independent of the initial flow rate estimate. 

Parametric Analysis 

Parametric analyses were made to provide insight to acceptable slot dimen- 
sions. For these analyses, it was assumed the slots were initiated 1 inch 
downstream of the injector, since it is possible the slots will be cut 
throughout the combustion zone. The slots are terminated 2 inches downstream 
of the throat since the liner strains in the nozzle portion of the chamber are 
much lower than the combustion zone and throat. For 1 inch at the beginning 
and end of each slot, the slot depth is tapered linearly to its maximum 
value. Therefore, from 12 inches upstream of the throat (2 inches from the 
injector) to 1 inch downstream of the throat, the slot depth is maintained 
constant. A constant width is used for the entire length of each slot. The 
axial location of the slots in the combustion zone could be tayloreed to the 
amount of strain or strain relief desired to extend the chamber life. Thus, 
while this study investigates slots that extend through most of the combustion 
zone, shorter length slots that are cut only in the throat region may provide 
all of the cycle life extension that is desired. However, the reduction of 
slot adiabatic wall temperature with slot length must be considered. 

The mass flux profile in the slots along the axial position of the chamber is 
shown in Fig. 35. For comparison, the main chamber hot gas mass flux profiel 
is also shown. An increase in mass flux means that gas is being added to the 
slot from the mainstream gas while the opposite is true when the mass flux is 
decreased (supersonic region). Since the heat transfer coefficient is 
directly related to the mass flux, the heat transfer in the slot will be high- 
est where the mass flux is highest (throat region) for a given gas temperature. 

The calculated stagnation temperature profile of the gas in the slots is given 
in Fig. 36. It should be noted that this is an averazge stagnation 
temperature . 

Using the two-dimensional thermal model, the liner wall temperature distribu- 
tion is calculated. Since the flow conditions in the slot are assumed to pos- 
sess one-dimensional effects, the gas temperature is constant around the 
periphery of the slot. Therefore, as shown in Fig. 35, the temperature 
increase over the mid-land is higher than over the mid-channel region. 
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Manufacturing and Structural Life Considerations 


A maximum usable life design would consist of a very narrow slot that would 
basically close up within a few ten-thousandths of an inch during mainstage 
hot-fire testing, negating the aerothermal activity. This would also provide 
the minimum tangential strain. The depth of the slot should be between 1.5 to 
2.0 times the thickness of the channel wall so any resultant geometrically 
oriented strain concentrations at the slot base would be located in a region 
of nearly zero strain. Therefore, the most desirable slot dimensions would be 
approximately 0.102 cm (0.040 inch) deep with a 0.025 cm (0.001 to 0.002 inch) 
width. However, from a manufacturing viewpoint, these slot dimensions are not 
attainable. Manufacturing experts contend that a minimum attainable slot 
width for a 0.102 cm (0.040 inch) depth is 0.025 cm (0.010 inch). 

Selected Mid-Land Slot Design 

The point design selected for this study is an advanced dual-enhanced design 
representing combined thermal and structural life enhancement benefits. The 
slotted mid-land life enhanced design is shown in Fig. 38. It utilizes a min- 
imum channel width of 0.030 inch and a minimum wall-thickness of 0.025 inch, 
identical to the increased number of channels (540) design for enhanced life. 
However, a basic land-width of 0.040 inch is used to accommodate a mid-land 
manufacturing limited minimum slot width of 0.010 inch. This combination 
results in 460 channels. 

The parametric thermal analyses show that a slot with 0.102 cm (0.040 inch) 
depth and 0.023 cm (0.009 inch) width during mainstage would increase the mid- 
channel wall temperature approximately 89 K (160 F) with a mid-land tempera- 
ture increase of over 111 K (200 F) • This temperature level is considered 
intolerable, even with a reduced biaxial strain level. Therefore, the sel- 
ected design assumes that either a physical or fluid barrier exists within the 
mid-land-slot to negate any adverse thermal effects. This can be attained by 
either film cooling at the upstream entrance or bonding a high temperature 
non-structural composite material at the slot base. A physical barrier gap of 
0.00127 cm (0.0005 inch) would allow transverse strain relaxation and negate 
any thermal effects. Adequate cooling of the physical barrier filler material 
would be provided by a braze interface at the slot base. The coolant channel 
geometry to accommodate this design and the 135% FPL hg thermal conditions are 
depicted in Fig. 39. 

This life enhanced design is considered to be an extremely advanced design 
requiring a high degree of technology. It is also combined with the most 
beneficial life enhancement design (increased number of channels design) and 
possesses the majority of the benefits of the 540 channel design. 
Considerable life enhancement can be obtained with this combined design, but 
much advanced technology is required. 
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Figure 39. Slotted Mid-Land Design Thermal Distribution 
for a Thermally Isolated Slot at 135% FPL hg 
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DISCUSSION 


STRUCTURAL ANALYSIS 
HISTORY AND DESIGN ANALYSIS APPROACH 
History of Channel Wall Combustors 

Rocketdyne has designed and tested several channel wall thrust chamber con- 
figurations. Typically, the structural analysis calculates the maximum strain 
range and minimum LCF life to occur in the hot gas wall over the center of the 
land. Hot fire testing of chambers shows a different mode of failure. The 
hot gas wall material at the center of the channel progressively thins from 
the inside of the channel and the material moves and thickens the adjacent 
material (Fig. AO). This phenomenon is a thermal ratchet condition termed 
’’cyclic creep.” In some cases the center of the channel thins sufficiently 
such that the wall fractures with an attendant local loss of coolant through 
the wall. Small cracks have occurred mid-land on the SSME as predicted by 
analysis, but they have not grown large enough to cause a problem. The cyclic 
creep mode of local channel cracking has usually been associated with local 
hot spots from injector effects, and has not been predictable by analysis. 

SSME History and Status 


Main combustion chambers (MCC) for the SSME have been hot-fired nearly 800 
times in static ground tests. The 25 MCCs used have accumulated approximately 
140,000 seconds of hot fire testing. Durations of these tests have ranged 
from less than 1.0 second to 832.0 seconds. Present operating range of the 
SSME is 65 to 1097o of rated power level (RPL), which encompases a combustion 
chamber pressure range of 13.4 MPa (1950 psia) to 22.5 MPa (3270 psia) P^. 

For certain test objec- tives, SSMEs have been tested up to 1117o of RPL, and 

as low as 507o of RPL. 

The current SSME MCC meets all design and operational requirements even though 
local injector effects may cause small areas of channel wall cracking. 
Techniques have been developed to meet those requirements with minor loss in 
performance . 

The initial main injector design had no provisions for film cooling of the 

MCC. After inception of hot-fire development testing, it was discovered that 
local circumferential areas of the combustion chamber had experienced 
over-heating. The overheating resulted in conditions ranging from blanching 
to surface roughening to surface erosion. Surface erosion was generally 
resolved by the addition of film coolant holes and larger face nuts in the 
outer region of the main injector faceplate. These modifications increased 
the fuel flow along the walls of the combustion chamber, resulting in 

reduction of the hot gas wall temperature. In extreme cases, film coolant 
holes are locally enlarged to provide an extra measure of protection for the 
combustion chamber hot gas wall. 
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Note cold wall cracks at corners of channels with hot 
wall through cracks. 


Also, note deformation of EDNi aligned with cracked channels 


P'igure 40 


Two Sections from MCC 9002 Approximately 3.175 x 10~2 
(1.25 in.) Forward of the Throat 



During the early stages of SSME testing, MCC hot gas wall roughening due to 
hot fire operation was not considered detrimental to engine operation. As 
testing progressed, it was discovered that roughening contributed to escalat- 
ing heat loads and further roughening and degradation of the hot gas wall. 
Polishing of the MCC liner between tests with fine grit sandpaper was chosen 
as the method for reducing the deterioration of the hot gas wall. The 
smoother hot gas wall surface is credited with reducing local distress, espec- 
ially on long duration tests. Certification testing of the MCC, up to and 
including full power level (FPL — 109% RPL) , has shown the component to have 
the capability of meeting a certification cycle without deterioration of the 
hot gas wall. This is accomplished with minor maintenance between tests; by 
polishing the hot gas wall if necessary and selectively increasing the film 
coolant flow over local hot spots that indicate signs of duress. Even with 
these modifications, the MCC has encountered minor cracking. This damage 
occurs only in local hot spots and is attributed to injector effects (ano- 
malies) rather than a general deficiency in the cooling capability or design 
of the chamber. 

No performance degradation has been observed in the SSME engine with as many 
as 20 thermal cracks or thermal pin holes in the hot gas wall. These cracks 
tend to close during mainstage when the material is hot and result in little 
or no leakage. Based on extrapolation of observed degradation of the MCC 
liner and resulting performance effects, the MCC will meet minimum performance 
requirements for the required cycle life. 

The degradation of the SSME combustion chambers with operating cycles and time 
is very consistent among components. Cyclic creep damage can be observed from 
initiation to development of a crack in the coolant channel. The consistency 
of data makes the MCC an excellent tool to develop improved analysis proce- 
dures. When this analysis is correlated to laboratory and subscale testing to 
verify improved design concepts for demonstrating the capability for fabricat- 
ing more efficient coolant channels, and as a tool for developing thermal bar- 
riers, the result will be new MCC technology. This technology could be 
applied to improved performance of future engines or uprated SSME 
configurations. 

Existing Analysis Methodology 

Channel wall combustion chambers are sized by hand solutions followed by local 
detailed finite element analysis of a channel cross section. The finite ele- 
ment analysis utilizes a duty cycle considering the prechill, start, mainstage 
and cutoff conditions for temperature , pressure and movement of any support 
jacket. The conditions that maximize the strain range are utilized in defin- 
ing the computerized duty cycle. Detailed temperatures and pressures are 
applied to the model and a maximum strain range calculated. In the 1970 to 
1975 time period when the SSME chamber was designed, the available finite ele- 
ment codes were limited to small deflection theory. Rocketdyne now uses codes 
that can account for geometric deformations as part of the incremental duty 
cycle analysis. 
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The chamber channel hot wall analysis is based on predicted maximum engine 
conditions. The environment is assumed axisymmetric and therefore local 
injector effects are not analyzed from either the heat transfer or structural 
analysis standpoint. It should be noted that for the SSME chamber design, 
temperatures are below the creep initiation temperature so creep relaxation or 
creep rupture are not a significant problem. 

Current Study Effort 

Rocketdyne's current effort in this study has been to analyze the SSME MCC 
using the methodology previously developed under LeRC sponsorship. This 
methodology consists of analyzing the chamber considering more increments in 
the duty cycle and geometric deformation effects. The geometric deformation 
formulation allows the analysis to adjust to the small geometry changes that 
occur during each duty cycle that are cumulative in nature. At higher than 
nominal temperatures, creep relaxation effects are also considered. In order 
to minimize computer time, an extrapolation procedure obtained from LeRC was 
utilized. To perform the analysis, five duty cycles are sequentially analyzed 
and the change in geometric shape is extrapolated to the deformation condi- 
tions 15 cycles further on, e.g., 5 cycles + 15 cycles = 20 cycle condition. 
Using the new geometry, additional duty cycles are analyzed and another extra- 
polated shape projection made. This technique is continued until failure 
occurs. A second part of this study was to evaluate temperature effects and 
design modifications to improve the life of the chamber. 

Once this analysis technique is calibrated to the SSME engine database, it can 
be utilized to analytically determine whether new designs are sensitive to the 
cyclic creep phenomena. These designs may have significant geometric differ- 
ences from the current channel configurations. At this point, the whole anal- 
ysis technique would be only partially validated in that it is calibrated only 
to an SSME MCC and no test demonstration would have been made of its suitabil- 
ity for geometrically different configurations. 
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SSME CONFIGUEiATION AND OPERATION 


The SSME MCC configuration and cross-section shown in Fig. A1 and A2 depict 
the basic geometry and component details. The design consists of an outer 
structural jacket forming the shape of the combustion chamber liner and carry- 
ing the internal pressure and external loads from interfacing components. The 
liner is attached to the jacket at the ends of the structure. It ducts hydro- 
gen coolant along axial channels and acts as a thermal barrier between the 
jacket and the combustion gases, as well as a heat exchanger for heating the 
hydrogen used to drive the low pressure fuel pump. The regeneratively cooled 
liner forms the desired combustion chamber profile. The slotted liner is 
closed out with EDCu and EDNi to form channels for the hydrogen coolant. The 
liner is supported by a high strength structural jacket, and is attached only 
at the ends of the jacket. Structurally, the liner is required to strain out 
to contact the structural jacket, to react the differential pressure load 
between the coolant and the combustion gases, and to accommodate the cyclic 
and racheting strain ranges from the extreme thermal operating environment. 
The structural jacket is required to provide external support for the liner, 
plus react the internal combustion pressure loads and the resulting thrust and 
gimbaling loads generated during hot-fire operation. Even though the liner is 
not attached all along the MCC, the liner motion is restricted during opera- 
tion by the motion of the external jacket. During steady-state operation, the 
liner hot-gas surface is nominally 867 K (1100 F) , and the back wall on the 
jacket side is typically 172 K (-150 F) . During start and cutoff, the com- 
plete liner temperature reaches approximately 33 K (-400 F) . This total cycle 
is a severe thermal cycle for the liner to accomodate. At the analysis loca- 
tion near the throat, the internal pressure of the coolant hydrogen is 43.4 
MPa (6300 psi), which is 30.0 MPa (4200 psi) greater than the hot wall chamber 
pressure of 14.5 MPa (2100 psi). 

Hot-Spot Cracking 

After a series of duty cycles, the geometry of the channel cross-section grad- 
ually distorts, thins along the channel hot wall, and cracks. Figure 40 shows 
typical cross-sections of SSME MCC hardware-through hot spots like the ones 
shown in Fig. 43. 

The geometric changes of the channels that result in cracks are a small por- 
tion of the total surface area of the MCC. The cracked channels are located 
typically in the center of the elliptical-shaped hot spots that occur just 
forward of the throat. The two cross-sections show that the deformation pat- 
tern occurs over a 5 to 10 channel width and each channel has its own geome- 
tric change. The liner consistently shows a thinning of the center of the 
coolant channel’s hot wall, plus a thickening and rounding of the corner form- 
ing the intersection of the channel and land. Detail measurement of channel 
geometries indicate that material volume has essentially been conserved; the 
thinned section is balanced by the thickened shape at the channel corner. 
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Figure 42. Main Combustion Chamber Liner Materials and Part Name 
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Figure 43. Typical SSME MCC Hardware- through Hot Spots 


This means that material has not been significantly removed from the surface 
by the hot gas environment. Other observations are that the channel thinning 
at the point of cracking is 5 to 10 % and that there is only small radial dis- 
tortion at the inner wall. This thinning varies from hot spot to hot spot, 
and liner to liner. Each combustion chamber has similar deformation patterns 
with large variations in the thinning and radial distortions (Ref. 1 and 2). 
In order to adequately analyze this type of structure, the standard structural 
considerations of strength, creep, and low cycle fatigue must be supplemented 
by a cyclic thermal ratchet analysis termed cyclic creep. 

Cyclic Creep Influence 

The cyclic creep phenomenon requires a biased load-structural duty cycle to 
cause an increase to the average maximum hoop strain in the mid channel for 
each duty cycle. Figure 44 summarizes the significant biased conditions that 
occur during the SSME duty cycle. The observed thinning and distortion of the 
channels show that the mid-channel thinning is a result of a local biased 
negative radial strain and a compatible increase in a local hoop strain. The 
reshaping of the corner is dominantly a large shear strain which resolves into 
plus and minus principal strains that thicken the corner zone. 

The dominant load is thermally induced, but the other biased conditions accen- 
tuate the cyclic creep. For the SSHE HCC, significant biased conditions 
include pressure, creep or relaxation effects, material strength effects, and 
the liner support. The differential pressure loading from the coolant varies 
throughout the duty cycle, and is significant only over a portion of the duty 
cycle. This variation in pressure may be part of the difference between the 
small radial distortion mid-channel and the thinning rate for this HCC vs 
results from smaller test chambers at Rocketdyne and NASA that maintain signi- 
ficant differential pressures throughout the test duty cycle. 

The maximum average liner channel wall temperature in a hot spot is approxi- 
mately 867 K (1100 F) . The maximum liner temperature of 1089 K (1500 F) in a 
local hot spot is well above the creep activation temperature for NARloy-Z. 
However, the temperature gradient through the hot wall is approximately 672 K 
(750 F) , which results in half of the hot wall being at temperatures below the 
creep activation temperature. The primary load-differential coolant pres- 
sure in the HCC can be easily supported by this low temperature portion of the 
wall. This essentially changes the load condition from a potential load- 
dependent creep problem to a local relaxation phenomenon. The stress state on 
the hotter portion of the wall reduces to a displacement-limited steady-state 
condition. The primary load is essentially shifted to the cooler half of the 
channel wall. Therefore, the creep/relaxation effects are a minor addition to 
the biased loading in this hardware, whereas it could be a more significant 
parameter if the total hot wall was above the creep activation temperature. 
This can occur in a plugged channel versus a hot-spot condition. 
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Figure 44. What's Needed to Demonstrate Cyclic Creep 





Plugged channels crack earlier than a typical hot-spot condition. Figure 45 
shows creep relaxation results for NARloy-Z from standard materials test 
specimens. It is evident that relaxation effects occur rapidly within less 
than 1 second of hold time and achieves an approximate steady-state condition 
within 10 to 20 seconds. For the SSME engine duty cycle, this means that 
steady-state creep relaxation occurs during every typical test firing of 50 to 
800 seconds. The creep effects are basically test-cycle dependent, not 
test-time dependent. 

Another parameter for the biased loading is the temperature-dependent cyclic 
strength of the NARloy-Z material. The hot-cold duty cycle causes a bias in 
the strain cycle because the elastic vs plastic strain distribution is differ- 
ent in the overall channel cross-section between steady state and the remain- 
ing portion of the duty cycle. 

The cyclic stability of the material can also be a significant parameter for 
some of the materials used for coolant liners. In this case, the NARloy-Z 
material used for the MCC is cyclic stable, i.e., it has a relatively small 
increase in load at constant cyclic strain range (Fig. 46). Therefore, cyclic 
strain hardening effects are not a significant consideration, and uniaxial 
stress-strain data were used in the analysis. 

The design concepts used for the MCC also bias the cyclic strain. The liner 
and jacket are similar in operation to a tube in a tire where the outer jacket 
limits the motion both radially and axially during the pressurized portion of 
the duty cycle. During the portion of the duty cycle where the liner is cold, 
the liner separates and pulls away from the jacket, which reduces the strain 
range, but biases the strain cycle. 

The other primary contributor to the cyclic creep mode is the geometric dis- 
tortion of the liner cross-sections. Previous SSME MCC analyses performed in 
1970 to 1975 were limited to small displacement analyses that utilize only the 
original geometry. This means that geometry distortions are not accounted for 
in developing the cyclic straining of the liner. The maximum strain or ten- 
sile instability-type failure mode typically requires the displacements and 
distortions of the structure to be considered as part of the stiffness-load- 
stress and strain conditions. The current analysis methodology incorporates 
those geometric changes that occur throughout the operating duty cycle. Anal- 
ysis results show that this is a significant parameter in changing the rate of 
distortion as the analysis progresses through the total duty-cycle analysis. 

The life analysis, therefore, consists of evaluation of the liner through a 
series of duty cycles to a point of liner cracking. From a low cycle failure 
evaluation standpoint, this reduces to a low cycle fatigue (LCF) failure mode 
that is related to the cyclic strain range vs a cyclic creep failure mode that 
is related to the total strain. The LCF is evaluated from an individual 
assessment of the cyclic strain range of local zones of the liner channel 
cross-section. The cyclic creep failure mode is evaluated as an average 
strain condition in the channel hot wall that is conveniently related to the 
channel wall thinning in the analysis methodology used. 
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Figure 45. NARloy-Z Thermal Fatigue Test - Load vs Cycles 
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Duty Cycle Analysis Procedure 


The duty-cycle analysis may require tens or hundreds of cycles prior to a 
failure condition. In cases where a stablized condition occurs in one or two 
cycles, the analysis can be terminated at this point, and the failure or 
cracking condition accurately extrapolated. For the structural response of 
the MCC, this is not possible, since the cyclic strain condition does not 
stabilize; and the gradual geometric distortions are significant in changing 
the rate of strain cycling and maximum strain. It is also impractical to con- 
sider analyzing every duty cycle. A compromise is to utilize the approximate 
methodology developed under NASA-funded efforts (Ref. 3 and 4), that analyze a 
series of duty cycles to determine the conditions consistent within a particu- 
lar portion of the distortion sequence, extrapolates the geometry to a geome- 
tric condition further into the deformation pattern, and then analyzes another 
series of duty cycles to determine local conditions at this point. It should 
be noted that this is a specialized analysis technique, not a general methodo- 
logy applicable to a large class of problems. The extrapolation technique 
requires that the geometric distortions further into the life cycle can be 
linearly approximated from a sequence of previous distortions. The analysis 
technique requires a balance of the cost in performing sufficient duty cycles 
to determine a local, stablized condition and the extent of the extrapolation 
further into the duty cycle where the rate of cyclic change is not signifi- 
cantly different. This study utilized a standardized approach of developing a 
5 duty-cycle detailed analysis followed by a 15 duty-cycle extrapolation of 
the geometric shape. The model was then rerun based on the distorted geometry 
for the next sequence of cycles. This procedure develops the distorted geome- 
try shape for the analysis, but does not track the strain history through the 
total history. The strain history is tracked only within each of the individ- 
ual 5 duty-cycle analyses. 
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STUDY METHODOLOGY 


Overview 


The SSME MCC is a three-dimensional shell structure. The detailed analysis of 
the MCC liner was performed using two separate, two-dimensional models one for 
the jacket boundary conditions and the other for the actual liner analysis. 
Figure 47 outlines the overall procedure used in this study for structural 
analyses. This effort consisted of an initial time-phased duty-cycle analysis 
for both the jacket and liner. Since the jacket loading and displacements 
were insensitive to the variations in thermal conditions or geometry changes 
in the liner, only one analysis is required to develop the boundary conditions 
for the liner model. The coolant liner loading and model varied as required 
in each task of this study. 

A summary of the analysis basis is furnished in Fig. 48. The two models and 
duty cycle are outlined, as well as the structural parameters that were con- 
sidered in the analysis. Loadings include temperatures, pressures, jacket 
restraints, and the initial fabrication gap between the liner and jacket. The 
models were an axisymmetric jacket model and a generalized plane strain cross- 
section model of the liner. The liner model analysis included kinematic 
hardening creep material characteristics and geometric distortion effects. 
The basic analysis of the liner for each configuration was a series of 5 duty- 
cycle analyses, followed by a 15-cycle extrapolation. 

Figure 49 is a block diagram of the liner analysis procedure for each geometry 
or thermal condition analyzed. The loads for each model or condition were 
defined and applied to the basic liner model. The 5 duty-cycle analysis was 
complete, results reviewed, and life predictions made as required. Next, the 
geometry was updated to represent the extrapolated shape at a point 15 cycles 
further into the analysis. This information was then combined with the basic 
input for the initial finite-element model, and the analysis repeated. 

Duty Cycle Definition 

At full power level (FPL) conditions, the 100% hg case represents the normal 
design conditions for the liner, and is representative of tempertures in the 
liner in over 95% of its surface area. In the local hot-spot regions, the 
temperatures are significantly hotter, and are represented by the 135 percent 
hg thermal analysis case. The 135% hg condition was chosen as the baseline 
duty cycle for establishing the analysis technique and verifying that the 
overall methodology would demonstrate the cyclic creep phenomenon. The duty 
cycle was defined and evualated, and divided into 11 discrete time slices to 
establish the analytical representation of the cycle. Figure 50 is a plot of 
hot-gas wall temperature, jacket temperature, Pfuel ^HGW time, and 
shows graphically the time history of the significant events for the start/ 
steady-state/ shutdown cycle, as well as the 11 discrete time-slice points. 
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Figure 47, Overall Structural Analysis 
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Figure 48. Analysis Approach 




Figure ^9. Liner Analysis Block Diagram 




Twg, TjackET K 


' JACKET 





Jacket Model and Duty Cycle Motion 


The MCC jacket model was run to establish radial displacement and axial strain 
influence coefficients for both combustion chamber pressure and average jacket 
temperature. Using heat transfer data generated from the 100% hg case FPL 
duty cycle run, compatible time slice values of jacket temperatures and com- 
bustion chamber pressure were established, which in turn were used to define 
the compatible radial displacements and axial strains of the jacket used for 
the liner cross-section model boundary conditions. 

A verbal description of the sequential events helps to understand the time- 
temperature-pressure relationships between the liner and jacket. The MCC as 
manufactured can have a maximum radial gap between the liner OD and structural 
jacket ID of 0.5 mm (0.020 inch). During the start prechill, the cold hydro- 
gen gas flowing through the liner coolant channels causes a thermal contrac- 
tion of the liner which increases the liner/jacket radial gap. During this 
time, there is a negligible cooling effect on the jacket due to the cold radi- 
ation from the liner. Then, upon ignition, the combustion pressure increases 
and the weak liner expands radially, like an inner tube in a tire, at a fairly 
low pressure until it contacts the jacket inner diameter. Once liner/jacket 
contact is made, the jacket mean temperature experiences a fairly rapid drop 
in temperature to approximately 200 K (-100 F> , while combustion pressure is 
increasing to its steady-state value. During this time, there is a net radial 
growth of the jacket from the positive radial displacement of the jacket due 
to combustion pressure. This displacement is partially offset by the negative 
radial motion of the jacket as it drops in temperature* After steady— state 
combustion pressure is reached at 4.9 seconds, the jacket continues to get 
colder over the next 15 seconds, resulting in a decrease in the jacket radial 
displacement and axial strain. 

At this point, pressure and thermal equilibrium loads in the jacket have been 
achieved. During shutdown, combustion pressure is throttled down to 65% of 
rated power level, and then the engine is shut down. During shutdown, the 
purge flow of hydrogen through the coolant channels rapidly chills the liner 
to 33 K (-400 F) , and the liner separates from the jacket. When combustion 
pressure reaches zero, the mean temperature of the jacket has warmed up to 
approximately 211 K (-80 F) . The jacket eventually warms back up to room 
temperature, and returns to essentially a line-to-line contact position with 
the jacket. 

Liner Cross-section Model 

The local generalized plane strain model for the current SSME configuration is 
shown in Fig. 51. This model has gap-type (compression only) elements at the 
liner-to-jacket interface. Only radial constraint that varies with the jacket 
displacements was considered at this interface. 
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This model is consistent with previous SSME MCC models and previous work per- 
formed for NASA-LeRC on cyclic creep analysis. All models analyzed used the 
same symmetry and boundary conditions. As previously stated, the jacket is 
the dominant structure, and the liner must follow its movements when the two 
are in contact. This model and loading technique is a good representation of 
the major portion of the liner that is unaffected by the local hot spots. 
When analyzing for the hot-spot condition (135% hg thermal profile), it is 
assumed that the hot spot occurs each duty cycle at the same position on the 
liner, the hot-spot covers a large enough zone, and other restraints are such 
that the jacket symmetry and boundary conditions are maintained in the center 
of the hot spot, as will be discussed later. Both analysis conditions were 
used for demonstrating the cyclic creep phenomenon and the 135% hg case in 
ranking the potential design improvements. 

Material Properties 

The four materials that make up the liner cross-section, INCO 718, EDNi, EDCu, 
and NARloy-Z; were represented in the structural model with temperature- 
dependent, elastic and inelastic properties (see the Materials Section for 
detailed materials data). The nonlinear behavior was represented with a 
bilinear stress-strain curve and the analysis used kinematic strain harden- 
ing. Since NARloy-Z in the condition used for the MCC is cyclically stable, 
i.e., it has a relatively small increase in load at constant cyclic strain 
range, strain hardening was not part of the analysis. Monotonic tensile 
stress-strain curves were used in developing the bilinear curves. Creep/ 
relaxation constants consistent with available data were used in the analysis. 

Large Displacement (Updated Geometry) 

A significant change from previous SSME MCC analysis is the inclusion of geo- 
metric updating of the structural model. This is accomplished analytically 
within the computer code analysis for each increment of the duty cycle anal- 
yzed, and through an extrapolation code which defines a distorted shape for 
restarting the finite-element analysis 15 duty cycles further into the 
analysis . 

The geometric distortions are accumulated throughout the analysis and total 
raid-channel thinning was used as a measure of the total (maximum) strain 
condition. 

The finite element analyses were performed with the APSAC code (Finite Element 
Axisymmetric and Planar Structural Analysis with Loading and Creep Duty 
Cycles, Ref. 5). The APSAC code was used for the jacket analysis where an 
adequate operational model was in existence and was rerun to develop the 
detailed radial and axial displacements for the incremental loading cycle. 
APSAC was also used for the liner hot-gas wall analysis. 
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APSAC is an in-house developed program for finite element analysis. Tempera- 
ture-dependent material properties, inelastic analysis using bilinear repre- 
sentation of stress-strain curves, and either kinematic or isotropic strain- 
hardening theory are incorporated in the program. 

Large displacement (updated geometry) analysis is available as an option. In 
addition to the quadrilateral axisymmetric or planar element, a compression- 
only gap element is available for uses such as the liner/jacket interaction. 
A secondary-type creep law, ey = Ao®, is implemented for creep or relaxa- 
tion effects. An efficient mesh generation scheme is available as a portion 
of the basic program. Grid and displacement plots, contour plots, and graph- 
ical plots are available through a postprocessing program. Information for 
geometric updates or geometry extrapolation can be performed using data-sets 
generated by the program. 

Rocketdyne has used APSAC and earlier versions of the program since 1969. 
This program was the primary analysis tool used for the many inelastic analy- 
ses performed on the various SSME components. The original MCC liner analysis 
used the APSAC code prior to implementation of the large-displacement option. 
Other Rockwell divisions use the program, and it is implemented and used by 
NASA-MSFC for SSME hardware analysis and other in-house analyses. APSAC has 
been verified against classical analyses; other codes such as ANSYS and a 
variety of specimen and test hardware with strain-gage verification, as well 
as production hardware usage experience. 

Displacement Extrapolation 

To reduce both engineering and computer costs, a displacement extrapolation 
program was used. A Rocketdyne code was developed for accomplishing this 
task, and is discussed in detail in Appendix C. The extrapolation technique 
used was a standard linear least-squares approximation for the computed rate 
of change of displacement. 

Analysis Strategy 

The completion of this study required a minimum of 36 runs consisting of a 
5-duty-cycle run plus an extrapolation run and many plot runs. In the course 
of the work, approximately 50 runs were made. Each 5-duty-cycle analysis con- 
sisted of 55 load increments with different temperatures, pressures, and 
boundary conditions applied in each increment. The temperatures and pressures 
were obtained from separate heat transfer studies using the REGEN and DEAR 
analysis codes. Considerable planning and strategy were required to develop 
techniques to track each run and increment, simplify input data required by 
the user, and minimize costs in manhours, computing, ^ printout, and graphical 
output. The intent was to develop techniques usable not only for this study, 
but applicable as a standard analysis methodology for advance design studies, 
as well as detail design variations and analysis. Details of these techniques 
are included in Appendix D as part of the example input data for the 135% hg 
case baseline analysis. 
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Creep Fatigue Analysis - Life Prediction 


General MethodoloRV . The method of analysis used included a summation of 
incremental damage caused by fatigue and creep mechanisms at progressive 
stages of thrust chamber operation (life). Fatigue is the dominant component 
in this study, but the general creep fatigue methodology is presented, since 
its proportion may vary from one MCC configuration to another. 

Creep fatigue damage is determined from the stress-strain-time history of the 
structure. Low cycle fatigue damage, and creep rupture damage, 
0c » are considered. The total damage, 0^, is represented in equation 
form as 


0t = A 0f + B 0c 

where A and B are damage interaction coefficients. 


Predicted life is found as the reciprocal of total damage. The method of 
analysis is shown schematically in Fig. 52. This is the manner originally 
proposed by Robinson (Ref. 6) and followed by Taira (Ref. 7), Hanson and 
Halford (Ref. 8), and Spera (Ref. 9). A similar procedure has been finalized 
and incorporated into the ASHE Pressure Vessel Code for creep-fatigue evalua- 
tion of nuclear power components. Code Case 1592, Appendix T. 

Channel Wall Fatigue Analysis . The fatigue analysis of the channel wall is 
performed using peak local strains and temperatures. The cyclic range and 
mean value of the strain were determined at progressive stages of thrust cham- 
ber life. The results from the duty-cycle load increments for the most cri- 
tically strained element are used in a formulation (Ref. 10) analogous to the 
von Hises criterion to find the strain range, c^, and the mean strain, 
®m" 


The values of the strain range components, Ac^, are calculated from the 
differences in each of the strain components, cj, for c^, and the sum 
of each of the strain components, cj, for Cj„ are the strain components 
at the duty cycle of interest. The values for Acj are determined from the 
maximum and minimum strain conditions in the cyclic analysis, as was previ- 
ously mentioned. The initial duty cycles adjust the computer model for the 

inelastic stress and strain state. The formulations for c*. and c„ are 
as follows: 



Acg)^ + (Ac^ - Acj)^ + 


(Acg - Ac^)^ + 3/2 (Ay^ 3 + Ay^^ + Ay^g) 
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Von Hises effective strain defines the uniaxial equivalent for the multiaxial 
strain state. 
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Creep Rupture Analysis . Creep damage for the total duration is determined by 
the summation of the ratios of time spent at a particular stress to the time 
to rupture at that stress. 


An illustration of the summation process is given in Fig. 52. If the tempera- 
ture varies simultaneously with the stress, the time-to-rupture curve for the 
specific temperature is used for each stress increment. The time-to-rupture 
data that are used are obtained from monotonic load-based creep rupture tests. 

As long as the creep rupture damage component is kept small, little inaccuracy 
is introduced using monotonic creep rupture data. Constants for the primary 
and secondary creep laws were used in the finite element solution to calculate 
wall creep deformation. 

Static Strength Analysis . The static strength of the thrust chamber channel 
walls was reduced by the thinning and bulging during the cyclic process. The 
wall-thickness adequacy to withstand tensile loads occurring during the duty 
cycle was verified as an additional check of life capability. 

Strain Limit 


The maximum strain capability of the NARloy-Z material is approached from the 
tensile instability failure criterion. The maximum hoop strain is related to 
the radial strain and the thinning of the channel wall. Tensile instability 
is defined in this case as where the incremental strain hardening of the 
deforming metal of the wall is less than the incremental increase of true 
stress due to local wall thinning. Once the instability condition is reached, 
further increase in load results in rapid failure of the channel wall. The 
criterion defines a critical strain limit for the material. The analysis 
method for determining the strain at the onset of instability (called critical 
strain) is based on the work of Svensson (Ref. 11). The MCC section between 
lands is modeled as a thin shell plate that has an axial strain of zero during 
the tensile instability portion of the duty cycle. This strain condition 
occurs during inelastic straining in a 2:1 (hoop to axial) stress field, which 
is consistent with tensile instability of a pressurized thin cylinder. This 
stress state is only a first order approximation of the true-stress state, 
since the actual hardware has a more complex stress ratio due to the overall 
duty-cycle response to thermal load variations. The maximum tensile stress 
state occurs during the cutoff portion of the duty cycle as the temperatures 
start to fall off and there is still pressure in the channel. 

The section of the hot-gas wall between the lands is shown in Fig. 53 with the 
stated stresses and strain condition. 
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Figure 53. Critical Conditions at Tensile Instability 


The critical tensile strain in the hoop direction at instability is equal to 
the material strain hardening coefficient, n, divided by 2. Once the critical 
tensile strain in the hot-gas wall is known, the corresponding change in wall 
thickness, or thinning, is determined from the radial strain, Cj.. 

The sum of the inelastic strain components is zero, based on a constant volume 
criterion. Since the axial strain = 0, and the tensile strain in the 

hoop direction cjj = n/2, at instability the radial strain is found by 

solving the following: 

cr + ch + ca = 0 

Cj. + n/2 + 0 = 0 
Gj. = -n/2 

n 
2 

or 



-n 


t 

cr 



2 

e 


where 

tgp = thickness of hot-gas wall for tensile instability 
tg = original wall thickness 
t 

cn 

In ~ — = definition of strain in the radial direction 
^o 
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Using this tensile instability methodology, the critical strain or thinning at 
failure can be calculated knowing the value of n, the exponent in Ludwig's 
Law. When stress-strain data are not available, a reasonable approximation 

for n can be determined using the tensile yield and ultimate values. Using 
the equation 


n 


.2 



Figure 54 plots the value of n versus temperature and thinning at instability 
for NARloy-Z, based on Boeing and Rocketdyne strength data plus the analysis 
criteria developed above. 


For the MCC, the maximum tensile stress occurs in the 590 to 700 K (600 to 
800 F) temperature range and generates a range of thickness reduction of 5 to 
8% at tensile instability. 

A review of typical cross-sections from hot-fired MCCs show that failure can 
occur at an estimated 5 to 10% thinning of the channel wall. There is a wide 
variation in this thinning for different MCC chambers data base, but a 5 to 8% 
thinning limit is a reasonable criterion for the minimum strain at failure for 
the SSHE. The variation in test results may be attributed to conditions in 
the heat input at hot spots, flow variation in a specific combustor and the 
sensitivity of instability analysis to small changes in material orthotropic 
properties. The 5 to 8% strain limit compared to the uniaxial strain limit of 
50% is consistent with a literature survey study Rocketdyne has previously 
made of available biaxial test results vs uniaxial test results. 


Ultimate elongations of 12 metals were obtained at room temperature in the 1:0 
(hoop/axial) , 1:1, and 2:1 tensile modes (Ref. 2 and 3). The specimens were 
made from steel, aluminum, and titanium sheet. Plots of the results (Fig. 55, 
56, and 57) indicate a large decrease in available uniaxial elongation due to 
the biaxial mode even for such ductile metals as 301 stainless steel. 
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ANALYSIS RESULTS 


As previously mentioned, the 135% hg thermal load case was used as a baseline, so 
that cyclic creep effects would be a sensitive influence in the structural 
analysis. Also, it more closely reflects conditions in a hot spot where known 
cracking occurs. The 100% hg condition — away from hot spots on the MCC — has not 
shown the channel cracking mode to date. Not knowing what to anticipate from the 
computer analysis, and desiring maximum sensitivity for evaluating the various 
recommended life enhancement designs, the hotter condition (135% hg) was the more 
logical choice for the baseline analysis, especially since cyclic creep was the 
primary interest of this study. 

To date, there is no detailed textbook explanation of how cyclic creep takes 
place or the relative importance of individually biased loads or geometric shape 
on cyclic creep. A simplistic sensitivity study was run on various geo- metric 
and property variations to try to better understand individual vari- ables. 
Whereas this technique was quite successful in the thermal analysis, it was 
inconclusive in the structural analysis. The structural boundary con- ditions 
from the operational duty cycle forces the entire hot wall of the liner well into 
the inelastic regime. This seems to require accurate shape changes or property 
variations that are truly representative of the geometric distortions and 
materials to furnish usable results; i.e., only computerized analysis utilizing 
updated geometries produced acceptable information for judging life enhancement. 

So, to obtain an overall perspective of cyclic creep, the basleline-analysis 
results have been studied in several ways. One class of studies was from an 
overall standpoint that is related to observable hardware-channel distortions and 
cracks. These results are: 

1. Overall distorted shape versus test cycle 

2. Thinning and cracking versus accumulated duty cycles 

3. Thinning versus radial distortion difference between mid- 
channel and midland 

The other class of studies related to trying to understand calculated stress- 
and-strain conditions and distortions from the standpoint of how to improve a 
design or evaluate a configuration before a detailed analysis. This effort 
required a detailed study of start, steady state, cutoff, and posttest incre- 
ments within a 5 duty-cycle analysis. Areas investigated included: 

1. Stress and strain distribution 

2. Critical strain locations for LCF analysis 

3. Distorted shape throughout the duty cycle 

4. Incremental stress-strain hysteresis loops for a critical 
point 
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5. Simplified correlation of biaxiality-strain distortion 

6. Definition of maximum gross tensile-stress condition in the 
duty cycle when cracking actually occurs 

The baseline analysis will first be discussed completely, followed by correla- 
tions with the remaining studies. 

Distorted Shape 

The shape change with accumulated duty cycles gives a perspective of how well 
the analysis procedure predicts the type of distortions observed from cross- 
sections of the SSME MCC. Figure 58 is a composite of the finite-element 
model distorted shape at the end of each 5 duty-cycle analysis (5, 25, 45, 65, 
85, 105). A gradual change in the cross-section is noted as the cycling pro- 
gresses. The midchannel point thins and the initial sharp corner at the land 
to channel gradually thickens and smooths. There is only a small permanent 
set inward (toward the centerline) at the midchannel point relative to the 
midland position. These geometric changes are typical of what has been 
observed on the SSME MCC cross-sections in the regions of local hot spots. 

Thinning and Cracking Versus Duty Cycles 


The maximum channel thinning versus accumulated duty cycles was used as the 
measure of the total average strain on the channel wall. Because of the 
approximation technique (updating the geometry, but losing the strain state) 
used for determining the restart condition after a 15 cycle extrapolation, the 
total strain in each element is not directly available. However, it is a 
rather simple task to track the accumulated thinning of the channel wall. The 
geometric updating is all part of the extrapolation code discussed in Appendix 
C). There was still additional judgments required to utilize this informa- 
tion. Figure 59 depicts the baseline 135 hg analysis results and shows the 
shakedown condition that occurred in each 5 duty cycles, and the method of 
adjusting the raw analysis results to account for the shakedown. Each config- 
uration and operating condition analyzed had similar shakedown conditions. 
For the first 5 duty cycles, the adjustment of the residual stress and strain 
conditions during the load cycling resulted in a realistic shakedown condi- 
tion. The first shakedown condition also included the gross distortion of the 
liner as it permanently yields through the 0.51 mm (0.020 inch) tolerance gap 
and touches the jacket. Thereafter, at the steady-state posttest condition, 
the liner remains virtually in line— to— line contact with the jacket, so the 
first shakedown is larger than the remaining ones. The other shakedown condi- 
tion, where there is a rapid change in thickness for the first 1 or 2 cycles, 
followed by a much slower change in cycles 3, 4, and 5, occur because of the 
approximations made. Each new 5 duty cycles start with a stress-and-strain- 
free distorted geometry rather than a shakedown continuation of the previous 
cycling. If the raw analysis results were used directly, a significant error 
in thinning would be obtained. 
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Figure 59. SSME MC Life Predict 


ANALYSIS METHODOLOGY 

5 DUTY CYCLE + 15 CYCLE EXTRAPOLATION 
RESTART WITH 

• DISTORTED GEOMETRY 

• VIRGIN MATERIAL(NO RESIDUAL STRAINS) 

• MCC LINER IN CONTACT WITH JACKET 

INITIAL "SHAKEDOWN" ONLY REALISTIC FOR 1-5 
DUTY CYCLES 

"SHAKEDOWN" EFFECT SUBTRACTED OUT 


- Baseline 135% hg 




The final methodology chosen to eliminate the shakedown effect was to assume 
that the distortions from the first cycles in the 5-cycle analysis were equiv- 
alent to a series of thinning cycles at the rate of change used for extrapola- 
tion to the restart condition. This is effectively obtained by extending the 
last 5-20-cycle extrapolation line (e.g., A) until it intersects with a hori- 
zontal line drawn from the end of the shakedown condition, B. At this point, 
the B slope applies until the C slope horizontal line intersects. This tech- 
nique was used to adjust all the raw analysis results in this study. As one 
would expect, it has the largest effect on slowly thinning analyses with 
larger shakedown conditions. The raw analysis results for all analyses are 
located in Appendix G. 

The adjusted results for the baseline configuration 135% hg and 100% hg and 
the strain limitations are shown in Fig. 60. The plot is developed as a per- 
cent of thinning versus cycles for two reasons. One is that various channel- 
wall thicknesses are used in the analyses, and a common baseline is desirable; 
and the other is the percent thinning is equal to the accumulated radial 
strain, and is within 10% of the accumulated effective strain. The results 
show, as expected, significant differences in thinning for the two conditions 
1600% hg and 135% hg) . The hot-spot condition (135% hg) failure life is 
between 80 and 125 duty cycles. SSME test history puts this at somewhere 
between 5 and 30, so there is still a significant difference between the 
cyclic-creep analysis and hot-fire test data. The most likely reason for this 
discrepancy is the approximation required to model the hot-spot condition as a 
plane strain model. The analysis assumes the hot spot occurs in a total cir- 
cumferential band around the MCC. In reality, the hot spot is a local condi- 
tion over a range of 5 to 10 channels. Appendix discusses in detail a 
refined approach to more accurately analyze a local hot spot and presents the 
results of a preliminary study. From a heat transfer standpoint, however, the 
size of a hot spot is negligible, as long as it covers several channels and is 
reasonably long, as has been experimentally evidenced. 

From a structural standpoint, relative to the remainder of the liner, the hot 
spot tends to act like a solid material for the start and cutoff conditions 
(when it is cool), and an elliptical hole filled with soft material (when it 
is not) during the steady-state firing conditions. This condition can signif- 
icantly increase both the out-of-plane axial strain condition and the hoop 
strain. The detail modeling of a hot spot requires a complex 3D inelastic 
analysis that is well beyond the scope of this study. The basic correlations 
and analysis, though, of this study’s methodology are useful for correlation 
between models and temperature changes, but are not considered to be exact 
answers. The 100% hg condition relates to a circumferentially uniform temper- 
ature condition, and the planar model is a more accurate representation. The 
calculated life of 290 to 445 duty cycles, for the lower temperature nominal 
SSME-MCC condition has not yet been verified by test data in that individual 
MCCs have not been tested more than 80 to 100 duty cycles. 
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DUTY CYCLES 

Figure 60. Midchannel Permanent Deformation (Thinning) vs 
Number of Duty Cycles 



Detail Duty Cycle Review 


This review presents detailed information and observations gleaned from an 
intense study of the 135% hg baseline analysis. Its intent is to furnish 
additional information to aid analysts in understanding possibly how cyclic 
creep occurs, how to improve a design, strain distributions, and how a typical 
S-cycle hysteresis loop looks. These results were specifically used in devel- 
oping the contoured life-enhancement design. 

For the review, the hot wall was separated into three zones: midland (ML), 
channel to land (LC), and midchannel (MC) . Figure 61 defines the characteris- 
tics, from a structural analogy, for the three zones. Figure 62 and Table 3 
describe the zonal thinning in a relative manner and list the thickness 
changes for the different key load incremenets per cycle. Figure 63 is a com- 
posite of critical time-slice distorted shapes. 


TABLE 3. SUMMARY OF CHANNEL WALL THICKNESS PER CYCLE 
FOR THREE KEY LOAD INCREMENTS (135% hg) 




CHANGE IN THICKNESS,* 

LOAD INCREMENT 

LOCATION 

mm 

INCH 

INCREMENT 6, 

ML 

0.0196 

+7.7 X 10-4 

STEADY STATE 

LC 

0.0218 

+8.6 X 10-4 


MC 

0.0183 

+7.2 X 10-4 

INCREMENT 10, 

ML 

-0.0053 

-2.1 X 10-4 

POST CHILL 

LC 

-0.0053 

-2.1 X 10-4 


MC 

-0.0091 

-3.6 X 10-4 

INCREMENT 11, 

ML 

-0.0017 

-0.67 X 10-4 

ROOM TEMPERATURE 

LC 

-0.0017 

-0.68 X 10-4 


MC 

-0.0053 

-2.1 X 10-4 

ML » MIDLAND 

LC = LAND/CHANNEL TRANSITION 
MC = HIDCHANNEL 

♦CHANGE IN THICKNESS FROM ORIGINAL GEOMETRY 
(ABSOLUTE) = AR(i)-AR(2) 


Table 3 shows that during steady state (increment 6), even though the hot wall 
is getting thicker due to aAt and compressive hoop load, the MC is getting 
less thick than ML or LC, or relatively speaking, MC becomes thinner than ML 
and LC. 
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DESIRABLE TO SEPARATE INTO 3 ZONES 


ML. LC, MC - MID-LAND, CHANNEL TO LAND & MID-CHANNEL 
EACH ZONE HAS DIFFERENT CHARACTERISTICS 

ML - 3D ZONE OF MATERIAL WITH FREE SURFACE (BIAXIAL RESTRAINT) 
(ANALOGY - THICK CYLINDER INSIDE SURFACE) 

MC - 2D ZONE OF MATERIAL WITH 2 FREE SURFACES (LOWEST RESTRAINT) 

(ANALOGY -THICK PLATE/BEAM WITH ZERO STRESS THROUGH THICKNESS) 

LC - NOTCH ZONE WHERE 2D TO 3D EFFECTS TRANSITION (HIGH SHEAR) 
(ANALOGY - (T-CONCENTRATION IN STEPPED BAR OR CYLINDER) 


HOT 

WALL STEADY STATE 

DEFORMED/UNDEFORMED PLOT 



Figure 61 


135% hg MCC Analysis Detail Review, Zone Definition 






© THINNING/RELATIVE THINNING/CYCLIC CREEP 

o relative THINNING APPARENT THROUGHOUT CYCLE 

• RATIO OF THICKNESSES AT ML, LC AND MC ZONES 

o DURING DUTY CYCLE 

• ABSOLUTE THICKNESSES INCREASE UNDER STEADY STATE 

. COMPRESSION LOADING 

• a AT GROWTH 

• ABSOLUTE THICKNESSES THIN DURING PRECHILL AND C/0 

• TENSION LOADING 

• a AT SHORTENING 

• RELATIVE THINNING - MC THINS MORE THAN LC OR ML THROUGHOUT 
DUTY CYCLE. FUNCTION OF 

. BIAXIAL RESTRAINT AT ML 

• NOTCH RESTRAINT AT LC 



Figure 62. 135% hg MCC Analysis Detail Review, Thinning 
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STEADY STATE 


e DISPLACED SHAPES INCLUDE THERMAL CHANGES 
TO GEOMETRY 



O DISPLACEMENTS MAGNIFIED TO ACCENTUATE THEIR 
EFFECTS 


Figure 63. Geometry-Scaled Distortions for a Typical Duty Cycle 
Basic Configuration at 135% hg 



This was a rather interesting result. If one ignores the biased loading 
effects,^ the large tens ion/compress ion thermal loading of a restrained beam 
could thin under tension and thicken under compression. The relative thinning 
under steady-state compression, time increment 6, is probably due to the 
land/channel transition notch effect and midland high biaxiality, plus the 
biased loading from the differential pressure. 

Figure 63 gives a distorted-shape perspective for the cycle, and shows that 
the general shapes are as expected, based on the hardware cross-sections. 

Table 4 presents an overall summary of the stresses and strains of the key 
load^ increments.^ This basically reinforces the general observations of the 
previous discussion on the biaxial restraints and shear conditions in the 
cross-section. It is interesting to note the large difference in mechanical 
strain across the channel wall from the 1090 K (1500 F) thermal gradients. 
This makes the hot-wall strain about twice that of the cold-wall strain. 
Figure 62 shows the yield stress versus temperature across the wall ranging 
from 20.68 MPa (3 ksi) at the hot surface, to 137.9 MPa (20 ksi) at the inner 
surface. This demonstrates the significant strength capability of the inner 
half of the channel wall that was alluded to as previously mentioned in the 
structural analysis. 

It also makes a hand-solution methodology, such as Ref. 12, more difficult 
with these kind of stress-strain gradients. This reference only assumes 
average^ temperatures and, therefore, average yield strengths are rapidly 
increasing. The channel remains in high tension during the remaining portion 
of the cutoff cycle. After test, steady-state conditions have the channel 
wall in a high tension stress in the axial direction with low stresses in 
plane in the hoop and radial directions. 

The primary purpose of Fig. 64 was to point out the stress state that occurs 
at the hot surface during steady state and what happens during creep relaxa- 
tion. The material yield strength at 1090 K (1500 F) is only slightly higher 
than the surface pressure stress (i.e., P^. = 2/3 FTY) . This means that the 
pressure stress is a significant part of the total stress state, and that the 
stress component in the radial direction must always equal the PjjQy value. 
At the steady-state condition, the top half of the channel wall goes through a 
creep relaxation process, which shifts the differential pressure primary load 
to the stronger (cooler), inner-half of the wall. At the hot wall, the stress 
state equilabrates to primarily a bulk compression equal to the Pugy pres- 
sure 14.54 MPa (2109 psi), plus a small component that develops a combined 
effective stress of 5.52 MPa (800 psi). This is a rather unusual stress con- 
dition in that rarely does a surface pressure equal twice the effective stress 
on the part. 

Figure 65 simply restates, in summary format, important observations made dur- 
ing this discussion. 
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TABLE 4. SUMMARY OF 135% hg APSAC PROGRAM RESULTS 


Steady State 


J 

1 / 




ML 


MC 



-£--X 


LC 


-170°F 
750° F 

1500°F 


CONDITION 


MID-CHANNEL (MC) 


CHANNEL-LAND 
INTERFACE (lc) 


MID-LAND (ML) 


General 

Comments 


2D Stress State 
2 Free Surfaces 
Plate/Beam Analogy 


2=-aAT =-.017 HM 
Range =-.0075 CW 

^x=1.2aAT, ^y=.4oAT c=1.3oAT 

^x=1.6aAT, ^y=.8aAT e=1.7aAT 


Notch Zone 
High Local -Shear 
Stepped Plate Analogy 


3D Stress State 
Free Surface Biaxial 
Restraint 


Steady State 
(INCR 6) 

Hotwall - 1500 F 
Coldwall - 750 F 


e 2 =aAT t< 2 ksi 


Cjj=1.5aAT, Gy=-.7aAT c=1.5aAT 

'x=1.2aAT, ^y=-.4cxAT, 
V=1.2aAT, e=1.73aAT 


e^=-aAT 


ey=e^ = -aAT e=1.8aAT 

Vo, ^j,=-otAT e=l.aAT 


Steady State 
to Cutoff 

HW-1047 F(INCR7 ) 
HW-739 F (INCR 8) 

HW-300 F(INCR9) 


Large Bending a, 
ave. axial a=-6 ksi 


Minimal Bending a 
ave. axial a=12 ksi 

No Bending a 
ave. axial a=28 ksi 


T < 2 ksi 

1=4. ksi 

0^=40 ksi 
T = 4. ksi J 


1 in. • 25.4 BUI 

1 psi * 6894.8 Pa 
•F .9^-460 


0 Notch 




TABLE 4. (Concluded) 


o 


-> Ev 


Steady State 



-170“F 

750”F 

ISOO^F 


CONDITION 

^ ■ ■ 

MID-CHANNEL (MC) 

CHANNEL-LAND 
INTERFACE (LC) 

MID-LAND (ML) 

Minimum Cutoff 
HW'-400®F(INCR1 ] 

®x*25 ksi. No Bending 
0 32 ksi 02=35 ksi 

No Bending Throughout, End 
(Constant Temp, , no Pressure 




V-’< 

ty=-.S% <'^^=1.4% 0 Notch 

c,.-4S 

11 

0 


• 8% C“l • 1 % 

c^=-1 .% 0 Notch e=1.3@NotcI 

Ey" • 2^ 

£ 2 = .35% Constant 


C2=.35» Constant 

Ey=.3" £ 2 =. 35*!; Constant 

e=.5% 0 Notch 


Strains 2X Mid-land or LC 
Interface except for Local 
Notch Zone 




Room Temp. 

No Bending Th 

•oughout End (70°F Constant, n 

) Pressure) 


INCR 13 


<si on Whole X-Section O En 

1 


70" F 

y:sr. 

c^=-.33% e^="7% 0 Notch 


IT =0 
xy 


C = 1% 

Cy=.3* -«• Cy=.7% 0 Notch 

ry=.2% 



^z=0 on Whole X-Section 

V,=.7% ^ y; =1.33% 0 Notch 

Aj Ajr 

e=.2% 

1 in. • 25.4 nn 

1 psi • 6894.8 Pa 




c-.6% g=1 .2% 0 Notch 


•F • ^ - 460 


/ 
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Figure 6A. 135% hg MCC Analysis Detail Review 
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• STEADY STATE RESULTS 

• CYCLES 1 AND 5 REVIEWED 

• SIMILAR RESULTS 

•DISPLACEMENTS AND STRAINS SHOW CONSISTENT RESULTS 

•MID-CHANNEL (MC) DISPLACEMENT AND STRAINS LESS THAN MIDLAND (ML) 

WHICH ARE LESS THAN THE LAND TO CHANNEL (LC) VALULS 

•AXIAL STRAINS (^z) UNIFORM AT POINTS ALONG THE HOOP DIRECTION 
•THINNING STRAINS (^x) “yz) 

•AT MC LOWEST RESTRAINT WITH LOW SHEAR RESULTING IN SMALLEST 

• AT ML, + H (“yz ~ 0) DEMONSTRATING HIGH BIAXIAL RESTRAINT 

•AT LC,“yz HIGH AND RESTRAINT IS MODERATE YIELDING THE LARGEST DISPLACEMENTS 
•CONTOUR, DIGITAL A GRAPHICAL PLOTS ALL CONEIRM IHIS CONDITION 

• THERMAL GRADIENT THROUGH THICKNESS LARGE (1500 F to 750 F) (1089 K to 672 K) 

•LARGE aAT GRADIENT - FACTOR OF 2 

•MAXIMUM e VARIES ALONG HOT WALL DUE TO VARYING RESTRAINT 

• HOT WALL CREEPS & RELAXES 10 a ^Oy , COMPONENTS ^ 

•STRESS GRADIENT THROUGH WALL TEMPERATURE DEPENDENI - a, Fjy 
NOT STRAIN DEPENDENT 

•GRADIENT FORCES MAXIMUM e ON HOT SURFACE 3% 


Figure 65. Observations of 135% hg Basic Geometry, SSME MCC 
APSAC Analysis Detail Review 



Figure 66 is a typical multiple duty cycle plot of the stress-and-strain his- 
tory of the element located midchannel at the hot wall. The efective stresses 
and strains are plotted with their sign obtained from the hoop-stress compo- 
nents (note: effective stress and strain are only magnitude quantities). 

The duty cycle increments and hot-wall temperatures are noted on the first 
cycle to help follow the duty cycle. The point first goes in tension during 
the prechill operation, followed by a rapid drop to a high compressive state, 
then reduces in stress magnitudes as the material strength decreases with tem- 
perature and, from creep/relaxation effects, to the maximum compressive strain 
(point 6). The stress decrease from increment 2 to 6 is mainly due to changes 
in yield stress with temperature. During the throttle back to 65% of thrust, 
the point moves to a tension stress state even though the total strain is a 
high compressive plastic strain. As cutoff starts, the element is driven to a 
high-tension, high-positive strain point. Finally, after the test duty cycle 
and the part is stabilized at room temperature, the element is in compression 
with a large tensile, plastic strain. 

Bach duty cycle progressively shifts to a more positive strain position. This 
is indicative of a thermal ratchet or cyclic-creep phenomenon and resultant 
channel-wall thinning. 
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Figure 66, Midchannel Hot-Gas Wall Multiple Duty Cycle 
Stress/Strain History 
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LIFE-ENHANCED DESIGNS 


General 


Task III of this study consisted of evaluating four life-improvement designs 
using the developed methodology. The requirements were to select two designs 
with minimum changes from the current MCC configuration; one best concept that 
could be a radical change, and one combination configuration or second-best 
concept configuration. In addition, concepts where coatings were required 
were declared outside the scope of the study. 

The total concept evaluation addressed manufacturing, thermal, structural, and 
operational considerations. Some of these are discussed in other sections of 
the report. The discussion herein is primarily limited to the structural 
considerations. The goals for the life enhancement were to: 

1. Reduce the cyclic creep deformation rate and extend the life to crack 
initiation 

2. Reduce the hot-wall surface and overall temperature. This lowers the 
thermal strains, increases the strength capability, and minimizes the 
sensitivity of the channels to creep or relaxation. This also mini- 
mizes the onset of surface-roughness changes that occur on the hot- 
gas-wall surface and increases the hot gas heat transfer rates. This 
condition further bootstraps local thermal conditions from a minor 
hot-spot to a significant local hot-spot that produces accelerated 
cyclic creep cracking 

3. Reduce the cyclic strain range and increase the low-cycle fatigue 
life of the part 

4. Maximize the available elongation or strain capability of the mate- 
rial by using designs with preferrential biaxial ratios, i.e., close 
to uniaxial stress states 

Figure 67 summarizes the four designs selected. The first design — contoured 
configuration — has essentially the same temperatures as the baseline configu- 
ration, but the thickness has been varied as a mirror image of the thickness 
changes that occur with cyclic creep. The thickness at the corner has been 
reduced by 0.038 mm (0.0015 inch) and thickened in the center by 0.025 mm 
(0.001 inch), making a total thickness change of approximately 10% of the nom- 
inal 0.71 mm (0.028 inch). The second design — increased number of channels — 
has reduced wall thickness and more channels. This design is similar to one 
of the proposed SSME uprated designs, and reduces the surface temperatures by 
over 83 K (150 F) . The third design — the keel rib — essentially adds a fin in 
the center of the channel to pull more heat from the hot wall. 
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Figure 67. Analyzed Channel Configuration 
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The channel height has also been increased to maintain the baseline 
crosssectional area. Structurally, this does not change the load path in the 
chan- nel wall markedly, since the rib is too narrow to allow the load to 
diffuse into the total thickness; but it does change the biaxiality strain 
picture and affect the rate of thinning. The lower temperatures for a keel 
rib design will also minimize the thermal strain state, but the channel still 
carries the same differential pressure-primary load. The fourth design 
concept— slotted mid-land— essentially interrupts the hoop continuity of the 
structure, thereby reducing the biaxial strains. The added notch is of 
concern, but its base is located in a cooler portion of the structure, which 
should reduce some of the local strain concentration. At this time, it is 
believed that a slot any thinner than 0.25 mm (0.010 inch) could not be 
fabricated on a repetitive basis with sufficient reliability. The channel 
width and hot-wall thickness of 0.75 mm (0.030 inch) and 0.63 mm (0.025 inch), 
respectively, provide opti- mum heat transfer capability, and with a 1.02 mm 
(0.040 inch) land width, the slotted ligament thickness of 0.38 mm (0.015 
inch) provides a thermally effi- cient channel/land ratio. The mid-land slot 
interrupts the hoop continuity of the hot wall, which will greatly reduce the 
thermally induced cyclic strain range in that direction. For flexibility, a 
slot depth of 1.02 mm (0.040 inch) was chosen. This creates a 0.38 x 0.38 mm 
(0.015 X 0.015 inch) ligament on each side of the channel to absorb the hot 
wall thermal displacements. During steady-state operation, the hot wall will 
expand and reduce the slot width by 0.025 mm (0.001 inch). 

In order to prevent high temperatures in the slot region, it was assumed that 
the slot was filled with a material such as felt metal or sintered screen 
mesh, which was attached by brazing at the bottom of the slot. A maximum 
operating slot/filler gap of 0.051 mm (0.002 inch) was assumed for the thermal 
analysis . 

The life enhancement designs should be viewed as a first-order approximation, 
since they have been optimized from only a thermal or structural standpoint. 
An actual design study would include further modifications and iterations to 
better optimize the life of the total design. 

Structural Analysis 

These concepts were evaluated using the analysis techniques that were devel- 
oped for the baseline analysis. The 135% hg case (hot-spot) was used for the 
evaluation, since it more accurately reflects the temperature conditions when 
channel cracking normally occurs. Three of the four concepts required addi- 
tional thermal analysis to optimize the geometry and furnished detail tempera- 
tures for the time slices of interest. The contoured channel design used the 
same temperatures as the baseline 135% hg condition. The jacket displacements 
and other boundary conditions were the same for all design analyses. Figure 
65 showed CRT pictures of the five designs studied, including the maximum 
surface temperatures during steady state from the heat transfer analysis. 
From a simplistic maximum temperature and resulting thermal strain condition, 
the keel rib, increased number of channels and slotted mid-land configurations 
have essentially the same hot-gas-wall temperature 987 to 1003 K (pi5 to 
1345 F) . The SSME baseline design and the contoured channel design are 
approximately 89 K (160 F) hotter at 1083 K (1489 F) . 
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The analysis procedure consisted of running the first 5 duty cycles plus a 15 
cycle extrapolation, resulting in detail plots for each design. These results 
were carefully reviewed to ensure modeling accuracy and general observations 
about which would be the best design. The rest of the analysis was then com- 
pleted. The primary methods of correlation of the designs were through cyclic 
creep life evaluation, LCF evaluation, maximum surface temperature, and a 
gross estimation of fabrication difficulty. Figure 68 furnishes a comparison 
of mid-channel wall thinning vs duty cycle for the various design concepts for 
the first 5 duty cycles, plus 15 cycle extrapolation. The life-enhancement 
designs were then extended to 120 duty cycles and beyond, by estimation, based 
on being bracketed by actual analysis data for the basic design 135% and 100% 
hg cases out to 120 and 100 duty cycles, respectively. These results are 
depicted in Fig. 69. The set of information available after the completion of 
the first five duty-cycle results for each life-enhancement design was util- 
ized as a test to determine how well the final results of the study could be 
extrapolated from this limited data. If good accuracy was obtainable, future 
sensitivity studies could potentially be made with only this information. The 
540-channel configuration and the contoured channel should give the best 
extrapolations, since they are only variants of the standard channel design. 
The keel rib and slotted configuration are more geometrically different, 
although the keel rib itself does not grossly affect the primary load path. 
Figure 69 shows the 120 duty cycle extrapolated results, and orders the con- 
figurations from the highest thinning to the lowest thinning as baseline 135% 
hg, slotted, keel rib, contour wall, 540 channels, and 100 percent hg base- 
line. The actual order of the life enhanced designs, as will be shown later, 
is only correct for the 540-slot design, the most similar design to the base- 
line configuration. Figure 68 shows another variable that was not completely 
assessed in this study. This is reflected by the rate of change of thinning 
which has not completely stabilized after 5 duty cycles for some of the con- 
figurations. The baseline 135% hg and the 540 channel design have the same 
rate of mid-channel deformation for between cycle 4 and 5, and therefore the 
extrapolation to 20. The other three design results are somewhat ratty, and 
would have required more duty cycles to completely stabilize. For this 
effort, this variation was judged relatively unimportant to the overall objec- 
tives of the study within the available time and cost restraints. Figure 70 
plots the peak strain at four critical locations in the models. The maximum 
strains across each configuration vary roughly as a function of maximum tem- 
perature, except for the slotted design, which does radically drop the strains 
due to the interruption of the hoop continuity. 

Two other plots were made to verify that the maximum rate of thinning was 
being used for the two designs that had variable-channel wall thicknesses 
(contoured and keel rib). Figure 71 is a plot of the thinning along the con- 
toured design channel wall, in a circumferential direction (Nodes 1 through 
9). As expected, the thinning rate is still the maximum at the center of the 
channel. Figure 72 is a similar plot for the keel rib design, and shows that 
the maximum thinning over a length that includes the keel rib and the channel- 
wall intersection with the rib. These results were also checked at the 100 to 
120 duty-cycle point and substantiated the methodology of using the same thin- 
ning concept for ranking the configurations. 
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Figure 70. SSME MCC Life Prediction Sensitivity Study Peak 
Strain vs Life Enhancement Configuration 
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Figure 71. Contoured Wall Thinning vs Nodal Location 
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Figure 72. Keel Rib Thinning vs Nodal Location 






The finite-element wall thinning analysis results of the five configurations 
through 120 duty cycles is presented in Fig. 73, The results are also extrap- 
olated beyond the 120-cycle point with dashed lines. The 5 duty cycle 
analysis thinning rates (Fig. 73) are significantly different than the 
120-cycle rates for all configurations. The ranking of the designs, by 

analysis, as a func- tion of thinning from maximum thinning to minimum 
thinning is the keel rib, baseline 135% hg, contoured, 5A0 channels, and the 
baseline 100% hg. Table 5 summarizes this information, plus other notable 
variables to allow assessment of all the configurations. The analyses were 
all made at a 135% hg, except for the baseline 100% hg. In all cases, the 
cyclic creep life is less than the typical LCF-life calculation. The keel-rib 
design and the baseline 135% hg are very close to each other from a cyclic 
creep standpoint, although the lower keel-rib temperatures would be less 
sensitive to hot-spots and surface roughing. The high thinning rates of the 
keel rib design are attributed to subtle changes in strain biaxiality, caused 
by the addition of the keel rib. 

The 80 to 115 duty-cycle life of the baseline 135% hg case is still consider- 
ably different than the actual cracking of 5 to 30 cycles. As previously men- 
tioned, the best estimation of the reason for this discrepancy is the poten- 
tial error in simulation of local strains in a hot-spot. Appendix discusses 
this in detail. The channel contouring seems to have a more significant 
effect for the first few duty cycles up to 60 or 80, but then has similar 
thinning rates to the baseline analysis. The 5A0-slot design is by far the 
best design of the more standard configurations. The slotted design has a 
high thinning rate, but because its biaxial restraint is closer to a uniaxial 
case, the allowable thinning for this configuration is assessed as twice the 
other configurations (10 to 16%). This configuration also has a significantly 
lower strain range. On the negative side, the slot-filler material and the 
heat transfer analysis accuracy for this model are in question. The filler 
material was added to eliminate unknowns associated with heat transfer charac- 
teristics within the slots. 

Copper and AMZIRC normally have more radial motion and thinning than is 
obtained with NARloy-Z. In some cases, such as the AOK chamber, NARloy-Z has 
shown large amounts of radial deformation and wall thinning prior to failure. 

As shown in Fig. 7A the amount of mid-channel radial deformation and mid- 
channel wall thinning is quite variable. The four basic channel wall cases- 
100% hg SSME, 135% hg SSME, increased number of channels, and contoured chan- 
nel wall — show similar and consistent solutions. The higher surface tempera- 
ture cases radially deform faster per cycle and the basic spread of data is 
reasonable. The slotted results furnish the most linear curve. The keel rib 
configuration has very small radial deformation, but one of the highest rates 
of thinning. 
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TABLE 5 


SSME MCC LIFE PREDICTION CONCEPT COMPARISON 



% Hg 

LCF LIFE, 
DUTY CYCLES 
(TYPICAL) 

CYCLIC 
CREEP LIFE, 
DUTY CYCLES 

MAXIMUM 

TEMPERATURE, 

F 

MAXIMUM 

TEMPERATURE. 

K 

FABRICATION 

DIFFICULTY 

COMMENT 

SSME 

100 

660 

290 TO 445 

1087 

859 

— 

BASELINE ANALYSIS 

SSME 

135 

300 

80 TO 115(1) 

1484 

1080 

— 

MAXIMUM HOT-SPOT CONDITION 

1 . CONTOURED 

135 

300 

100 TO 140 

1484 

1080 

LOW 


2. 540 CHANNELS 

135 

400 

170 TO 260 

1313 

985 

MODERATE 

SIMILAR TO SSME UPRATED 

3. KEEL RIB 

135 

400 

75 TO 110 

1337 

998 

MODERATE 

SIMILAR TO PROPOSED OTHER 
HARDWARE 

4. SLOTTED 

135 

550 

210 TO 320(2) 

1346 

1003 

HIGH 

NEW 

(1) 

ACTUAL 

CRACKING 10 

TO 30 CYCLES 





(2) 

TWICE THE ALLOWABLE STRAIN LIMIT USED 

FOR LOW BIAXIAL 

RESTRAINT CONDITION 
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WALL THICKNESS, t(%) 





It is interesting to note the highly nonlinear shape of the curve with 
cycling. During the initial sensitivity study previously discussed, the dis- 
torted results from the first 5 duty-cycle analysis was extrapolated to a 107» 
thinning condition. This required an accumulated duty cycle of approximately 
400 cycles, and a radial distortion of 3.81 mm (0.15 inch). As previously 
mentioned, trying to use estimated geometries, etc., led to conditions both in 
cyclic life and geometric distortions that were grossly inconsistent with test 
hardware history and/or detailed duty-cycle analysis. 

This type of plot is of interest for two other reasons. NASA has previously 
used these plot parameters for evaluation of measurements from test hardware 
(Ref. 2), and a closed-form solution technique for channel wall analysis 
(Ref. 12) uses radial deformation to relate to thinning. A markedly nonlinear 
relationship is not projected by the closed-form methodology, since the total 
analysis is based on a 1/2 duty-cycle consideration with no shakedown effects. 
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TECHNOLOGY IDENTIFICATION 


Recommended technology items to be considered for follow-on work are summar- 
ized below. These suggested areas requiring experimental investigation are 
based on the results of this analytical study. 

PROPOSED FOLLOW-ON EFFORTS 

1. HCC Analysis Improvements (Strain History and Extrapolation 
Methodology) 

2. Materials Laboratory Testing 

3. Hot Spot Analysis Development 

4. Midland Hot Wall Slot Manufacturing Verification 

5. Hot-Fire Test Program 
MCC AI4ALYSIS IMPROVEMENTS 

The current SSHE MCC life prediction study has demonstrated that cyclic creep 
and the use of a biaxial strain failure mode can be used to predict channel 
cracking using currently available analytical tools, i.e., FE (finite element) 
codes such as APSAC, extrapolation procedures, and appropriate failure cri- 
teria. This effort is a worthwhile step toward accurate MCC life prediction, 
but there are several areas where the life analysis procedure and methodology 
can be improved. Specific items include the following: 

1. Improved extrapolation methodology 

2. More acurate representation of material properties in FE 
codes 

3. Better postprocessing to mechanize tracking of thinning, 
strain, and stress history, etc. 

4. Investigation of the cold-wall fracture mode that can occur 
after hot-wall cracking 

5. Improved hand analysis methodology 

The current SSME MCC cyclic-creep analysis is an implementation of methodology 
developed in other LeRC supported work. The procedure includes a 5-duty-cycle 
FE analysis followed by an extrapolation of the geometry to the shape pro- 
jected at 15 additional duty cycles. The new shape is then used as a basis 
for additional analysis using the 5+15 additional duty cycles concept. The 
procedure is continued until failure occurs based on low-cycle fatigue or 
biaxial strain limits. 
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The resultant thinning for the basic SSME configuration is shown in Fig. 75. 
The basic technii}ue loses the strain history in the extrapolation procedure, 
since only the geometry shape is projected. This causes the rounded stepped 
shape at the first 1 of 2 cycles of the detailed duty-cycle analysis. The 
material gradually increases in maximum strain from geometric changes and 
biased loading during the duty cycle. This shakedown effect can be removed 
from the results by subtracting out the bias, but more accurate methods are 
recommended. 

An improved extrapolation methodology would add the strain history as part of 
the procedure. This can be accomplished in at least two ways. One is to 
extrapolate the strains as well as the geometry. The other is to use the mod- 
ified geometry as a rigid distorted shape relative to the original geometry 
and calculate the resulting strains. This loading would be a prestart condi- 
tion for a succeeding 5-duty-cycle analysis. Either methodology requires 
careful attention to the mechanized procedure in the FE code for accounting 
for the distorted shape throughout the total part history. This latter pro- 
cedure will be tested in the next month on Rocketdyne funds to partially 
prove-out the technique. 

The analysis methodology also can be improved by using cyclic stress-strain 
curves and more accurate creep relaxation information (discussed in a separate 
section) as well as a more accurate representation of the stress-strain 
curve. Currently, only a bilinear representation is utilized for the stress- 
strain curves. To demonstrate a necking or true tensile instability phenom- 
ena, the stress-strain curve representation must have a maximum point. There- 
fore, addition of a multilinear representation is required to the APSAC code. 
Addition of cyclic strain hardening would be desirable, but is not a critical 
item because the cyclic strain hardening for NARloy-2 is relatively small 
(<20%) , and its effect will be partially considered by using cyclic stress- 
strain curves. 

Another area for improving the evaluation and understanding of the computer- 
ized analysis would be some additions to the postprocessing plot and output 
code. Additional features are recommended to summarize data across the 5 duty 
cycles of information. Currently, all processing is done on a single incre- 
ment and the code has no features for collecting data for multiple increments. 

Plotting the element strain cycling loop, summarizing nodal displacements, 
thinning and element strain information are additional needs to minimize hand 
plotting and allow a more thorough understanding of results. 
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Another area where improvements are needed is the incorporation of cyclic 
creep into the hand analysis methodology. Rocketdyne has hand solution 
methods considering hot-wall temperature gradients, pressure, etc., for deter- 
mining factor of safety and LCF life. The cyclic creep methodology (Ref. 12) 
NASA CR-165585, "The Development of a Simplified Procedure for Thrust Chamber 
Life Prediction," can be upgraded to consider such items as hot-wall tempera- 
ture gradients, associated strains, and material property variations. This 
methodology could then be included as part of the hand solution for MCC chan- 
nel wall life and factor of safety. It can be used for initial sizing and 
sensitivity studies. 

MATERIALS LABORATORY TESTS TO IMPROVE THE 
ANALYSIS EFFORT 

Information on the cyclic stress-strain response, ductility, and creep/relaxa- 
tion of NARloy-Z needs to be developed to verify and calibrate the life analy- 
sis of the SSME MCC channel wall liner for the cyclic creep failure mode. The 
cyclic response of the material would be investigated through the determina- 
tion of cyclic stress-strain curves. Large mean strain effects should be con- 
sidered to relate to the hardware strains. Availability of ductility would be 
investigated through the testing of multiaxial specimens to obtain elongation 
at onset of tensile instability. Creep relaxation data would be obtained from 
conventional test specimens. Information is required over a temperature range 
of -400 to 1500 F. Some background on the development of the desired informa- 
tion follows. 

Cyclic stress-strain curves would be obtained over a wide temperature range to 
a maximum cyclic strain range dictated by analysis. Consideration would be 
given to the strain rate corresponding to hardware operation. The program 
testing would be limited to isothermal uniaxially loaded, deflection control- 
led tests. The curves would be run to get information on the material strain 
hardening exponent used in instability analysis for determining biaxial strain 
limits and to provide elastic/inelastic stress/strain data for the finite- 
element analysis of the main combustion chamber channel wall. 

Ductility defined by elongation would be measured under isothermal conditions 
using multiaxial monotonically loaded specimens. A review of the literature 
indicates biaxial loading through the use of cylindrical pressure vessels or 
cross specimens is feasible (Ref. 13). The influence of cycling and noniso- 
thermal conditions is recognized, but investigation would not be part of this 
program. The available elongation obtained from test would be used to evalu- 
ate elongation at onset of tensile instability. 

Creep/relaxation data would be used to support life analysis of combustion 
chamber channel wall hot spots. At the present time, no information exists on 
NARloy-Z creep behavior in the temperature regime to 1089 K (1500 F) . Data on 
creep rupture ductility, creep deformation, and relaxation would be obtained 
from tests of uniaxially loaded specimens. A theoretical adjustment to the 
data would have to be made to account for material strain cycling and 
multiaxiality. 
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DEVELOPMENT OF HOT SPOT ANALYSIS TECHNIQUE 


The 135% hg thermal load case being used to evaluate and compare the various 
MCC life enhancement design concepts represents an operating environment that 
generates a 1089 K (1500 F) hot-gas wall temperature 167 K (300 F) above nomi- 
op®J^^bing conditions). This condition typically occurs due to injector 
anomalies, and on the actual MCC hardware shows up as local hot spots that 
generally have a 2:1 to 4:1 elliptical pattern. 

The current analysis technique used for the finite element channel/land/close- 
out model treats the overtemperature condition as occurring axisymmetrically 
around the entire ID circumference. This generates more total circumferential 
growth of the hot wall that is restrained by the structural jacket and the 
EDNi closeout; however, the generated strain is shared uniformly by all the 
channels around the circumference. In reality, the local hot-spot area chan- 
nels are much weaker than the remaining normally heated hot-wall channels 
around the circumference. Although there is less total circumferential growth 
to be restrained, the hot spot is a weaker material that acts like a hole in 
the total liner with a #/£ concentration factor causing local distortions 
that force more strain into the local weak zone. 

A preliminary APSAC analysis of a plate with a 4:1 elliptical hot spot (135% 
hg hot-wall temperature) showed peak strains and strain ranges that were a 
factor of 2 to 4 times higher than the remainder of the plate (100% hg hotwall 
temperature) . 

The axisymmetric modeling technique only shows an increase in peak cyclic 
strains of 15 to 20% when comparing the 100% hg and 135% hg thermal load 
cases. Therefore, the local hot spot appears to provide a mechanism that 

could generate significantly higher strain and midchannel deformation rates 
that would cause the cyclic creep failure mode to occur in a number of cycles 
that would more closely match the hardware failures observed to date. 

The optimum hot-spot analysis technique would utilize a three-dimensional pro- 
gram with inelastic, creep, and incremental load case capabilities. The cost 
and difficulty involved in completing this type of analysis would not be prac- 
tical as a standard analytical tool, and this follow-on effort would be 
directed toward the development of a simplified analysis approach that would 
simulate the increased localized hot-spot effects. 

The current 2-D model of the channel cross section would be supplemented by a 
second 2-D plate model that would generate the membrane hot-spot effects. 
Simplistically, this could be an APSAC plate model of just the hot-gas wall, 
or could be a more involved ANSYS plate model incorporating multiple layers 
that would simulate the hot wall, lands, and cold wall. The plate model hot- 
spot effects would be imposed on the channel cross-section model in the form 
of displacement boundary conditions or as pseudo material properties such as 
E, etc., and then the model would be subjected to the typical 5— cycle 
PSAC multi-increment duty cycle analysis plus the 15-cycle extrapolation to 
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establish the cyclic strain range and channel deformation rate (tensile insta- 
bility mode). To verify the technical adequacy of the simplified approach, it 
is recommended that a one-time 3-D model be generated. The 3-D model would be 
used to compare with the 2-D model(s) results to ensure that the hot-spot 
effects have been properly simulated. Once good correlation is achieved 
between the 2— D and 3-D models, the simplified approach could be used as a 
standard analytical tool for determining life of combustor hot-wall channels 
that are subjected to local hot-spot injector anomalies. 

HOT-FIRE TEST PROGRAM 

A hot-fire test program is proposed using 20 K size test hardware to substan- 
tiate life prediction techniques and evaluate life-enhancement channel rede- 
signs. The proposed program consists of five technical tasks: 

J. ~ Establish measuring techniques to define channel geometry history 

(deformation/thinning) of test hardware. 

- Design and fabrication of test hardware. A minimum of two life- 
enhanced channel concepts will be utilized and compared to the baseline SSME 
MCC design. 

Task — III — Hot- fire testing, including facility build and/or modification and 
installation of test hardware. 

Task IV - Data correlation and empirical cyclic life relationships established, 

^ - Analytical life cycle model upgrading culminating in a versatile com- 
puter code (to NASA). 

A new combustion chamber with axisymmetric coolant passages would be designed 
and fabricated for the hot— fire testing. The combustion chamber would be 
designed with different coolant passages in three or more sectors of the cham- 
ber. Each sector would be independently manifolded and instrumented to meas- 
ure the total section heat load. Thermocouples would be installed, at 
selected locations at the throat plane to measure the hot gas side surface 
temperature and the temperature at the midpoint of the lands. The test matrix 
would include varied operating points and injector configurations to allow 
investigation of the effects of local mixture ratio zones on the wall, chamber 
pressure, established during Task I, coolant flowrates, and coolant fluids. 
Posttest hardware inspections would be performed to allow correlation of chan- 

life as a function of the accumulated exposure of thermal and structural 
loading conditions. The proposed program effort would colminate in the devel- 
opment of a versatile computer-code capable for cyclic life prediction of any 
convectively cooled channel-configured heat exchanger. 
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LIFE ENHANCEMENT MIDLAND HOT-WALL SLOTS (LIFE ENHANCEMENT 
CONCEPT NO. 4) MANUFACTURING VERIFICATION 

A manufacturing feasibility demonstration is proposed to support life enhance- 
ment Concept No. 4 (Fig. 75). This concept relaxes the biaxial strain in the 
tangential direction, thereby increasing the MCC cyclic life. This life 
enhancement concept creates a substantial manufacturing challenge. Manufac- 
turing areas of concern are: (1) machining narrow slots of 0.010-inch width, 
(2) identification of midland locations and numerical indexing, and (3) 
machining on an internal contoured surface (ID of the MCC). The latter 
requires a special machine or adaptation to current types of milling machines. 

The proposed experimental program consists of four Tasks: 

Task I 


Verification of slot dimensional control and cutter life. This effort would 
be conducted on inexpensive flat plate NARloy-Z or AMZIRC (OFHC copper-base 
alloys) samples. 

Task II 


Search for cutting machine (or adaptation of existing machines) capable of 
cutting slots on the MCC contoured internal hot-wall surface (MCC ID). The 
milling cutter needs to be affixed to an arm with a minimum reach of several 
feet. It must also be capable of contour tracing and indexing with extreme 
accuracy. 

Task III 


Demonstrate machining of midland hot-wall slots on existing SSME design full 
size hardware (Fig. 76). As a minimum, the axial location of the slots will 
be from 4 inches upstream of the MCC throat to 1 inch downstream of the geo- 
metric throat. Inspection techniques will be established to identify slot 
dimensions and location relative to midland and adjacent coolant channels. 
Quality control will be emphasized. 

Task IV 


Evaluation of thermal barrier slot filler materials. Concepts to be evaluated 
will be low-density flexible filler materials, such as felt metal construc- 
tion. Usable surface temperature range will be >2480 K (4000 F) and accom- 
modate a LOX/H 2 combustion gas environment. The purpose of the slot filler 
is to retard or negate convective heating of the slot side walls. Bonding 
(brazing) of the selected slot filler materials at the slot-base is to be 
established during this task effort. 
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CONCLUSIONS 


1. Cyclic-creep in the form of channel wall thinning is the dominating factor 
in defining the reusable life of channel configured combustion chambers. 

2. The predicted cyclic-creep life is typically between 30 and 50% of the 

predicted low-cycle-fatigue (LCF) life. 

3. The analysis methodology used in this study predicts the correct failure 
mode, including midchannel cyclic-creep deformation patterns. 

4. The predicted midchannel wall thinning rate is less than indicated by 

hot-fire test results. However, a preliminary 3-D hot-spot analysis indi- 
cates a potential for generating higher local strain levels and an accel- 
erated thinning rate. 

5. The most feasible approach to increasing the SSME-MCC usable life is 

increasing the number of coolant channels. Increasing the number of chan- 
nels from 390 to 540 doubles the usable life. Channel dimensions (channel 

width, land width, and wall thickness) are typical of those used for 

advanced engine design studies. 

6. A maximum usable SSME-MCC life benefit of a factor of three was shown for 
the maximum risk design which incorporated two life enhancement features. 
This included increasing the number of channels from 390 to 460/and slot- 
ting the midland hot-gas wall to minimize the transverse strain. 
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APPENDIX A 


SSME-MCC EXPERIMENTAL BACKGROUND 

Substantial hot-fire testing has been conducted on the SSME and 40K main com- 
bustion chamber with similar configurations. As a result of this testing a 
considerable number of thermal cycles has been accumulated. This section 

"'“‘"I 

SPACE SHUTTLE MAIN ENGINE 

The following case histories are presented indicating the hot-gas wall thermal 
crack histones as a result of Space Shuttle Main Engine testing. 

MCC 0008, Engine 0009 

Main combustion chamber 0008 accumulated 21 tests and 9000 seconds total oper- 
llllU »^he time of its first noted hot-gas wall crack. The crack was 

located approximately 2.54 x 10~2 m (1 inch) above the throat centerline and 

rathe“r"Th^^ perpendicular to the direction of hot-gas flow 

endured initiation, that area of the ch^ber had 

endured discoloration and blanching, with a general trend of increasing rough- 
ness, throughout the chamber’s test history. Roughness readings of the area 
increased from the 7.62 x lO"! to 1.016 x 10"6 m (30 to 40 pin.) range 

aLlV^tfr * *"®sts immedi- 

in Wh' crack initiation. Radiographic examination following the test 
in which the crack was noted revealed no signs of obvious contamination. 

Following two additional 520-second flight profile tests, two more longitudi- 
1 cracks were reported in the same general plane as the first, each down- 
other elements. All three cracks were 5.08 x 10-5 to 1.27 x 10-4 m 

i^th^testine°°^‘T"''^^H-T^^^ (0/030-inch) long at this point 

P Vrol^"the ^cttVaT^^arLTwV^^^^^ 5o1 S? 

sln"u'cooi°an't channeT""" horizontal crack in a 

next fe°sf'‘8T3^\^n” T <^racks. Following the 

Ki 823 seconds, the chamber had a total of 10,900 seconds with five 

th/ oYh.'r <='■“''' ">= “"<=»*"6ed i„ s-AVl how! 

in-2* M . uf . f existing cracks had increased to approximately 2.54 x 
ieiehJnr- length. In addition, one of the cracks had developed a 

neighbor, approximately two or three channels away at a slightly different 
eig . Liner roughness in the immediate area was only 2.0324 x lO"® m (80 

cool ^ locally combust-i.on zone helped 

MCC 0007, Engine 0008 . 

tes^‘"^??Ar hot-gas wall crack very early in its test history (5 

seconds) and none since. The crack developed from an area which 
exhibited only minor discoloration after the first main stage test (third 
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hot-fire test) and relatively rough surface following the second main stage 
test. On the following test, the fifth (third main stage), there was a small 
eyelid-like channel rupture approximately 8.89 x 10“^ m (0.035. inch) Jong in 
the converging radius of the combustor zone. The surrounding area was very 
rough, approximately 5 x 10”^ m (200 yin.). Radiographic examination 
1.27 X 10"^ m (5 inches) fore and aft of the throat revealed no obstructions 
in the cooling channels. 

The general area of the crack has exhibited intermittent overheating since 
crack generation. This has apparently not been entirely eliminated by 
enlargement of the three overhead film coolant hoyes in the main injector 
primary faceplate. 

Another area of the chamber exhibited higher than normal roughness subsequent 
to the second main stage test, and finally reached severe enough proportions 
to warrant coolant hole enlargement. The one hot fire test since shows the 
area to have apparently been cured of overheating. To date, the chamber has 
accumulated 21 tests and 6750 seconds hot-fire duration with more than 1400 
seconds at 107% RPL or higher. 

MCC 2007. Engine 2008 (2018) . 

The first instance of hot-gas wall distress on this MCC was noted after its 
seventh test (fifth main stage) and 1020 seconds. Not coincidentally, this 
was the first long-duration (520 seconds) test and first test in excess of RPL 
(102%). Following two additional long-duration tests, the chamber exhibited 
two areas of severe roughness. Four film coolant holes were enlarged, which 
totally relieved the surface distress. No further signs of distress were 
noted during the remainder of the engine’s 102% RPL certification series. 

Following rebuild of engine 2008 as 2108, the engine was employed for a series 
of FPL tests. Two added distressed areas required additional film cooling. 
In one area, the remedy worked totally, however, in the second area, the dis- 
tress appeared even rougher after film coolant hole enlargement and exhibits 
fine surface cracking. 

MCC 0005. Engine 0005 

Numerous blanched areas appeared randomly dispersed throughout the convergent 
and throat sections of the chamber. Five areas of parent metal erosion 
appeared on the second hot-fire test on. the. engine. Maximum roughness in 
these areas was 5.08 x 10~^ m (110 yin.). 

Surface cracking appeared after the next test in two of the areas eroded on 
the previous test. One crack originated in the convergent region, and one 
just upstream of the throat. The fourth test resulted in rough blanching 
beneath 12 elements; two elements have two blanched areas each. Two more 
cracks appeared in and just above the throat beneath the convergent region 
crack. All three cracks penetrated the hot-gas wall, and all three were in 
different, adjacent channels. 
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Two tests later, there were two additional areas showing surface cracking. 
The next four tests resulted in no apparent further degradation. The fourth 
test ended in a main injector failure which eroded the hot-gas wall and coated 
it with flame spray and slag. 

MCC 0005. Engine 0006 

The first two tests on this engine resulted in no reported anomalies. After 
the third test, there were two areas of moderate to severe blanching. The 
next test resulted in a new through-crack underneath one element. There was 
no obvious change to the other cracks. Flame spray on the chamber came from 
FPB wall burn-through. 

Two tests later, there were five areas of light to moderate blanching, and two 
additional tests later, it was the same. Three tests later, there was one 
moderate to severe blanched spot. Following ten additional tests, primarily 
relatively short— duration, resulted in no more distressed areas, and all 
cracks appeared unchanged. 

RECENT SSME MCC TEST HISTORY 

Main combustion chamber (MCC) hot-gas wall failure in the form of channel 
cracks in the more recent phase of the test program is essentially the same as 
in the earlier phases and subscale testing. A noticeable trend during full 
power level (FPL) testing has been for cracking to occur faster in terms of 
number of engine tests and test durations (Tables 6 and 7). Cracks have 
developed in as little as three tests on a chamber during the FPL test phase. 
It is presumed that the crack failure modes are similar to those experienced 
on subscale chambers and earlier testing of the SSME, but no destructive anal- 
ysis has been performed on an FPL MCC. 

Circumferentially, liner distress most frequently occurs downstream of main 
injector baffle elements due to decreased mass flux in these areas. In 
attempt to minimize this distress, it has become standard practice to enlarge 
the 10 main injector boundary layer coolant holes adjacent to the baffles fol- 
lowing the first or second test of an FPL engine. 

One new form of hot gas wall distress, known as pimpling, has been noted on 
several recent MC^ operating at FPL. Circular areas of raised material up to 
about 1*27 X 10“ m (0.050 inch) base diameter and 5.08 x 10“^ m (0.020 
inch) high appear in hot spots where blanching or roughness has already 
occurred. A small pinhole may open at the peak of these pimples, producing 
what resembles a tiny volcano. Enlarging appropriate main injector film 
coolant holes seems to help prevent pinholes from developing and generally 
relieves this distress similarly to the case of blanching and channel cracks. 

The following case histories are presented indicating the hot-gas wall thermal 
distress histories of some recent high-time SSME main combustion chambers. 
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TABLE 6. MCC INITIAL CRACK HISTORY-RECENT ENGINES AT 
FULL POWER LEVEL (FPL) 


MCC/U/N 

ENGINE NO. 

TESTS BEFORE 
FIRST CRACK 

TEST TIME 
(SECONDS) 
NOMINAL/FPL 

CIRCUMFERENTIAL CRACK LOCATION 
(DOWNSTREAM OF OUTER ROW 
INJECTOR ELEMENT NO.) 

0010 

0010/0107 

6 

735/440 

36 (BAFFLE) 

9004 

0204 

7 

750/440 

56/57 

2011 

2010 

7 

2600/1890 

76 (HIGH VELOCITY AREA) 

9005 

2013 

3 

600/760 

36 AND 68 (BAFFLES AND 
LONGITUDINAL WELD) 

2015 

2014 

6 

1910/1310 

36 (BAFFLE) 

0007 

0008/2308 

5 

1145/945 

68 (BAFFLE AND LONGITUDINAL 
WELD) 

2015 

2011 

4 

128/0 

53 

2017 

2012 

6 

802/390 

76 (HIGH VELOCITY AREA) 

4001 

2017 

4( PINHOLE) 

1094/0 

72 

2014 

2015 

NO CRACKS 



*85 INJECTOR ELEMENTS IN OUTER ROW 
TRANSFER DUCT (COUNTING CLOCKWISE 

WHERE NO. 1 IS AT FUEL TURBOPUHP HOT-GAS 
LOOKING FORWARD). 


TABLE 7. 
MANNED 

MCC INITIAL CRACK HISTORY-RECENT ENGINES, FIRST 
ORBITAL FLIGHT (FMOF) ENGINES AT RPL OPERATION 

MCC/U/N 

ENGINE NO. 

TESTS BEFORE 
FIRST CRACK 

TEST TIME 
(SECONDS) 
NOMINAL 

CIRCUMFERENTIAL CRACK LOCATION 
(DOWNSTREAM OF OUTER ROW 
INJECTOR ELEMENT NO.) 

0006 

2004 

6 

220 

21 

2005 


NO CRACKS 

3840 


2004 

■■ 

NO CRACKS 

3760 


0009 

2007 

NO CRACKS 

4210 



13A 















MCC 2015. Engine 2014 (FPL) 


The 10 film coolant holes adjacent to baffles were enlarged as a precautionary 

<^®sts on this engine 

resulted in blanching which required the enlargement of nine additional film 
coolant holes. A surface crack developed after the fifth test in a previously 

PI ^ ^2.3 X 0.011 inch) by the sixth 

test. Blanching increased such that two more film coolant holes were enlarged 
after the seventh test. 

seconds, blisters appesred In a blanched 
throuEh cracks developed after the ninth test, the next test, one 
the through cracks grew and two merged and extended. Pinholes and pimples 

anrhactL f ^he 16th test (3868.1 seconds) in a previously blanched Lea 
and became four small holes with a crack following the 20th test. A channel 

lit ^ "‘°‘^® Blanching has been evident 

throughout this MCC s test history to date (36 tests, 11590.31 seconds) 

MCC Unit 2011. Engine 2010 (FPL) 

Blanching (measuring lOOpinch surface roughness) first occurred after the 
third test and 601.5 seconds. Blanching increased in five areas such that 
film coolant holes were enlarged by the fifth test. The next test revealed 

rf blanched areas became rougher and one area became 
eroded and cracked because of element contamination. Subsequent tests 
revealed more blanching and the need to open more film coolant holes. A 

Ki ^ u (1/2-inch) long crack developed in an area that was previously 

blanched after the seventh test. Five tests later (seven tests, 2601.5 

•" * * 0.01 inch) developed 

and stabilized. A new crack developed adjacent to the previous crack six 
tests later. 


Six cracks and two pinholes were noted after the 17th test. The cracks 
extended longitudinally and the pinholes stabilized after six tests. Pimples 
appeared after the 27th test. A crack connected the pinholes after the 37th 

seconds. Two new cracks occurred after the 40th test, one 
at the throat and one in the upper convergent region. 

MCC 0006, Engine 2004 (FMOF) 


instance of hot-gas wall distress, which was minor, was noted early 
® history (four tests, 5.6 seconds). A surface crack at 

the noticed following the fifth test. Following the eighth test, 

^ surface crack, when viewed under magnification, 

showed molten Narloy-Z and slight erosion. 

After six additional tests (14 tests, 321.4 
revealed a channel with a small 7.94 x 10“^ 


5.08 X 10 


-2 


seconds) , another hot spot 
m (1/32 inch) long rupture 


m (2 inches) above the throat centerline and a total of five 

a blTshL'^^^^ "* channel cracks. Blanching with 
LxLV this area occurred during the next test. After the 16th test, 
here was some disturbance at a longitudinal weld, a 0.51 mm (0.02 inch) hole, 
and some cracking was visible in this area. 
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Seven test later, the cracks seemed to stabilize. A new 7.94 mm (5/16 inch) 
long crack was discovered after 30 test and 4587.5 seconds. This crack propa- 
gated to 9.53 mm (3/8 inch) long. A blanched area with high roughness (300 
yin.) developed a small crack after the 44th test. Blanching was evident 
between the 44th and 49th tests. 

SSME SUBSCALE THRUST CHAMBER 

A 1.8143 X 104 kilograms (40,000 pound) thrust subscale SSME thrust chamber 
was fabricated to demonstrate the cyclic life of its larger counterpart. The 
thrust chamber is designed to operate at 2.0684 x 107 pascal (3000 psia) with 
liquid oxygen and fuel-rich O 2 /H 2 combustion products for fuel. The com- 
bustion chamber is copper alloy with a nickel electro-formed closeout to com- 
plete the channel wall flow passages. These coolant passages duplicate the 
geometry of the full scale thrust chamber to obtain realistic cycle life data. 

A staged combustion cycle was employed to simulate accurately injector com- 
bustion chamber heat transfer relationship. The primary goals of the program 
were to: (1) accumulate 100 main stage thermal cycles on a hydrogen cooled 
combustion chamber, and (2) demonstrate the equivalent thermal strain damage 
of 60 SSME missions on the main injector oxidizer posts. 

The combustion system assembly (Fig. 77A) consisted of a single preburner 
assembly, a hot-gas transfer manifold, a main injector, and a main combustion 
chamber. The preburner assembly and hot-gas transfer manifold designs were 
workhorse in nature; however, the main injector and main combustion chamber 
incorporated key design features of the flight-type SSME hardware. 

The main injector assembly (Fig. 77B) consisted of 58 recessed coaxial main 
injection elements, 3 baffle injection elements, and a transpiration cooled, 
dual faceplate configuration. The injection elements and faceplates dupli- 
cated the full scale injector design relative to propellant injection param- 
eters and cooling. Demonstration of the injector body and oxidizer post 
material to withstand the thermal cycle strains induced by ducting both high- 
temperature (hot-gas) and cryogenic temperature (liquid oxidizer) fluids 
through the injector was one of the two design verification items included in 
the subscale program. Main combustion chamber test hardware includeded both 
water- and hydrogen-cooled units and are shown in Fig. 77C and 77D, 
respectively. 

The 1.8143 x 10 kilogram (40K pounds) subscale primary goals were met by 
accumulating more than 100 main stage thermal cycles on the regeneratively 
cooled chamber. 

Analysis of Test Results and Posttest Hardware Evaluation 

The subscale program* s primary goals were met by accumulating more than 100 
main stage thermal cycles on the first regenerative chamber (U/N 2). 

The service life of the SSME main combustion chamber, as demonstrated by 40K 
subscale chamber hot-fire testing, is 103 starts and cutoffs. This is based 
on the demonstrated 109 thermal cycles of the subscale chamber (U/N 2) prior 
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LED CHAMBER 


to attaining a low cycle fatigue crack in the NARloy-Z liner. The life anal- 
ysis utilized nominal heat transfer operating conditions of the regenerative 
chamber, where the average heat load during the test series was established to 
be 95 percent of the calorimeter chamber heat load. 

Regenerative chamber thermal escalation was noted by an increase in total heat 
load during cyclic testing. It was determined that the heat load could be 
reduced to a nominal value by polishing the hot-gas wall surface with a crocis 
cloth. It was also observed that local regions located 1,016 x 10“^ m (4 
inches) downstream of the injector and at the throat were roughened more than 
other areas between the injector and throat. All areas downstream of the 
throat were unaffected. Posttest hot-gas wall surface finish values for the 
regenerative chamber U/N 2 are summarized below: 


REGION 

SURFACE FINISH 
METERS 

MICROINCHES 

NOZZLE 

5.08 X 10"^ 

20 

THROAT 

2.54 X 10"® - 4 X 10"® 

35 TO 45 

COMBUSTION ZONE 

2.54 X 10"® - 4 X 10"^ 

35 TO 45 

LOCAL HOT SPOTS 

2.54 X 10"® - 4 X 10"® 

35 TO 45 

OVER MID-CHANNEL AND MID-LAND 

6.35 X 10"® - 8.89 x 10"® 

250 TO 350 

OVER CHANNEL TO LAND CORNERS 

1.524 X 10"® 

60 


The difference in surface finish over land to channel corners in the highly 
roughened regions, is attributed to polishing since these surfaces were 
slightly deformed toward the hot-gas. 

The LCF history of the two chambers is noted below: 


CHAMBER NO. 

TOTAL NUMBER 
OF TESTS 

LOCATION 
LFC CRACKS 

CHANNEL NO.* 

THERMAL 
CYCLE NO. 

002 

118 

10 INCH 

98 

109 



THROAT 

98 

112 



THROAT 

121 

112 

001 

57 

THROAT 

121 

53 


On the 40K chamber U/N 2 (Fig. 78) two hot-gas wall cracks appeared early in 
the test program, which were related to a manufacturing defect and a material 
defect; noted after 10 and 29 thermal cycles, respectively. Continued testing 
resulted in fatigue cracking after 109 and 112 thermal cycles. Further cycl- 
ing resulted in slight growth of existing cracks before the chamber was 
retired after a 118 thermal cycle test program (Fig, 79 and 80). 
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Figure 79 . 40K Chamber Section Through a Hot Spot 

Near the Throat 



Figure 80 . 40K Chamber Section Showing Typical Channel 

Deformation Near the Throat 
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On the 40K chamber U/N 1 a low-cycle fatigue through crack appeared in the 
throat region after 53 thermal cycles. The chamber was then tested for four 
more additional thermal cycles. The crack grew from approximately 1.27 x 
10"^ m to 1.91 X 10“^ m (0.5 to 0.75 inch) in length. The 57 thermal 
cycles completed the test program on the 40K chamber. 

3.3K Program 

The 3.3K chamber (Fig. 81) underwent 1013 hot-fire test cycles prior to leak- 
age through the hot-gas wall. The failure analysis summary, MPR 74-904, 
attributes failure to fatigue initiating from the coolant side of the hot-gas 
wall surface. The distortion and cracking are symptomatic of significant 
temperature variations between identical axial locations. 
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APPENDIX B 


MATERIAL PROPERTIES 

SSME MCC (REF. lA) MATERIALS DEFINITION 

NARloy-Z is a copper base, silver, zirconium-containing alloy developed by 
Rocketdyne specifically for high heat flux combustion chamber applications. 
The material has excellent thermal conductivity and can develop good strength 
and ductility up to 922 K (1200 F) by a simple heat treatment. The slotted 
liner requires a material which has good strength at elevated temperatures to 
resist coolant pressures and superior thermal fatigue resistance to accommo- 
date the severe thermal strains generated during hot firing. NARloy-Z has 
demonstrated that it offers the best combination of characteristics to satisfy 
these requirements. 

Electrodeposition of nickel (EDNi) having an initial barrier of electro- 
deposited copper (EDCu) for hydrogen environment protection was selected as 
the method for closing out the coolant channels in the slotted liner. The 
deposited structure has excellent bonding actions with the NARloy-Z liner 
material and makes a good structural closeout. The properties of electro- 
deposited nickel are good making it suitable for containing the coolant pres- 
sures in the channels while providing a transition metal for welding the jac- 
ket to the liner. 

EDCu is deposited initially to provide a barrier to the embrittling effects on 
EDNi of high-pressure hydrogen in the coolant passages. This barrier has 
imparted excellent low cycle fatigue life to EDNi in the presence of hydrogen. 

The material requirements to support the MCC life predictions for this study 
are summarized in Table 8. The materials considered are NARloy-Z, electro- 
deposited-nickel (EDNi) and electrodeposited-copper (EDCu). Available test 
data were compiled (Ref. 15) and the material properties curves extrapolated 
to the ranges of interest. In the absence of test data, various prediction 
methods were used. 

1. Thermal-Physical Properties 

The elastic .modulus and Poisson* s ratio data for NARloy-Z had to be extrapo- 
lated from test data between 978 to 1089 K (1300 to 1500 F) . For EDNi and 
EDCu, the elastic modulus and Poisson’s ratio was extrapolated between 29A to 
78.6 K (70 to -320 F) from test data. Thermal expansion data was available 
for all three alloys in the temperature range of interest. 

2. Structural Properties 

Tensile test data were available for all three materials. However, stress- 
strain curves at various temperatures were also needed. MIL-HDBK-5 suggests 
two methods for determining stress-strain curves (Ref. 16). They are the 
Strain-Departure Method and the Ramberg-Osgood Method. The elastic modulus 
and the yield strength at the required temperature plus stress-strain curves 
at other temperatures to determine a representative shape are needed. 


143 



144 


TABLE 8. MATERIAL PROPERTIES SUMMARY 


MATERIAL 

PROPERTY 

TEMPERATURE, 

K 

TEMPERATURE , 
F 

NARloy Z 

ELASTIC MODULUS 

78 TO 1089 

-320 to 1500 


POISSON'S RATIO 

78 TO 1089 

-320 TO 1500 


THERMAL EXPANSION 

78 TO 1089 

-320 TO 1500 


TENSILE PROPERTIES 

78 TO 1089 

-320 TO 1500 


STRESS-STRAIN DIAGRAM 

78 TO 1089 

-320 TO 1500 


LOW CYCLE FATIGUE 

78 TO 1089 

-320 TO 1500 


STRESS EXPONENT n 

644 TO 1089 

700 TO 1500 


MATERIAL CONSTANT A 

644 TO 1089 

700 TO 1500 


CREEP-RUPTURE 

694 TO 1089 

700 TO 1500 

EDNi 

ELASTIC MODULUS 

78 TO 422 

-320 TO 300 


POISSON'S RATIO 

78 TO 422 

-320 TO 300 


THERMAL EXPANSION 

78 TO 422 

-320 TO 300 


TENSILE PROPERTIES 

78 TO 422 

-320 TO 300 


STRESS- STRAIN DIAGRAM 

78 TO 422 

-320 TO 300 


LOW CYCLE FATIGUE 

78 TO 422 

-320 TO 300 

EDCu 

LASTIC MODULES 

78 TO 478 

-320 TO 400 


POISSON'S RATIO 

78 TO 478 

-320 TO 400 


THERMAL EXPANSION 

78 TO 478 

-320 TO 400 


TENSILE PROPERTIES 

78 TO 478 

-320 TO 400 


STRESS-STRAIN DIAGRAM 

78 TO 478 

-320 TO 400 


LOW CYCLE FATIGUE 

294 TO 617 

70 TO 617 


FIGURE NO. 


8? 

83 

84 

85, 86 
87 


91 

92 

93 

94. 95 

96 

97 
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Figure 82. NARloy-Z Elastic Modulus As a Function of Temperature 
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Figure 83. NARloy-Z Poisson’s Ratio as a Function of Temperature 
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Figure 84. NARloy-2 Thermal Expansion as a Function of Temperature 
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Figure 85. 

NARloy-Z Tensile Strength as a 

Function 

of 

Temperat 
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Figure 86. NARloy-Z Ductility as a Function of Temperatue 
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Figure 89. 


NARloy-Z Stress Exponent 'e and Material 
Constant A as a Function of Temperature 
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Figure 90. NARloy-Z Stress Rupture 
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Figure 91. EDNi Elastic Modulus as a Fucntion of Temperature 
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Figure 93. EDNi Thermal Expansion 
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Figure 96. EDNi Stress-Strain Diagram 
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Figure 98. EDCu Elastic Modulus 
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Figure 100. BDCu Thermal Expansion 
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Figure 101. EDCu Tensile Strength 
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Figure 102. EDCu Stress-Strain Diagram 
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Usually, when nominal materials properties are evaluated, the concepts of 
engineering stress (s) and engineering strain (e) are employed. However, 
since dimensions do change as the load is applied, calculations of engineering 
stress and strain are subject to errors. To provide a more accurate measure 
of stress and strain, the quantities true stress (o) and true strain (c) 
have been defined in the following way: True stress S is the actual stress 

based on the actual area A corresponding at every instant to the current 
value of load P. 

The relationship between true stress and strain and engineering stress and 
strain can be expressed as follows: 

o = S (1 + e) 

(valid up until necking) 

G = n (1 + e) 

Using the above relationship, a comparison of the true stress/strain curve 
versus the engineering stress/strain curve showed only a small difference for 
NARloy Z (Fig. 87). Engineering stress/strain curves were used for the 
structural analysis. 

3. Life Prediction 


a. Low Cycle Fatigue 

The low cycle fatigue life was predicted for NARloy-Z in the STA condi- 
tion at 978 K (1300 F) and 1089 K (1500 F) based on Hanson's Universal 
Slopes (MUS) equation. However, service at such a high temperature is 
well into the creep range of the material. Also, tensile properties for 
all the temperatures of interest are not available since 1500 F is 
close to the solution temperature of NARloy-Z. Therefore, extrapolated 
values were used. But even if tensile properties were available, ten- 
sile strength and ductility is very sensitive to strain rate in the 
creep range and their values would depend on the rate at which the 
strain is applied in the tensile test. Therefore, it is reasonable to 
assume that some time factor would have to be introduced into the anal- 
ysis. Hanson's Universal Slopes method does not take this into account. 

In general, attempts to predict high temperature fatigue behavior in 
the creep range by using tensile test properties at the temperature of 
interest result in predictions of unconservative lives that are 
obtained by fatigue test. Part of this discrepancy is due to the large 
creep component at high temperatures. As an attempt to account for 
creep-fatigue interactions, Hanson suggests a 10% rule at high tempera- 
tures (Ref. 17). LCF life is taken as 10% of the predicted value. But 
for NARloy-Z, HUS yields an unconservative slope which crosses test 
data. Available test data from 294 to 866 K (70 to 1100 F) show an 
increasing slope with temperature for a log percent total strain range 
versus log cycles to failure plot (Fig. 88). This is probably due to 
an increase in the amount of creep damage at the higher temperatures. 
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Another approach to high temperature life prediction is based on linear sum- 
mation of time dependent creep damage and cycle dependent fatigue damage, with 
failure anticipated when the total damage, equals a given value (in many cases 
unity for convenience (Ref* 18). This model assumes that fatigue and creep 
are linearly interactive* But for high temperature life prediction, the 

linear damage summation rule must be used carefully since it depends upon the 
amount of mixing of the damage components, as well as the order of damage for 
fit given material and temperature. Plumtree and Abdel-Raouf studied the 
effects of first imposing one type of damage (e.g., fatigue) followed by the 
other type of damage (e.g., creep-type) to failure (Ref. 19). It was found 
that the amount of damage which occurred during the first operative damage 
mechanism had no influence on the amount of damage accumulated in the 
following block if the two mechanisms were different. For OFHC copper at 
811 K (1000 F), linear summation to failure may vary from less than 0.4 for 
maximum interaction to greater than unity for little or no interaction of the 
damage mechanisms. 

Hirschberg-Halford-Manson predicted life during high temperature fatigue by 
partitioning the creep and plastic strains in each hysteresis loop and line- 
arly summing their respective damage contribution (Ref. 20). An attempt at 
using this method yielded appproximately the same results as MUS. Also, there 
was not enough data to predict all the damage fractions. But this, too, is an 
idealized method, because only damage associated with plastic strains are 
summed. 

Since none of the life prediction methodologies can accurately account for 
creep-fatigue interactions, the NARloy-Z prediction at 978 K (1300 F) and 
1089 K (1500 F) are merely a best-guess estimate taking into account the slope 
determined from test data. 

b. Creep 

Creep in many materials at intermediate stress levels can be described by 
a power law creep expression (Ref. 21 - 23). 

c = K exp (-Q /RT) (1) 

a fi^pp 

where Cg the steady state creep rate, K is a material constant, cf^ is 
the applied stress, n is the stress exponent for creep, Qa^ apparent 

activation energy for creep, R is the gas constant, and T is the absolute 
temperature. At a given temperature, the above equation can be reduced to the 
following assuming K to be independent of creep strain: 

Cg = A (2) 

where A is a constant at a given temperature 

. * . A = K exp. ^ (3) 

R 

A plot of log creep rate versus log stress will generate a straight line with 
the slope being the value of the stress exponent if Eq. 2 holds. 
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Values of the stress exponent n and the material constant A of NARloy Z for 
the temperature range 644 to 922 K (700 to 1200 F) were taken from Ref. 24. 
The values were based on test data. Values at 1089 K (1500 F) were estimated. 

Exponent n at 1089 K (1500 F) is based on extrapolation of the stress exponent 
versus temperature data given in Ref. 24. Extrapolation resulted in a value 
of 4.6 at 1089 K (1500 F) . 

The constant A at a given temperature can be calculated from the following 
relationship: 


A = K exp =1 

where , 


K, for NARloy Z ~ 1.6 x lO^^ 

Q, activation energy = 110 x 10^ cal/mole 

R, Universal Gas Constant = 1.98 cal/mole K 
T, Temperature = 1089 K (1500 F) 

Substitution of these values in the above equation yields the value of A as 
1.10 X 10“^. 
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APPENDIX C 


EXTRAPOLATION IMPLEMENTATION* 

The APSAC analysis system accounts for deformed geometry and inelasticity 
through a total Lagrangian approach that utilizes the geometry changes from 
the previous loadings and the historical information related to inelastic 
strains and yielded surface. The overall capability already existed for 
sequential analysis of a duty cycle through the use of a restart dataset, so 
the incorporation of the extrapolated geometry into the analysis was rather 
straightforward. The most important point was to define the overall extra- 
polation requirements and determine how to best implement them into the avail- 
able programs. The requirements for this task included: 

1. Being able to output a displacement dataset at the end of specified 
duty cycles 

2. Combining and merging separate displacement datasets for use in the 
extrapolation 

3. Extrapolation from a series of duty cycles to a specified future duty 
cycle 

A. Plotting extrapolated geometry 

5. Generation of an equivalent dataset for use with the APSAC code to 
restart the next analysis sequence 

After reviewing the available codes and how they interact, the block diagram 
strategy was developed for the extrapolation, as shown in Fig. 104. The top 
half of the page outlines the current applicable codes and their interaction. 
Basically, a dataset was output from APSAC that included all information for 
future identification of the job (input dataset), all necessary arrays for 
restarting the job for the next load increment, and output results for plot- 
ting. This dataset also has other uses such as the input to the STAMP code to 
generate geometric data for the DEAP thermal code. This methodology was used 
for this analysis so that the thermal and structural models were identical, 
and the temperatures could be properly accessed from the DEAP datatape. 

The output dataset is currently accessed by APSAC for two separate purposes: 
one (Tape 11) is for restarting a job for additional loading increments or 
changing the geometry, and the other (Tape 13) is for reading in a field of 
displacements that can be used as fixed boundary displacements. This dataset 
is also used in the plot code as part of the displacement accumulation. 

The simplest method for accomplishing the extrapolation technique was to 
incorporate it as an option of the plot program. This code already read and 
searched the dataset for displacement information, and would be readily 
changed to output displacements and generate a modified dataset for plotting 
and for generating a restart dataset with the proper displacements. 


*This effort was funded by Rocketdyne. 
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FOR 

ACCUMULATING 

DISPLACEMENTS 


Figure 104. Block Diagram of Strategy for Incorporating 
Extrapolation Technique 
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Figure 105 shows the concept of the duty cycle analysis for five increments, 
an extrapolation and the restart for the next series of increments* Two 
options are listed, but only option 1 has been used for the contract analy- 
sis. The accumulated displacements from the extrapolation tape were added to 
the initial geometry, and the stresses and strains on the part were initiali- 
zed to zero. Essentially, a new analysis with a distorted geometry and no 
strain history. 

Option 2 is for experimental use to evaluate whether a forced distorted shape 
with or without the history of strains from the last increment prior to the 
extrapolation can be used to minimize or eliminate the restart shakedown 
observed in the analysis. A future option could be to add extrapolation of 
inelastic strains, etc., to get a more accurate restart condition. 

The extrapolation code evolved through the addition of several options to 
allow flexibility in positioning the model in the appropriate radial displaced 
position. Working through the first model from 1 to 120 duty cycles identi- 
fied the need for this flexibility to account for the initial liner gap and 
the liner movement away from the jacket during multiple duty cycles. These 
considerations were required since the liner is not attached to the jacket, 
except at its ends. At the start of the analysis, there is an initial 0.020 
tolerance gap. During the first duty cycle, the liner cross-section yields 
and radially displaces to meet the jacket motion. At the end of the first 
duty cycle, the liner is firmly compressed against the jacket. During the 
succeeding cycles, the liner has additional relaxation that reduces the 
compression against the jacket. This was essentially a small rigid body 
motion of the total model. Even though this motion is small, it was signifi- 
cant compared to the displacements of the rest of the geometry. During the 
extrapolation process, both this rigid body motion, as well as local distor- 
tions, were extrapolated and pulled the liner away from the jacket. The 
actual hardware, though, only relaxes the liner to some stable point — where 
the gross yielding of the liner is balanced by the relaxation effects. The 
actual liner motion is nonlinear, whereas the basic extrapolation technique is 
linear. Options have been added to the extrapolation code to subtract out 
rigid body motion, move the model into position, and adjust fixed boundary 
node positions to properly adjust for the true hardware motions. Essentially 
what was done was to allow the liner to yield up against the jacket in the 
first 1 to 5 duty cycles. Then, as part of the extrapolation technique, a 
rigid body motion was added to the model, such that one of the nodes in the 
liner that touches the gap element was held against the jacket. The gap ele- 
ment node positions were also adjusted to a compatible position. These dis- 
placements were then used to update the geometry for the 21st duty cycle. The 
same procedure was used as part of the next extrapolation with the addition of 
the accumulation of the current extrapolated displacements and the displace- 
ments that were used to start the 21st cycle. The latter technique is repeti- 
tively applied for the remainder of the analysis to develop the distorted 
geometry. This procedure is pictorially shown in Fig. 106. 
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OPTION 1 :• NEW GEOMETRY ~ NO HISTORY (VIRGIN PART) 

• NEW ANALYSIS, SAME AS 1-6 INCREMENT 

OPTION 2: •UPDATED GEOMETRY 

HISTORY THROUGH LAST CALCULATED CYCLE 

• SPECIAL RESTART CYCLE 

• EVERY NODE HAS FORCED BOUNDARY 
CONDITION OF DISPLACEMENTS IN 

2 DIRECTIONS 

• UTI LIZE EXTRAPOLATED GEOMETRY 
AND HISTORY INFORMATION 

• NO OTHER APPLIED LOADING 


Figure 105. Duty Cycle Analysis Concept 
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0.508 mm 
0.020 MOTION 


RADIAL 



GAP 

ELEMENTS 


INITIAL 

CONDITIONS 



END OF FIRST 
CYCLE 



END OF FOURTH 
CYCLE 


0.020 IN. 



CYCLE 



NODE-B 


EXTRAPOLATED 
GEOMETRY AT 
END OF 20TH 
CYCLE 


RESET POSITION 
OF GAP ELEMENT NODES 




0.508 mm 
0.020 IN. 


• MOVE LINER THROUGH RIGID BODY MOTION SO 
THAT NODE-B = 0.020 RADIALLY 

• ADJUST GAP ELEMENT BOUNDARY TO 0.020 TO 
WHERE A IS THE LENGTHOF THE GAP ELEMENTS 


ON FUTURE RESTARTS 
(41, 61, ETC.) ACCUMULATE 
DISPLACEMENTS FROM PREVIOUS 
EXTRAPOLATION 


Figure 106. Extrapolation Detail Model Runs 







pe plot program input revisions to add the extrapolation 
included in the following pages. These can be run as a separate 
mixed with a standard plot data deck. For traceability, it is i 
the extrapolation as a separate computer job. 

The details of the data card inputs follow. 


technique are 
job or inter- 
simpler to run 
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CARD 1: 


TITLE CARD 


NO. REQUIRED: 


FUNCTION: 

FORMAT: 


One card per load increment of APSAC analysis. 
Title and terminator card must be present whether 
or not increment is considered for extrapolation 
or plotting. 

Title used on each CRT within the current load 
increment 

18A4 


— 




1 


— * — ' — ' — ' — 1 — 1 — ■- ■ ■ . 


13 


— ' — • — « — I — 1 — . ■ . 


25 


— - . .T. ijl.L . E . 


37 


' * ' 1 • ■ 


49 

1 


' ' 1 J 1 . 

— 


[ — 

1 • 

' ' ' 1 L__l j , 5 


identification . 73 ] .A 


CARD 1; TITLE CARD 


lELi" be present on the title card for the first 

increment (first card in data deck) to write displacements on Tape 4 
for later use in extrapolation procedure. 
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EXTRAPOLATION 


CARD 2: 


CONTROL CARD 


NO. REQUIRED: 
FORMAT : 


One card for each increment from which 
displacements are to be written to TapeM 

(in, 211, 12, 13, 4E6.0, 9A4) 


! 

Column 

! Format 

I 

Variable 

! 

Description 

! 

! 

! 

1 

1 

B 

1 

14 

I 

! 

! 

LDEL 

! 

! U 7 Ref. 
I G-1. 

node as described in Table 

! 

! 

1 

! 

1 

5 


! 

! 

11 

! 

! 

IGRID 

! 

! IGRID=5 

- Displacements for current 

I 

! 

! 

! 



! 

! 


1 

1 


! 

! 

increment are written on 
Tape4. 

! 

! 

! 



! 


1 


! =6 

- Extrapolate displacements 

! 

! 

! 



! 

I 


1 

! 


! 

! 

to Nth increment, using 
this and the previous 

! 

1 

! 



! 


1 


I 

increments from Tape4 

i 

1 



1 


! 


1 =7 

- Extrapolate to Nth incre- 

1 

1 



! 


! 


! 

ment using this and the 

t 

! 



! 


! 


I 

previous increments from 

! 

! 



! 


! 


! 

Tape4. Subtract the rigid 

j 

! 



1 


! 


! 

body displacement 

! 

! 


IGRIE 

= 

6, 7, 

or 8, read 

1 =8 

- Extrapolate to Nth 

j 

1 


in card 

2 



1 

increment using this and 

! 

1 



! 


! 


! 

the previous increments 

! 

1 



! 


1 


! 

from Tape4. 

! 

I 



! 


1 


1 

Options include: 

! 

! 



I 


I 


! 

1)node no. to keep to a 

! 

! 

! 



1 

1 


! 

! 


1 

! 

specified displacement 
(DISP) 

! 

! 

! 



1 


! 


1 

2}nodal displacements to 

! 

1 



! 


! 


! 

directly replace the 

! 

1 



! 


! 


! 

extrapolated value. 

! 

1 

1 

6 


! 

1 

11 

! 

! 

ITYPE 

! 

1 ITyPE=value=order of extrapolation 

! 

! 

1 



! 


! 


1 for IGRID=6, 7. or 8 

! 

I 

! 

7 

- 8 

! 

! 

12 

! 

1 

IFUN 

1 

! 

Not Used 

! 

I 

1 

! 

9 

- 12 

1 

! 

14 

1 

I 

NNC 

1 

1 NNC=value=increment identification 

! 

! 

1 



1 


1 


I for Tape4 

! 

! 

! 

li 

- 36 

I 

1 


1 

1 


! 

! 

Not Used 

1 

! 

1 

! 

2L 

- 72 

1 

I 

9A4 

1 

1 

LABEL 

! 

! 

Title Printed 

! 

! 
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CARD 3: 


EXTRAPOLATED CYCLE CARD 


FUNCTION: Immediately follows any extrapolation card where 

IGRID=6, 7» or 8, to define the number of the cycle 
where extrapolated data is desired. 


d 


N.I .N .C .R 1 .M.I . N.C R 


In 


. N.0_D. E .N.D; I.S.P 


E 


. .U.X.D I I S.P . . 


E 


. .U.X D ,r.S .P . _ 


E 


..F.L.A.G 


E 


1 I 1 1 1 1 « « « A 


[fOENTIFICATIO* T3 ^0 


DISPLACEMENT data from the current increment and displacement data 
from TAPE4, either from previous increments of same job, or saved 
from a previously run job, and the selected order of extrapolation 
are used to calculate the displacements at "NINCR." These may be 
saved on ''TAPES" and used for restart purposes. "TAPES" has the 
same format as "TAPEIO," generated by APSA or APSAC. TAPES then has 
accumulated displacements up to the duty cycle at the end of the 
current extrapolation. Only 1 increment is generated on TAPES*. 

When accessing this tape for restart, request the first increment. 


•TAPES has the stress and strains from the local increment from 
wnich the extrapolation is generated, plus the extrapolated 
displacements . 
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REQUIRED 

FOR 

IGRID 



180 


CARD 3 


! 

VARIABLE ! DESCRIPTION 

! 

! 

NINCR lINCREMENT WHERE EXTRAPOL- 
lATED DATA IS DESIRED 

! 

! 

MINOR lINCREMENT TO BE USED AS 
! REFERENCE ZERO FOR RIGID 
IBODY DISPLACEMENT REMOVAL 
!(MINCR=0, USE AVERAGE OF 
lALL PREVIOUS INCREMENTS) 

I 

! 

NODE !IGRID=7, NODE NUMBER TO BE 
lUSED AS REFERENCE ZERO FOR 
IRIGID BODY DISPLACEMENT 
I REMOVAL, IGRID=8, NODE 
INUMBER TO BE KEPT TO A 
ISPECIFIED DISPLACEMENT 
KUXDISP, UYDISP) 

! 

. j 

NDISP INUMBER OF DISPLACEMENTS TO 

IBE READ IN AND DIRECTLY 
I REPLACE THE EXTRAPOLATED 
IVALUES (SEE NEXT CARD) 

I 

I 

IFIXED DISPLACEMENT 
UXDISP IVALUES WHICH WILL BE 

I ASSIGNED TO NODE "NODE” 

I CAUSING THE REST OF I 

UYDISP ITHE DISPLACEMENTS TO BE I 

I TRANSLATED BEFORE I 

I EXTRAPOLATION I 

I ! 

I 

IFLAG lIFLAGrI, READ IN 

I DISPLACEMENTS FROM TAPE? : 
I (OLD TAPES) AND ADD TO. 

I CURRENT EMRAPOLATED 
I DISPLACEMENTS 
I 



CARD 4 


DISPLACEMENT CARD 


NUMBER REQUIRED: "NDISP" 

FUNCTION: To replace displacement values after 

extrapolation with the specified input value 



CARD 4 


II! ! 

! COLUMN I FORMAT ! VARIABLE ! DESCRIPTION 

I ! ! I 



I I 
NODE INODE NUMBER FOR REPLACING THE DISPLACE! 

IMENT UX (OR UR) AND UV (OR UZ) AFTER ! 
I EXTRAPOLATION ! 
I ! 








Increment Termination 


Card b: 

NO. REQUIRED; 
FUN CTION; 


Increment End Card for APSAC 

One to end each load increment 

Terminates a load increment. Terminates the 
repeat plot card set. Specifies that the 
repeat plot set is to be used in this load 
increnent 


I Column I Format I Variable 1 


Descri ption 



t ! 

! fa - 3b 1 
! ! 


48 1 A12 




! ! 

IGRID ! IGRID = 2 Terminate the load increment 

! ! 


ITYPE 

1 

! 

ITYPE Non 

Blank 

! 

, Plot the repeat ! 


! 

1 

dataset in the 

current increment ! 

! 


! 

1 

i 

Not used 



1 

! 

1 

"END PLOT 

SET" 

specified on the end ! 

1 

1 

card used 

to terminate the repeat 

1 

1 

! 

1 

plot card 

data 

set 

! 

! 

1 

! 

1 

1 

Not used 




INCREMENT END CARD FOR APSAC 
IGRID ITYPE 


NUMBER DESCRIPTION 


IDEMTIFICAHON 

















APPENDIX D 


FINITE ELEMENT MODEL 


INTRODUCTION 

The APSAC — Axi symmetric and Planar Structural Analysis (APSAC) with Creep code 
and associated system of programs were utilized for this study. Ths code has 
the capability (options) for accounting for deformed geometry and inelasticity 
through use of a total Lagrangian approach that utilizes the geometric changes 
through a duty cycle, as well as a kinematic strain hardening flow rule for 
tracking the stress and strain history of the analysis. The code is part of 
an integrated set of codes for structural and thermal analyses that allowed 
optimization of the entire analysis. This was an important factor in being 
able to complete the study within cost and time constraints. The total analy- 
sis requires a large number of individual model runs, plus extrapolations. In 
order to accomplish these objectives, the overall study was reviewed, and a 
strategy developed to efficiently set up, track, and process the models. The 
basic philosophy was to organize the study so that it could be used, not only 
for this effort, but for future use on advance design or production analysis. 
This appendix discusses the details of the finite element model (FEM) and the 
strategy for processing the jobs. The thermal model is closely related to the 
structural model in that the same finite element model is used to develop the 
nodes, connectors, volumes, etc. for the thermal analysis. This results in an 
efficient one-to-one correspondence between nodes of the two models, and 
allows the structural model to assign temperatures to nodes from specified 
time slices in the thermal duty-cycle analysis. The one-to-one node match was 
used in an APSAC code modification to allow the thermal field to track the 
structural model as it distorts. 

Structural Model 


The primary structural model of concern is the 2D generalized plane strain 
model of the coolant liner. The jacket model will not be discussed in detail, 
since it is only used to furnish boundary displacements and the out-of-plane 
total strain for the various duty-cycle increments. The basic model and anal- 
ysis considerations are summarized in Fig. 107. The FEM model is one-half of 
a channel cross-section (Fig. 108). The model has symmetry boundaries on both 
circumferential edges. The basic model subtends an angle of 0.46154 degrees, 
since it is one-half of one of the 390 channels that comprise the 360 degrees 
total circumference (1/2 x 360/390). The inner edge (inside diameter) of the 
model is pressurized at chamber pressure, P2, and the coolant channel is pres- 
surized to hydrogen coolant pressure, PI. The pressures are consistent with a 
location that is 1.2 inch forward of the MCC throat. Pressure is denoted by a 
V on the model pointing in the direction of positive pressure. The B and E 
signify the sliding boundary condition definition, and the C signifies the 
structural jacket boundary condition. 

The liner has three materials — NARloy-Z (1), EDNi (2), and EDCu (3), — note the 
( ) values are the material nomenclature used in the code. A fourth material 
(4) is used to specify the gap element identified as IS0=4 material in the 
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code. This material is stiff in compression but weak, in tension and shear, 
which allows proper simulation of the liner-jacket interface. The model node 
and element numbers and material numbers are shown in Fig. 109. The materials 
used are isotropic, with bilinear stress-strain curve representation, which 
require temperature-dependent properties over the range of temperatures the 
model sees throughout the duty cycle. Each material has up to 8 sets of input 
specifying temperature, modulus, Poisson's ratio, coefficient of expansion, 
yield stress and tangent modulus ratio. The temperature increments are sel- 
ected to best approximate the total set of properties that will be linearly 
interpretted for actual values of elements at intermediate temperatures. The 
bilinear stress-strain curves used are furnished in Fig. 110 and 111. 

Where applicable, secondary creep constants are input for the materials. Only 
NARloy-Z has temperatures above the creep activation temperature, and has 
creep constants specified. Again, the creep properties are input at discrete 
temperatures to properly control the temperature interpolation. The material 
constants and creep relaxation data are obtained from the material curves 
included in this study. Table 9 shows the details in setting up the NARloy-Z 
creep constants in the format used by the code. 

The rest of the model geometry definition consists of defining an (I, J) 
topology numbering system to relate the model grid and coordinates. Input to 
define materials and boundary conditions are then specified, and the code 
develops the node and element sequencing for the above information. Other 
inputs include the out-of-plane strain, displacements, pressures, and tempera- 
tures. For all runs, the temperatures are recovered from a thermal DEAF code 
datatape for every time slice, except the room temperature posttest increment 
11. Each analysis run consisted of 55 increments used to define the 5 duty 
cycles (11 increments per duty cycle). Table 10 summarizes the duty cycle 
loading on the structure, except for thermal conditions. 

Title and Control cards are input to specify analysis options. For these 
runs, the large deflection codes to control the analysis are specified as 
shown in Table 11. The codes define how the large deflection option is to be 
processed on each increment of load, and for restarting from a datatape. The 
creep increments are handled slightly different from the load-only increments, 
due to a code limitation that both creep and large deflections are not allowed 
within the time-dependent portion of the total increment. This is not a sig- 
nificant limitation in that the geometry is updated at the end of the creep 
increments — 5 and 6. In the other loading increments, the geometry is updated 
throughout the calculation process. 

The model was optimized and thoroughly checked as part of the initial study. 
This optimization included varying the model refinement to achieve acceptable 
accuracy in stresses and strains, as well as temperatures. The model, as was 
shown in Fig. 108 has variable size elements in the x-direction to pick up 
representative surface strains and strains in the channel wall to land-notch 
area. The same type variation was included in the EOCU notch zone to properly 
represent the inelasticity that occurs during the duty cycle. The model was 
exercised through a series of small deflection duty cycles with and without 
creep, and large deflection duty cycles with and without creep, to build up 
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Figure 109. Nodal Network and Materials Identification 
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TABLE 9. STRESS RELAXATION (CREEP) CONSTANTS 




E 

YOUNG'S MODULUS 

n 


■ 

TEMPERATURE 


MATERIAL 


A 


K 

F 

MPTR 

70-902 

PROPERTIES 

MANUAL 


A 

MPTR70-902 

A 

ADJUSTED BY E* 

22 

-420 


20.0 X 10^ 

19.2 

1.0 X 10“®° 

1.0 X 10“®° 

644 

700. 


16.2 X 10^ 

16.1 

2.9 X 10"^’ ' 

2.9 X 10“^’ 

811 

1000. 

22.9 X 10^ 

14.9 X 10^ 

12.5 

6.2 X 10"^® 

(6.2 X 10“^®) = 9.529 X 10-^® 

922. 

1200. 

15.7 X 10^ 

14.0 X 10® 

8.7 

4.5 X 10“^ 

14.0 ^ ) = 5.046 X 10 ^ 

1089. 

1500. 

10.0 X 10^ 

11.5 X 10® 

5.08 

2.073 X 10“^ 

^^5 (2.073 X 10 ^) = 1.803 x 10“^ 

1 

THE STRESS RELAXATION FORMULA = [AE(n-l)t + p"'’ SHOWS THE DEPENDENCY OF THE CONSTANT, A, 

UPON THE VALUE USED FOR YOUNG’S MODULUS. FOR THIS ANALYSIS, A WAS ADJUSTED TO REFLECT THE TRUE VALUE OF E. 



TABLE 10a. SSME DUTY CYCLE INCREMENTS THERMAL, PRESSURE 
BOUNDARY CONDITIONS DATA FOR ANALYSIS 
(ENGLISH) 


LOAD 

INCREMENT 

DEAP TIME 
SLICE NO. 

TIME, 

SECONDS 

^HGW’ 

F 

^OACK* 

F 

•"hot WALL’ 
PSIA 

^^COOLANT’ 

WALL 

PSIA 

JACKET BOUNDARY 
CONDITIONS 

^R’ 

INCH 

®Z’ 

INCH/INCH 

1 

15 

1.4 

-365 

R.T. 

40 

525 

0.00016 

-0.000013 

2 

19 

1.6 

404 

R.T. 

384 

909 

0.00149 

-0.000122 

3 

35 

2.4 

610 

R.T. 

518 

2150 

0.00202 

-0.000165 

4 

55 

3.4 

1015 

-23 

1436 

4035 

0.00310 

-0.000404 

5 

58 

4.9 

1306 

-100 

2109 

6270 

0.00370 

-0.000578 

6 

60 

20.5 

1324 

-175 

2109 

6270 

0.00175 

-0.000540 

7 

79 

395.0* 

1047 

-175 

1267 

3629 

-0.00149 

-0.000275 

8 

85 

395.6 

739 

-100 

537 

1930 

-0.00234 

-0.000084 

9 

100 

397.1 

-297 

-84 

109 

1028 

-0.00354 

0.0000403 

10 

129 

400.0 

-405 

-82 

9 

0 

-0.00386 

0.0000730 

n 



70 

70 

0 

0 

0 

0 


*DEAP DATATAPE TEMPERATURE INPUT ONLY CHANGE IN INPUT LOADS FROM MODEL TO MODEL 



TABLE 10b. SSME DUTY CYCLE INCREMENTS THERMAL, PRESSURE 
BOUNDARY CONDITIONS DATA FOR ANALYSIS 
(METRIC) 


LOAD 

INCREMENT 

DEAP TIME 
SLICE NO. 

TIME, 

SECONDS 

^HGW* 

K 

^JACK' 

K 

^HOT WALL’ 
MPa 

^COOLANT’ 

WALL 

MPa 

JACKET BOUNDARY 
CONDITIONS 

"r* 

(mm) 

cz. 

(mm) 

1 

15 

1.4 

53 

R.T. 

0.2758 

3.6197 

0.004064 

-0.00033 

2 

19 

1.6 

480 

R.T. 

2.6476 

6.267 

-0.037846 

-0.003099 

3 

35 

2.4 

594 

R.T. 

3.5646 

14.824 

0.051308 

-0.004191 

4 

55 

3.4 

546 

243 

9.9009 

27.820 

0.07874 

-0.010262 

5 

58 

4.9 

819 

200 

14.541 

43.230 

0.09398 

-0.014681 

6 

60 

20.5 

991 

158 

14.541 

43.230 

0.04445 

-0.013716 

7 

79 

395.0* 

837 

158 

8.7356 

25.021 

-0.037846 

-0.006985 

8 

85 

395.6 

666 

200 

3.7024 

13.307 

-0.059436 

-0.002134 

9 

100 

397.1 

90 

209 

0.7515 

7.0878 

-0.089916 

0.001024 

10 

129 

400.0 

30 

210 

0.0621 

0 

0.98044 

0.001854 

11 



R.T. 

R.T. 

0 

0 

0 

0 


*DEAP OATATAPE TEMPERATURE INPUT ONLY CHANGE IN INPUT LOADS FROM MODEL TO MODEL 




TABLE 11. SUMMARY OF APSAC CONTROL CODES FOR MCC 
LARGE DEFLECTION ANALYSIS 



CYCLE 

INCREMENT 

NEWRZ 

LGDFL 

AO 

TOLLD 

LDOPT 

LDRST 

CREEP 

Initial 5-Cycle 

1 

1 


-1 

.7 

.0005 

2 

1 


Analysis 


2 



.7 

.0005 

2 

3 




3 



.5 

.0005 

2 

3 




4 



.7 

.0005 

2 

3 




5. 6 





1 

2 

1 



7-n 



.7 

.0005 

2 

3 



2 

1 



.7 

.0005 

2 

3 




2 -n 


Same as Cycle 1 




3-5 

Repeat Cycle 2 

Restart for 

21 

1 

1 

-1 

.7 

.0005 

2 

1 


Succeeding 5- Duty 
Cycl es 

Etc. 



APSAC 


LGDFL / 0 Large deflection analysis for some increment of this job. Read 
in additional card to define type of analysis (to be performed 
this iteration) and start geometry. 

< 0 Change temperature coordinates read from DEAR datatape for this 
increment to nodal coordinates for this increment. There must 
be an exact correlation between nodes. 


NEWRZ > 0 

Displacements from a previous duty cycle are added to basic 
nodal geometry. Used to restart after first 5 cycles and 
extrapolation to cycle 20. 

AO 

Factor to limit amount of geometry change per iteration in 
large deflection analysis. 

TOLLD 

Large deflection iteration convergence number 

LDOPT 

Type of analysis to be performed this increment 
=1 elastic-plastic 

=2 elastic-plastic and large deflection simultaneously 
=3 elastic and large deflection simultaneously 

LDRST 

Geometry to start this increment 
=1 start geometry is basic geometry 

=2 start geometry is basic geometry + U from previous increment 

=3 start geometry is basic geometry + U from previous increment, 
U* equals U' from previous increment 

«4 start geometry is basic geometry + U from previous increment, 
U' equals U' from previous increment 

Previous increment is stored on tape 11, increment number = 
INREADI 




overall confidence in the total methodology. The initial duty cycle had been 
divided into 13 increments, but two increments were eliminated when model 
testing deomonstrated that the two additional increments in the steady— state 
portion of the cycle were not significant in defining the stress and strain 
and displaced conditions. Elimination of these increments reduced computer 
run costs by about 10%. 

This model checkout also showed that the compression only elements must be 
carefully sized when used with this type of large deflection analysis. The 
APSAC code performs an iteration process, and the geometry distortion is 
accounted for during these iterations. If too large a load increment is used 
as the liner displaces toward the jacket boundary, the model can overshoot 
during the iteration; and the boundary elements can turn inside-out, which is 
not allowable in the code. To minimize this problem, the boundary elements 
were sized to be almost square. Usually for small displacement analysis, thin 
elements are used; but the size is inconsequential to the analysis, except for 
a code restraint that the longest side must be normal ho the stiff direction 
of the element. The problem of too large a displacement within a load incre- 
ment, as the liner deformed through the 0.508 mm (0.020 inch) initial gap and 
touched the jacket, occurred in two of the models. This was readily elimin- 
ated through adjusting the displacement relaxation factor, AO. Input data- 
card^ image listing of the 1—5, 21—25, and 41—45 duty cycles and extrapolations 
are included for examples of the procedure used to run all the models. Other 
geometries are only minor variations of these runs. These listings show all 
the detailed code input required, and the methodology used for control, track- 
ing and changing the input from the terminal. 

Model Stragegy 

The model strategy used for this study was to standardize the input so that 
only a minimum number of changes were required to submit individual jobs and 
to develop a naming convention for datasets saved and processed in the analy- 
sis sequence. The other major technique was to utilize the terminal and edi- 
tor for making all job sulmittals and input data changes. This procedure 
effectively eliminated job submittal errors and reduced the time to review and 
process a 5 duty-cycle analysis to less than 4 manhours effort. 

The decks were setup so that only two job-control cards, the title cards for 
each increment and one or two other cards, need modification. The changes on 
the cards are also standardized so that all changes to the data deck can be 
made with a few (less than 10) editor commands. This is accomplished with the 
naming convention of datasets and identification on title cards. The dataset- 
naming conventions are as shown in Table 12. For instance, using the editor 
feature of replacing character sets, by replacing AOOl by A021 changes all the 
dataset names in a submittal from the first-cycle input convention to the 
21st-cycle input convention. Since the rest of the duty-cycle definition and 
the temperature dataset did not change, the only other change dealt with 
adjusting large deflection control flags from an initial submittal to a con- 
tinuing analysis, and adjusting the displacements in the 21st cycle load 
increments to account for the rigid body motion of the structure of 0.508 mm 
(0.020 inch) that occurs between the first and 21st duty cycles. 
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TABLE 12. DATASET NAMING CONVENTIONS USED FOR COMPUTER RUNS 


JOB NAME 

TITLE 

PERMANENT 
FILE NAME 

FILE USAGE 

QMCCA 

AOONOX 

TAOONOX 

TAPEIO 



PAOONOX 

PRINT 



CAOONOX 

CRTS 


WHERE A IS THE CODE DESCRIBING THE TYPE OF ANALYSIS BEING RUN - 

• BASIC SSME ANALYSIS lOOA HG CASE 

• BASIC SSME ANALYSIS 135A HG CASE 

• CONTOURED GEOMETRY 

• SLOTTED 

• KFEL RIB 

• 540 CHANNELS 


N IS THE INITIAL CYCLE OF EACH ANALYSIS, AND X IS THE INCRE- 
MENT WITHIN THAT CYCLE CYCLE THAT THE JOB IS INITIATED. THE 
RESTART, PRINT, AND CRT FILES ARE CODED CONSISTENTLY TO 
SIMPLIFY JOB TRACKING. 
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For the standard analysis, all results were obtained from the postprocessing 
plotting code, either as in graphical or pictorial format, or as printouts of 
selected displacements, stresses or strains. Other detailed information from 
job runs were scanned at the terminal from saved datasets. 

This overall strategy ended up saving considerable costs from computing, man- 
hours, and overall lapsed time. As previously mentioned, it also virtually 
eliminated job reruns from data errors. 

The remaining pages of this section furnish examples of the card images for 
running typical APSAC and APSAPLOT/EXTRAPOLATION codes. The first set of card 
images show the setup for the first 5 duty cycle analysis that accounts for 
the initial liner to structural jacket gap. The next set of card images show 
the setup of the analysis runs where the liner is in contact with the jacket 
at the start of the duty cycle. The extrapolation card images for the first 
and second extrapolation cores are next. This set of card images was used to 
obtain plots and output summaries. 
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2 

THIDlOf BCC LINER t.39 HQ CYCLEIOI ♦3.INCR-06 T. 20.5 4 

I * 

.THtOtOlMCC LINER 1 .39 HO CYCLEIOI »3.INCR>09 T«399. S 

<2 

TH1010I BCC LINER 1.35 no CTCLEtOI ♦3.INCR>10 T«395.S 1 

2 

TH10101 BCC LINER 1 .39 HO CYCLEIOI -fa.INCRall T-397.1 $ 

a 

:THt010I BCC LT11ER'f:3T Ha CTCLC10!^^3.INCRi12 T.400, c 

2 

tHIOIOl BCC LINER 1 .39 HO CYCLEIOI ♦3.1NCRa13 T- INF. « 

1 2 4 

JS”— ■ - - 

1 1 9 9 -1 72 0 94 

9 4 

1 

tmOIOt BCC LtNER t. 35 no CYCLEIOI ♦4,INCR*01 T« ' T.4 % 

2 

jTHtOlOl mrc LINER 1.39 HQ CYCLEIOI 4-4.INCR>02 T* 1.9 S 

i * 

ITH10101 BCC LINEN t.39 nO CYCLEIOI 44,TNCR*03 T* 2.4 5 

I a 

{YH10101 BCC LINER 1.39 HO CYCLEIOI «4.lNCR*04 T- 3.4 $ 

I a 

I THIOIOt BCC LINEN t.39 HO CYCLEIOI ♦4.lwm«09 Ti 4.9 S 

2 

ITHIOIOI BCC LINER 1.39 HO CYCLEIOI 44.INCR>09 T> 20.5 $ 

2 

itHIOlbt BCC LINEN t.39 HO CYCLEIOI 44,INCN«09 T*399. S 

2 

THIOIOt BCC LINEN t.39 HO CYCLEIOI ♦4,IMCR«10 Ta395.6 $ 

2 

'tHIOIOt BcC IIHEN t.39 HO CTCLEIOt 44,INCN«11 T«39T.T S 

2 
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t 

9 

9 

• 

f 





--- 38 



i 

130 0 

227 


1 

4 4 




» 

I 


f 

1 8 

8 

- 

9 

1 8.1 8 

.40 

0.0 


>2 *18 8.2 9 

.40 

0.0 


23 *18 9.2 9 

.40 

0.0 

81 

9 ^ 



20 

217 

fO 

.00 

>18 

.0000 


0.0 

218 

.0000 


0.0 

220 

.0000 


.. 0.0 

221 

.COOP ' 


0.0 

222 

.000 0 


0.0 

223 

.000 0 


0.0 

224 

.000 0 


0.0 

229 

.000 0 


0.0 

228 

.000 0 


0.0 

»T 

.000 0 


0.0 


2 

SB 

0.20 

0.20 

0.20 


39 



220 


OPLTS, T7( . 

USER(0 - I 

CHARGE ( * ) 

ACCT(HAVEN . 000454St27*055' 

PURGE, PRD101/NA. 

RETURN, PRDIOi. 

DEFINE, PR0101. 
GET,PROCPMR/UNaQOB28. 
attach .TAPE 1 0«TD10 tOt . 

BEGIN, APSAC, PROCPMR.NAPSAC. 
rewind , OUTPUT. 

COPYEI . OUTPUT, PR01 01. 

COPYEl , OUTPUT, PRD1 01. 


PLOTMCC 


1 .35HG 


CARD IMAGES FOR APSAPLOT FOR PLOTS AND OUTPUT SUMMARIES 



21 

010 

21 

010 

21 

DIG 

21 

010 

3 

3 


1 

2 

3 

A 

22 

It 


21 

1401 

1402 

1403 
14 4 

1422 
14 1 1 
14 12 

1413 

1414 

1423 
1421 

1 5 

8 

8 

8 

I 5 

II 


MuC 

3 4 

LINER 

1.35HG 

CYCLEIOI 

•INCR*01 T« 

1.4 

TMAX.-360 

$ 

22 


1 


2 

3 


4 


1 

3 

8 

8 

-t 

56 

0 49 




23 


11 



13 


14 

0 






END >>LOT SET 




MCC 

LINER 

t .35HG 

. CYCLEIOI 

,INCR*02 T« 

1.6 

TMAX> 488 

$ 

MCC 

LINER 

1 .35HG 

CYCLEIOI 

,INCR«03 T» 

2.4 

TMAX« 718 

$ 

MCC 

LINER 

1 .35HG 

CYCLEIOI 

»INCR^04 Ta 

3.4 

TMAX.1157 

s 

MCC 

LINER 

1 .35HG 

CYCLElOt 

tiNCRaOS T« 

4.9 

TMAX»1465 

$ 

MCC 

liner 

1.35HG 

CYCLEIOI 

tINCR*06 T* 

20.5 ’ 

TMAX.1484 

$ 

1 

1 

5*2 

5.40 

0.0 

0.20 


‘ ' 



-18 

5.2 

5.4 

0.0 

0.20 





1 -18 

5.2 

5.4 

0.0 

0.20 

• 




( -18 

5.2 

5.4 

0.0 

0.20 





1 -18 

5.2 

5.4 

0.0 

0.20 





! -18 

5.2 

5.4 

0.0 

0.20 





-18 

5.2 

5.4 

0.0 

0.20 





» -18 

5.2 

5.4 

0.0 

0.20 





1 -18 

5.2 

5.4 

0.0 

0.20 





-18 

5.2 

5.4 

0.0 

0.20 





t -18 

5.2 

5.4 

0.0 . 

0.20 





-16 

5.2 

5.4 

0.0 

0.20 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

,06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

-.02 

.06 






5.24 

5.32 

--02 . 

.06 






5.24 

5.32 

-.02 

.06 





16 









1 


• : 2 


3 

4 


22 


1 

0 

49 

a 

2 

0 50 

8 

3 0 

51 

. 4 

0 

52 

0 

5 

0 53 

8 

6 0 

54 

7 ' 

. Q 

55 

8 

8 

0 56 




11 

4 


. 12 


13 

14 


23 

0 


1 


34 
1 

35 


-1 


72 


64 


ttN£R 1.35HG 


21 

TD1010 1 MCC LINER 1 .35HG 
21 

TDIOlOl MCC 
21 

TD1010 1 MCC LINER 1 .35HG 
21 

TD10101 MCC LINEH 1*35HG.- 
1 5 16 

1 3 


CYCLEIOI ,INCR«09 Tr395. 

CVCLE101 -►ItlMCRslO Ta395«6 

CYCLEIOI ♦INCRall T*397.1 

CYCLEIOI ^lMai*t2 T 400* - 

3 4 


TMAX;. 1191 
TMAXs 858 
. TMAX*-297 . 
-.TMAXs-405- 
22 


S 

$ 

S 


222 



8 

8 

8 

1 5 
1 1 


1 

t 

21 


4 

7 

11 

4 

34 
1 

35 


Q 

0 

0 


49 

52 

55 

12 


8 

8 

8 


-1 


2 

5 

6 

13 


72 


0 

0 

0 


50 

53 

5S 

14 


64 


8 

8 


3 

6 


23 


0 

0 


51 

54 


T010101 MCC LINER 


1 
1 

1 ,-18 
-10 
•18 


5.2 


5.2 

5.2 

5.2 


21 

1401 

1402 

1403 

1404 

1422 
14 1 1 

1412 

1413 

1414 

1423 
1421 

1 5 


2 

3 

4 -10 5.2 

22 -IB 5.2 
-18 
-18 
-18 

14 -18 

23 -18 5.2 
18 5.2 

5.24 


11 

12 

13 


5.2 

5.2 

5.2 

5.2 


16 


5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5«24 


1 .35HO 
5.40 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.4 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 

5.32 


CYCLE101 

0.0 0.20 


tlNCR-13 AMBIENT TMAX» 70 


0 

$ 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 

-.02 


0.20 

0,20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 


6 

8 

8 

1 

5 

1 

1 

4 

7 

j. 

0 49 

0 S3 

0 55 

8 

8 

8 . 

3 

2 b 

5 • • 0 

8 0 

4 

50 

53 

56 

8 

8 

22 

3 0 

6 0 

51 

54 

1 

11 

1 4 

fa 


13 

14 


23 

0 


34 








1 

1 

1 

1 

9 9 

-1 

72 0 

64 



- 


35 







0 

21 










T010101 

2 

MCC 

LINER 

1 .35HG 

CVCLE101 

+1.INCR.01 

T« 

1.4 

TMAX.-360 

s 

TD10101 

2 

MCC 

liner 

1 .35HG 

CYCLE101 

+1.INCR=02 

T* 

1.6 

TMAX* 488 

s 

TD10101 

2 

MCC 

LINER 

1 .35HG 

CYCLEIOI 

+1,INCR=03 


2.4 

TMAX* 718 

$ 

T010101 

2 

MCC 

liner 

1 .35HG 

CYCLE 101 

+1.INCR.04 

Ts 

3.4 

TMAX.1157 

s 

TD10101 

2 

MCC 

LINER 

1 .35HG 

CYCLEIOI 

+1.INCR.05 

T» 

4.9 

TMAX»1485 

. $ 

TD10101 

2 

MCC 

LINER 

J, ,35HG... 

.CYCLEIOI. 

+1,INCR=06 

T = 

20.5 

TMAX.1484 . 


TD10101 

MCC 

LINER 

1 .35HG 

CYCLEIOI 

•M,INCR = 00 

T=395. 

TMAX=1191 

$ 


223 



2 

T010101 MCC LINER 1.35HG 
2 

TOIOIOi MCC LINER 1.35HG 
2 

TOIOIOI ftICC LINER 1 .3SHG 
2 

TOIOIOI ««CC LINER t .35HQ 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER l.35HG^ 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1 .35HO 
2 

TOIOIOI MCC LINER 1 .35HG 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1 .35HG 
2 

TOIOIOI MCC LINER 1.35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

T010101 MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1.35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1.35HG 
2 

TOIOIOI MCC LINER 1 .35HG 
21 

TOIOIOI MCC LINER 1.35HG 

3 1 5.2 5.40 

22 -18 5.2 5.40 

23 -18 5.2 5.40. 


CYCLElOl +t,INCR»10 T»395.6 

CTCLE101 -m/inCR*!! T»397.1 

CYCLElOl +1.INCR«12 T 400. 

CYCLElOl -f1,INCR>13 AMBIENT 

CYCLElOl +2.INCR*01 T« 1.4 

CYCLElOl +2,INCR«02 T« 1.6 

CYCLElOl +2,INCR*03 T» 2.4 

CYCLElOl +2,INCR«04 T« 3.4 

CYCLElOl +2,INCR=05 T» 4.9 

CYCLElOl -fb.INCR-Oe T« 20.5 

CYCLElOl +2,INCR=09 T*395. 

CYCLElOl +2.INCR310 T-39S.6 

CYCLElOl +2.INCR=11 T*397.1 

CYCLElOl +2,INCR=12 T 400. 

CYCLElOl +2,INCR=>13 AMBIENT 

CYCLElOl +3,INCR*01 T* 1.4 

CYCLElOl +3,1NCR=02 T» 1.6 

CYCLElOl -f3.1NCR-03 Ta 2.', 

CYCLElOl +3,INCR=04 T« 3,4 

CYCLElOl +3.INCH*05 T= 4.9 

CYOLE101 4-3.1NCR»06 T» 20.5 

CYCLElOl ■f3,INCR>09 Ta395. 

CYCLElOl ■•■3.INCRX10 Ts395.6 

CYCLElOl +3,INCR«11 T*397.1 

CYCLElOl +3,INCR*12 T 400. 

CYCLElOl +3,INCR*13 AMBIENT 

CYCLElOl +4.1NCR»01 T« 1.4 

CYCLElOl +4,INCR»02 T« 1.6 
0.0 0.20 
0.0 0 . 20 ' 

0.0 0.20 


.21 

TOIOIOI MCC LINER 1 .35HG CYCLElOl +4,INCR*63 Ts 2.4 


■ TMAX= 858 • $ 

TMAX«-297 $ 

TMAX.s-405 . $ 

TMAXa 70 $ 

TMAX*-360 S 

TMAX- 488 $ 

TMAX* 718 $ 

TMAX»1157 $ 

TMAX*1465 S 

TMAX a 1484 S 

TMAXa 11 91 S 

TMAX* 858 S ■ 

TMAXa-297 $ 

TMAXa-405 $ 

TMAXa 70 $ 

TMAXa-360 ■ $ 

TMAXa 488 . • S 

TMAXa 718 $ ■ 

TMAXa1157 $ 

TMAX a 1465 S 

TMAX =1484 . $■ 

TMAXa 11 91 . $ . 

TMAXa 858 ,• $ 

TMAX *-297 S 

TMAX=-405 $ 

TMAXa 70 $ 

TMAX=-360 $ 

TMAXa 480 • $ 

tMAXa' 718' $■ * 


224 



T010101 hcC LINER 1.35HG CTCLE101 44,1NCR*04 T« 3.4 TMAX=1157 $ 

21 

T010101 MCC LINER 1.3SHG CYCLE101 +4,1NCR*05 T« 4.9 TMAX=146S S 

21 . 

T010101 MCC LINER 1.35HG CYCLE101 44.1NCRs06 T> 20.5 TMAX*1494 . S 


3 1 

3 1 

5.2 

5.40 

0.0 

0.20 

1 *10 

5.2 

5.4 

0.0 

0.20 

2 -18 

5.2 

5.4 

0.0 

0.20 

3 -18 

5.2 

5.4 

0.0 

0.20 

A -18 

5.2 

5.4 

0.0 

0.20 

22 -10 

5.2 

5.4 

0.0 

0.20 

11 -10 

5.2 

5.4 

0.0 

0.20 

12 -10 

5.2 

5.4 

0.0 

0.20 

13 -10 

5.2 

5.4 

0.0 

0.20 

14 -10 

5.2 

5.4 

0.0 

0.20 

23 -18 

5.2 

5.4 

0.0 

0.20 

21 -10 

5.2 

5.4 

0.0 

0.20 

1401 

5.24 

5.32 

-.02 

.06 

1402 

5.24 

5.32 

-.02 

.06 

1403 

5.24 

5.32 

-.02 

.06 

1404 

5.24 

5.32 

-.02 

.06 

1422 

5.24 

5.32 

-.02 

.06 

1411 

5.24 

5.32 

-.02 

:'06 

1417 

5.24 

5.32 

-.02 

.06 

1413 

5.24 

5.32 

-.02 

.06 

14 14 

5.24 

5.32 

-.02 

.06 

1423 

5.24 

5.32 

-.02 

.06 

1421 

5.24 

5.32 

-.02 

.06 


1 5 16 



1 2 


3 

4 

22 


6 

1 0 49 

8 

2 0 

50 e 

3 0 

51 

8 

4 0 52 

8 

5 0 

53 8 

6 0 

. 54 

8 

7 0 55 

8 

8 0 

56 



1 

5 





0 


11 12 


13 

14 

23 • 

. 

1 

1 4 







34 






1 

1 9 9 

-1 . 

72 0 

64 



1 

1 





0 


35 






21 







TD1010 1 

MCC LINER 1.35HG 

CYCLE101 

•*-4,lNCR=09 

Tr395. 

TMAXa 1 191 

$ 

21 






TD1010 1 

MCC Liner i.35hg 

CYCLE101 

+4.INCRa10 

Ts395.6 

TMAX* 858 

$ 

21 






T010101 

MCC LINER 1.35HG 

CYCLE101 

+4,INCR=11 

T=397.1 

TMAXs-297 

$ 

21 






TD10101 

MCC LINER 1.35HG 

CYCLE101 

+ 4,INCR*12. 

T 400. 

TMAX*-405 

$ 

3 

1 5.2 5.40 

0.0 0.20 




3 

1 






\ 

5 16 






8 

1 2 


3 . 

4 

22 


1 0 49 

8 

2 0 

50 6 

3 0 

51 

8 

4 0 52 

8 

5 0 

53 8 

6 0 

54 

8 

7 0 55 

8 

8 0 

56 



1 

5 





0 


11 12 


13 

14 

23 


1 

1 4 





^ 

„_34 






1 

1 9 9 


72 0 

64 ’ 

; 7 


1 

1 





0 
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j5 

21 


TOIOtOI MCC 

LINER 

1 . 35 HG 

CYCLE 101 +4 

4 -1 

1.0 

3 1 

5.2 

S .40 

0.0 

0.20 

3 . 1 

1 -18 

5.2 

5.4 

0.0 

0.20 

2 -18 

5.2 

5.4 

0.0 

0.20 

3 -18 

5.2 

5.4 

0.0 

0,20 

4 -18 

5.2 

5.4 

0.0 

0.20 

22 -18 

5.2 

5.4 

0.0 

0.20 

11 -18 

5.2 

5.4 

0.0 

0.20 

12 -18 

5.2 

5.4 

0.0 

0.20 

13 -18 

5.2 

5.4 

0.0 

0.20 

14 -18 

5.2 

5.4 

0.0 

0.20 

23 -18 

5.2 

5.4 

0.0 

0.20 

21 -18 

5.2 

5.4 

0.0 

0.20 

1401 

5.24 

5.32 

-.02 

.06 

1402 

5.24 

5.32 

-.02 

.06 

1403 

5.24 

5.32 

-.02 

..06 

1404 

5.24 

5.32 

-.02 

.06 

1422 

5.24 

5.32 

-.02 

.06 

14 1 1 

5.24 

5.32 

-.02 

.06 

1412 

5.24 

5.32 

-.02 

.06 

14 13 

5.24 

5.32 

-.02 

.06 

1414 

5.24 

5.32 

-.02 

.06 

1423 

5.24 

5.32 

-.02 

.06 

1421 

5.24 

5.32 

-.02 

.06 

1 5 16 

1 

8 1 

0 

2 

49 

8 

3 

2 

8 4 

0 

52 

8 

5 

0 7 

0 

55 

8 

B 

1 5 

11 


12 


13 

1 1 4 

34 

1 1 

9 

9 

-1 

72 

1 1 




35 

35 


INCR-13 AMBIENT TMAX- 70 $ 


. ^ 22 

0 50 8 3 0 51 

0 53 ‘8 G 0 54 

0 56 

■ 0 

14 23 


b 64 

0 
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'gPXST, T7l . 

USER(Q , 

CHARGE ( » » 

ACCTIHAVEN 000 7949t27*05Sv _ )_ 

CET,PR0CPMR/UN>0CB28. 

ATTACH, TAPEI O-THOO lOt. 

DEFINE , tapes «EH20X 13/CT.PU. PXSTAftT 

BEGIN, APSAC, PROCPNR.NAPSAC. 


PXSTART - CARD IMAGES FOR APSAPLOT USED TO 
£XIBARQiAIE 3RD, 4TH, AND 5TH DUTY-CYCLE 
RESULTS TO 20TH CYCLE CONFIGURATION 



'^Hobioi I»^,C LINEM 1.35 HO CTClIoOI ,IHCR »01 T« 1.4 
2 

TH00101 WCC LINER 1*35 MO CYCLE001 *INCR<02 T* 1.6 

TH00101 MCC LINER 1.35 HO CYCLE001 ♦iWCR^Oa T« 2.4 

2 

TH00101 WCC LINER 1.35 HO CTCLEOQI tlNCR«04 T« 3.4 
2 

TH00101 KCC LINER 1 .35 HG CYCL EOOi ilNCR*05 T*‘~4V9 
2 

TH00101 NCC LINER 1.35 HO CYCLE001 ,1NCR«06 T* 20.5 

2 

TH001QI MCC LINER 1.35 MO CYCLEOOI 4lNCR»09 fr395. 

2 

tHOOIOI WCC LINER 1.35 HG CYCLEOOI .iNCRalO T«395*6 

2 

tHooioi MCC Liner 1.35 hg cycleooi TxHCR*ri“f* 39771 
2 

THOOIOI MCC LINER 1.35 HO CYCLEOOI ,1NCR«1S T«400. 

2 

tMOCIOl MCC LINER 1.35 HQ CYCLEOOI 7lNCR«13 T*1NF. 

12 4 

34 35 

^ , V 9 9 .-1 73 a 64 

THOOIOI MCC LINER 1:35 HO CYCLEOOI -M.1NCR.01 T« 1.4 
2 

THOOIOI MCC L*NER 1 .35 HJ1 CYCLE0014 1 .INCR.02 T* 1.6 

THOOIOI MCC LINER 1.35 HG CYCLEOOI -M.INCR.03 T« 2.4 

2 

THOOIOI MCC LINER_1[.35 HG CYCLEOOI •M.INCR.04 T« 3.4 

2 - -- _ , 

THOOIOI MCC LINER 1.35 HQ CYCLEOOI -ft.lNCRsOS T* 4.9 

THOOIOI MCC LINER 1.35 HG CYCL EOOI -M . INCR *06 T* 20.5 

THOOIOI MCC LINER 1.35 HG CYCLEOOI ^1.lNCRx09 T«395. 

2 

THOOIOI MCC LINER 1.35 HG CYCLEOOI -f1.INCR«10 T::395.6 

THOOIOI MCC LINER 1.35 HQ CYCLEOOI 4>1«INCRs11 T*397.l 

2 

THOOIOI BfICC LINER 1.35 HO CYCLEOOI >1,lNCR>12 T>400. 

THOOIOI MCC LINER 1.35 HO CYCLEOOI +1.1NCR-13 T.INE. 

12 4 

34 35 

1 9 • 9 -I 72 0 64- 

THOOIOI MCC LINER 1.35 HG CYCLEOOI -f2.INCR*01 T* 1.4 

THOOIOI WCC LINER 1.35 HO CYCLEOOI +2, INCR*02 T« " T.8 

THOOIOI MCC LINER 1.35 HG CYCLEOOI +2,INCR«C3 T« 2.4 

THOOIOI MCC LINER t.35 HO CYCLEOOI +2.1NCR«04 ' 3.4 ‘ — 

THOOIOI MCC LINER 1.35 HO CYCLEOOI +2,INCR«05 T« 4.9 

^T H OOIOI wee li n e r > .39 hu ucleoui y 2 .i i xn « w t « 2 0 . 9 
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APPENDIX E 


LOCAL HOT SPOT ANALYSIS REFINEMENT STUDY 

The primary effort of this study was analysis of the basic SSME MCC liner and 
evaluation of life-enhancement designs. An additional task was evaluation of 
increasing the local heating to account for local hot-spot effects This 
local heating was accounted for in the modeling by simply increasing the local 
heat transfer rate, obtaining the resultant temperatures, and analyzing them 
in the 2D plane strain model. This methodology inherently assumes that the 
local heating occurs over a significant band of material axially and complete- 
ly around the hoop direction. For the thermal analysis, this may be an appro- 
priate methodology for obtaining maximum temperatures; but it is only an esti- 

stress and strain picture. In reality, and previously 
noted, the mid-channel cracking predominantly occurs in hot-spot areas. These 
hot-spots roughly take on the shape of a 3:1 or A:1 ellipse, with the long 
widths^^°"^ meridian, and the short axis covering from 5- to 10-channel 


A pctorial view of a developed liner with a hot spot near the throat through- 
out a duty cycle is shown in Fig. 112. The liner has two axisymmetric snr- 
aces with axial land members. This results in an orthotropic stiffness path 
that makes a hot-spot analysis more difficult. During the prechill portion of 
the duty cycle and up to say 600 or 700 F, the liner wall temperatures and 
relatively uniform with little variation from local injector 
effects. As the temperature rises, the hot-spot markedly affects the stiff- 
nesses and thermal loads. At steady state, the hot-spot is essentially an 

elliptical hole to the overall liner since its stiffness is much less than the 
surrounding zone of material. A local hole will result in a distorted shape 

location and strain concentrations at its edges. 
Within the hot-spot, the hole distortion is additive to the local hot-spot 

strains. As the engine is throttled back, and goes through the ensuing cutoff 
conditions, the liner cools and again becomes essentially a stiff membrane 
that locks in a portion of the strain distortions from the hot-spot. This 
hot-spot duty cycle has the potential to both add in large local strains in 
the axial and hoop directions over an axisymmetric model analysis, and add a 
biased strain to each duty cycle, 

A detail model of a total liner to develop this hot-spot duty-cycle strain 
history is a large study in manhours and computer dollars. A first approxima- 
tion, though, was made using a 2D plane strain analysis of one-quarter of a 

anddSt^’c^cle symmetry assumed. Figure 113 pictorially shows the model 

An approximately 8-inch square model was used. For this first-cut analysis 
-^^® modeled, although a refinement of the model would be to 

add in a variable thickness to simulate added stiffness from the land. The 

modeled as a stepped variation in x, y dimensions to approximate 

ellipse. A 4 X 1-inch hot-spot was simulated. Symmetry boundary 
conditions were applied along x=0 and y=0. Fixed displacements were added tl 
the other boundaries to simulate the hoop and axial displacements. 
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Figure 112. 3D Hot-Spot Effects 












The duty cycle was reduced to seven load increments, as noted, and constant 
temperatures applied in two zones — one inside the hot-spot and the remaining 
area on average liner temperature. Only one duty cycle was run as a test. 
Figures 114 through 116 are plots of the steady-state condition; and Fig. 117 
and 118 shows the residual strains after cooldown at 70 F condition. At both 
steady state and the cooldown condition, the peak hoop strain in the hot spot 
zone is about 6 or 7 times the strain level away from the hot spot. The axial 
strain is about the same throughout the model. From a strain range basis, the 
peak strain range is 3.4%, which is greater than twice the strain range away 
from the hot-spot. A bias of almost -3% occurs in the maximum strain condi- 
tion after cooldown (Table 13). 

Some sensitivity studies were made early in the project on adjusting axial 
strains and hoop effects by using pseudo AT effects. These results showed 
some reduction in life or more rapid thinning, but the results were not con- 
sistent with the strain bias or large factors herein shown on hoop strains. 

Future efforts on understanding local hot-spot effects should be slanted 
toward more accurate 3D representations, and possibly more 2D simplified anal- 
ysis may well explain the difference in actual local hot-spot life in test vs 
the current axisymmetric (total hot spot) analysis. 

TABLE 13. COMPARISON OF HOT-SPOT STRAINS TO 
NOMINAL CONDITION STRAINS 


CONDITION 

RATIO 

PEAK TO AVERAGE c 

PEAK 

AVERAGE 

HOOP, Y 

AXIAL, X 

HOOP, % 

AXIAL, % 

HOOP, % 

AXIAL, % 

STEADY-STATE 

6 

1 .25 

-6.2 

-1 .25 

-1 

-1 

294 K (70 F) 

-7 

.8 

-2.8 

0.65 

0.4 

0.8 

DIFFERENCE 



-3.4 

-1.9 

-1.4 

-1.8 
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APPENDIX F 


GEOMETRY OF MCC COOLANT LINER DESIGNS 

The MCC coolant liner geometry is presented for each of the designs analyzed 
ratio length. The hot-gas wall contour is shown by the area- 
neso length. The coolant geometry (channel, land, and wall thick- 

^ presented as a function of axial location. The con- 

bas^ SSMe"mcc life design is not shown because it is the same as the 

oasic SSME-MCC design, except as previously noted within the text. 
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TABLE 14. BASELINE SSHE MCC COOLING CHANNEL DIMENSIONS 
(390 CHANNELS) IN SI UNITS 



"j 

LAND 

WIDTH 

HEIGHT 

WALL 

(cm) 


(cm) 

(cm) 

(cm) 

(cm) 

24.71 

5.0 

0.3155 

0,1575 

0.6096 

0.08 8 9 

20.32 

4.052 

0.2685 

0.1575 

0.6096 

0.0889 

17.78 

3.538 

0.2405 

0.1575 

0.6096 

0.0889 

15.24 

3.066 

0.2131 

0.1575 

0.5817 

0.0889 

12.70 

2.608 

0.1844 

0.1575 

0.5476 

0.0889 

10.16 

2.182 

0.1554 

0.1575 

0.5100 

0.08 8 9 

7.62 

1.782 

0.1603 

0.1227 

0.5004 

0.0889 

5.08 ' 

1.426 

0.1516 

0.1016 

0.4953 

0.08 8 9 

2.54 

1.103 

0.1209 

0.1016 

0.4318 

0.0711 

j 0.0 

1.000 

0.1105 

0.1016 

0.2870 

0.0711 

j -1.27 

lx}008 

0.1113 

0.1016 

0.2741 

0. 0711 

j -2.54 

1.038 

0.1143 

0.1016 

0.2609 

0.0711 

1 -3.0A8 

1.055 

0.1261 

0.1016 

0.2558 

0.0711 

-3.81 

1,090 

0.1196 

0.1016 

0.2477 

0. 0711 

-5.08 

1.169 

0.1275 

0.1016 

0.2352 

0.0711 

-6.35 

1.259 

0.1361 

0.1016 

0.2217 

0.0711 

-7.62 

1.369 

0.1463 

0.1016 

0.2139 

0.0711 

-8.89 

1.481 

0.1562 

0.1016 

0.2337 

0.0711 

-10.16 

1.590 

0.1653 

0.1016 

0.2535 

0.0711 

-12.70 

1.831 

0.1847 

0.1016 

0.2845 

0.0711 

-15.24 

2.085 

0.1844 

0.1217 

0.2619 

0.0889 

-17.78 

2.362 

0.1748 

0,1509 

0.2512 

_P-5889_ 

-20.32 

2,628 

0.1857 

0.1575 

0.2489 

0.0889 

-22.86 

2.812 

0.1976 

0.1575 

0.2489 

0.0889 

-25.40 

2.924 

0.2045 

0.1575 

0.2489 

0,0889 

-27.94 

2.962 

0.2068 

0.1575 

0.2489 

0.0889 

-30.48 

2.962 

0.2068 

0.1575 

0.2489 

0.0889 

-32.41 

2,962 

0.2068 

0.1575 

0.2489 

0.0889 

-34.29 

2.962 

0.2068 

0.1575 

0.5842 

0.0889 

-35.56 

2.962 

0.2068 

0,1575 

0.635 

0.0889 
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TABLE 14. (Concluded) 

BASELINE SSME MCC COOLING CHANNEL DIMENSIONS 
(390 CHANNELS) IN ENGLISH UNITS 




LAND 

WIDTH 

HEIGHT 

mALl 

(in) 

(in) 

(in) , 

Cln) 

(In) 

9.73 

5.0 

0.1242 

0.062 

0.24 

0*o45 

8.0 

4.052 

0.1057 

0.062 

0.24 

0.035 

7.0 

3,538 

0.0947 

0.062 

0.24 

0.035 

6.0 

3.066 

0.0839 

0.062 

0.229 

0.035 

5.0 

2.608 

0.0726 

0.062 

0.2156 

0.035 

4.0 

2.182 

0.0612 

0.062 

0.2007 

0.035 

3.0 

1.782 

0.0631 

0.0483 

0.1970 

0.035 

2.0 

1.426 

0.0597 

0.04 

0.1950 

0.035 

1.0 

1.103 

0.047 6 

0.04 

0.17 

0.028 

0.0 

1.000 

0.0435 

0.04 

0.113 

0.028 

-0.5 

1.008 

0.0438 

0.04 

0.1079 

0.028 

-1.0 

1.038 

0.0450 

0.04 

0.1027 

0.028 

-1.2 

1.055 

0.0457 

0.04 

0.1007 

0.028 

-1.5 

1.090 

0.0471 

0.04 

0.0975 

0.028 

-2.0 

1.169 

0.0502 

0.04 

0.0926 

0.028 

-2.5 

1.259 

0.0536 

0.04 

0.0873 

0.028 

-3.0 

1.369 

0.0576 

0.04 

0.0842 

0.028 

-3.5 

1.481 

0.0615 

0.04 

0.0920 

0.028 

-4.0 

1.590 

0.0651 

0.04 

0.0998 

0.028 

-5.0 

1.831 

0.0727 

0.04 

0.1120 

0.028 

-6.0 

2.085 

0.0726 

0.0479 

0.1031 

0.035 

-7.0 

2.362 

0.0688 

0.0594 

0.0989 

0.035 

-8.0 

2.628 

0.0731 

0.062 

0.098 

0.035 

-9.0 

2.812 

0.0778 

0.062 

0.098 

0.035 

-10.0 

2.924 

0.0805 

0.062 

0.098 

0.035 

-11.0 

2.962 

0.0814 

0.062 

0.098 

0.035 

-12.0 

2.962 

0.0814 

0.062 1 

0. 098 

0.035 

-12.76 

2.962 

0.0814 

0.062 

0.098 

0.035 

-13.50 

2.962 

0.0814 

0.062 

0.23 

0.035 

-14.0 

2.962 

0.0814 

0.062 

0.25 

0.035 
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TABLE 15. INCREASED NUMBER OF COOLANT CHANNELS DIMENSIONS 
(540 CHANNELS) IN SI UNITS 


xJ 

(cm) 

Ej 

LAND 

(cm) 

WIDTH 

(cm) 

HEIGHT 

(cm) 

WALL 

(cm) 

24.71 

5.0 

0.200 

0.142 

0.610 

“07(^13 

20.32 

4.052 

0.165 

0.142 

0.610 

0.0813 

17.78 

3.538 

0.145 

0.142 

0.610 

0.0813 

15.24 

3.066 

0.137 

0.130 

0.610 

0.0813 

12.70 

2.608 

0,129 

0.118 

0.555 

0,0813 

10.16 

2.182 

0,120 

0,106 

0.500 

0.0813 

7.62 

1.782 

0.110 

0.094 

0,446 

0.0813 

5.08 

1.426 

0.101 

0.082 

0.391 

0.0813 

2.54 

1.103 

0.084 

0.076 

0.337 

0.0635 

0.0 

1,000 ^ 

0.077 

0.076 

0.282 

0.0635 

-1.27 

1.008 

0.077 

0.076 

0.269 

0.0635 

-2.54 

1.038 

0.080 

0.076 

0.257 

0.0635 

-3.048 

1.055 

0.081 

0.076 

0.257 

0.0635 

-3.81 

1.090 

0.084 

0.076 

0.257 

0.0635 

-5.08 

1.169 

0.089 

0.076 

0.257 

0.0635 

-6.35 

1.259 

0.096 

0.076 

0.257 

0.0635 

-7.62 

1.369 

0.103 

0.076 

0.257 

0.0635 

-8.89 

1.481 

0.110 

0.076 

0.269 

0.0635 

-10.16 

1.590 

0.117 

0,076 

0.282 

0.0635 

-12.70 

1.831 

0.131 

0.076 

0.307 

0.0635 

-15.24 

2.085 

0.126 

0.095 

0.287 

0.0762 

-17.78 

2.362 

0.122 

0,113 

0.26? 

0.0762 

-20.32 

2.628 

0.131 

0.117 

0.257 

0.0762 

-22.86 

2,812 

0.139 i 

0.117 

0.257 

0.0762 

-25.40 

2.924. 

0.145 

0.117 

0.257 1 

0.076 

-27.94 

2.962 

0.146 

0.117 

0.257 

0.0762 

-30.45 

2.962 

0.146 

0.117 

0.320 

0.0762 

-32.42 

2.962 

0.146 

0.117 

0.369 

0,0762 

-34.29 

2.962 

0.146 

0,117 

0.584 

0.0762 

-35.56 

2.962 

0.146 

0.117 

0.635 

0.0762 


NOMENCLATURE: 


xj - AXIAL DISTANCE FROM THROAT 

Ej - CHAMBER AREA RATIO 

LAND - LAND WIDTH BETWEEN CHANNELS 

WIDTH - CHANNEL WIDTH 

HEIGHT - CHANNEL HEIGHT 

WALL - HOT GAS WALL THICKNESS 




TABLE 15. (Concluded) 


INCREASED NUMBER OF COOLANT CHANNELS DIMENSION 
(540 CHANNELS) IN ENGLISH UNITS 


Xj 

( In ) 

EJ 

LAND 

( In ) 

WIDTH 

( in ) 

HEIGTH 

( In ) 

WALL 

( In ) 

9.73 

5.0 

0.0784 

0.056 

0,24 

0.032 

8.0 

4.052 

0.0651 

0.056 

0.24 

0.032 

7.0 

3.538 

0.0571 

0.056 

0.24 

0.032 

6.0 

3.066 

0.0541 

0.0513 

0,24 

0.032 

5.0 

2.608 

0.0507 

0.0465 

0.2185 

0.032 

4.0 

2.182 

0.0471 

0.0418 

0.1970 

0.032 

3.0 

1.782 

0.0433 

0.0371 

0.1755 

0.032 

2.0 

1.426 

0.0396 

0.0324 

0.1540 

0.032 

1.0 

1.103 

0.0332 

0.030 

0.1325 

0.025 

0.0 

1.000 

0.0302 

0.030 

0.1110 

0.025 

- 0.5 

1.008 

0.0305 

0.030 

0.1060 

0.025 

- 1.0 

1.038 

0.0314 

0.030 

0.1010 

0.025 

- 1.2 

1.055 

0.0319 

0.030 

0.1010 

0.025 

- 1.5 

1.090 

0.0329 

0.030 

0.1010 

0.025 

- 2.0 

1.169 

0.0351 

0.030 

0,1010 

0.025 

- 2.5 

1.259 

0.0376 

0.030 

0.1010 

0.025 

- 3.0 

1.369 

0.404 

0.030 

0.1010 

0.025 

- 3.5 

1.481 

0.0433 

0.030 

0.1060 

0.025 

- 4.0 

1.590 

0.0459 

0.030 

0.1110 

0.025 

- 5.0 

1.831 

0.0514 

0.030 

0.1210 

0.025 

- 6.0 

2.085 

0.0497 

0,0373 

0.1130 

0.030 

- 7.0 

2.362 

0.0480 

0.0445 

0.103 

0.030 

—8 . 0 

2.628 

0. 0515 

0.04 60 

0.1010 

0.030 

- 9.0 

2.812 

0.0549 

0.046 

0.1010 

0.030 

- 10.0 

2.924 

0.0569 

0.046 

0.1010 

0.030 

- 11.0 

2.962 

0.0575 

0.046 

0.101 

0.030 

- 12.0 

2.962 

0.0575 

0.046 

0.1260 

0.030 

- 12.76 

.2962 

0.0575 

0.046 

0.1451 

0.030 

- 13.50 

2.962 

0 . 05*'5 

0.046 

0.23 

0.030 

- 14.0 

2.962 

0.0575 

0.046 

0.25 

0.030 


NOMENCLATURE; 

xj 

- AXIAL DISTANCE FROM THROAT 


Ej 

- CHAMBER AREA RATIO ; 


LAND 

- LAND WIDTH BETWEEN CHANNELS 


WIDTH 

- CHANNEL WIDTH 


HEIGHT 

- CHANNEL HEIGHT 


WALL 

« HOT GAS WALL THICKNESS 




TABLE 16. KEEL-RIB CHANNEL DIMENSION 
(390 CHANNELS) IN SI UNITS 



(cm) 


LAND 

(cm) 

WIDTH 

(cm) 

HEIGHT 

(cm) 

WALL 

(cm) 

F.H./C.H. 

RATIO 

FIN BASE 
(cm) 


2A.71 

5.0 

0.3155 

0.1575 

0.6096 

0.0889 

N.A, 

N.A. 


20.32 

4.052 

0.2685 

0.1575 

0,6096 

0.08 8 9 

N.A. 

^N.A. 


17.78 

3.538 

0.2405 

0.1575 

0.6096 

0.08 8 9 

N.A. 

N.A. 


15.24 

3.066 

0.2131 

0.1575 

0.5817 

0.08 8 9 

N.A. 

N.A. 


12.70 

2.606 

0.1844 

0.1575 

0.5474 

0.0889 

N.A. 

N.'A. 


10.16 

2.182 

0.1554 

0.1575 

0.5093 

0.0889 

N.A. 

N.A. 


7.62 

1.782 

0.1603 

0.1227 

0.5352 

0.0889 

N.A. 

N.A. 


5.08 

1,426 

0.1516 

0.1016 

0.4953 

0,0889 

N.A. 

N.A. 


2.54 

1.103 

0.1209 

0.1016 

0.4318 

0.0711 

N.A. 

N.A. 


0.0 

1.000 

0.1105 

0.1016 

0.3442 

0.0711 

0.3 

0.0339 


-1.27 

1.008 

0.1U3 

0.1016 

0.3307 

0.0711 

0.3 

0.0339 


-2.54 

1.038 

0.1143 

0.1016 

0,3147 

0.0711 

0.3 

0.0339 


-3.048 

1.055 

0.1261 

0.1016 

0.3117 

0.0711 

0,3 

0.0339 


-3.81 

1.090 

0.1196 

0.1016 

0.2982 

0.0711 

0.3 

0.0339 


-5.08 

1.169 

0.1275 

0.1016 

0.2977 

0.0711 

0.3 

0.0339 


-6.35 

1.259 

0.1361 

0.1016 

0,2812 

0.0711 

0.3 

0.0339 


-7.62 

1.369 

0,1463 

0.1016 

0.2728 

0.0711 

0.3 

0.0339 


-8,89 

1.481 

0.1562 

0.1016 

0.3987 

0.0711 

0.3 

0.0339 


-10.16 

1.590 

0.1653 

0.1016 

0.3231 

0.0711 

0.3 

0.0339 


-12,70 

1.831 

0.1847 

0.1016 

0.4338 

0.0711 

0.3 

0.0339 


-15.24 

2.085 

0.1844 

0.1217 

0.3792 

0.0889 

0.3 

0.0406 

-17.78 

2.362 

0.1748 

0.1509 

0.3284 

0.0889 

0.3 

0.0503 

-20.32 

2.628 

0.1857 

0.1575 

0.3175 

0.0889 

0.3 

0.0525 

-22.86 

2.812 

0.1976 

0.1575 

0.3175 

0.0889 

0.3 

0.0525 

-25.40 

2.924 

0.2045 

0.1575 

0.3175 

0.08 8 9 

0.3 

0.0525 

-27.94 

2.962 

0.2068 

0.1575 

0.3175 

0.08 8 9 

0.3 

0.0525 

-30.48 

2.962 

0.2068 

0.1575 

0.2489 

0.0889 

N.A. 

N.A. 

-32.41 

2.962 

0,2068 

0.1575 

0.2489 

0,0889 

N.A. 

N.A. 

-34.29 

2.962 

0.2068 

0.1575 

0.5842 

0.0889 

N.A. 

N.A. 

-35.56 

2.962 

0.2068 

0.1575 

0.635 

0.0889 

N.A. 

N.A. 
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TABLE 16. (Concluded) 


KEEL-RIB CHANNEL DIMENSIONS 
(390 CHANNELS) IN ENGLISH UNITS 


(In) 


LAND 

(In) 

WIDTH 

(In) 

HEIGHT 

(In) 

WALL 

(In) 

F.H./C.H. 

RATIO 

FIN BASE 
(in) 

9.73 

5.0 

0.1242 

0.062 

0.24 

0.035 

N.A. 

N.A. 

8.0 

4.052 

0.1057 

0.062 

0.24 

0.035 

N.A. 

N.A. 

7.0 

3.538 

0.0947 

0.062 

0.24 

0.035 

N.A. 

N.A. 

6.0 

3.066 

0.0839 

0.062 

0.229 

0.035 

N.A. 

N.A. 

5.0 

2.068 

0.0726 

0.062 

0.2155 

0.035 

N.A. 

N.A. 

4.0 

2.182 

0.0612 

0.062 

0.2005 

0.035 

N.A. 

N.A. 

3.0 

1.782 

0.0631 

0.0483 

0.2107 

0.035 

N.A. 

N.A. 

2.0 

1.426 

0.0597 

0.04 

0.1950 

0.035 

N.A. 

N.A. 

1.0 

1.103 

0.0476 

0.04 

0.1700 

0.028 

N.A. 

N.A. 

0.0 

1.000 

0.0435 

0.04 

0.1355 

0.028 

0.3 

0.0133 

-0.5 

1.008 

0.0^38 

0.04 

0.1302 

0.028 

0.3 

0.0133 

-1.0 

1.038 

0.0450 

0.04 

0.1239 

0.028 

0.3 

0.0133 

-1.2 

1.055 

0.0457 

0.04 

0.1227 

0.028 

0.3 

0.0133 

-1.5 

1.090 

0.0471 

0.04 

0.1174 

0.028 

0.3 

0.0133 

-2.0 

1.169 

0.0502 

0.04 

0.1172 

0.028 

0.3 

0.0133 

-2.5 

1.259 

0.0536 

0.04 

0.1107 

0.028 

0.3 

0.0133 

-3.0 

1.369 

0.0576 

0.04 

0.1074 

0.028 

0.3 

0.0133 

-3.5 

1.481 

0.0615 

0.04 

0.1176 

0.028 

0.3 

0.0133 

-4.0 

1.590 

0.0651 

0.04 

0.1272 

0.028 

0.3 

0. 0133 

-5.0 

1.831 

0.0727 

0.04 

0.1708 

0.028 

0.3 

0.0133 

-6.0 

2.085 

0.0726 

0.047 9 

0.1493 

0.035 

0.3 

0.0160 

1 

O 

2,362 

0.0688 

0.0594 

0.1293 

0.035 

0.3 

0.0198 

-8.0 

2.628 

0.0731 

0.062 

0.1250 

0.035 

0.3 

0.0207 

-10.0 

2.924 

0.0805 

0.062 

0.1250 

0.035 

0.3 

0.0207 

-11.0 

2.962 

0.0814 

0.062 

0.1250 

0.035 

0.3 

0.0207 

0 

CM 

H 

1 

2.962 

0.0814 

0.082 

0.0980 

0.035 

N.A. 

N.A. 

-12.76 

2,962 

0.0814 

0.062 

0.0980 

0.035 

N.A. 

N.A. 

-13.50 

2.962 

0.0814 

0.062 

0.23 

0.035 

N.A. 

N.A. 

-14.0 

2.962 

0.0814 

0.062 

0.25 

0.035 

1 N.A. 

N.A. 
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TABLE 17 . SLOTTED HOT-GAS WALL COOLANT CHANNEL DIMENSION 
(460 CHANNELS) IN SI UNITS 


1 

(cm) 

Ej 

LAND 

(cm) 

WIDTH 

(cm) 

HEIGHT 

(cm) 

WALL 

(cm) 

2A.71 

5.0 

0.259 

0.142 

0,610 

0.0813 

20.32 

4.052 

0.219 

0.142 

0,610 

0.0813 

17.78 

3.538 

0.195 

0.142 

0.610 

0.0813 

15.24 

3.066 

0.184 

0.130 

0.610 

0.0813 

12.70 

2.608 

0.172 

0.118 

0.563 

0.0813 

10.16 

2.182 

0,159 

0.106 

0.516 

0.0813 

7.62 

1.782 

0.146 

0.094 

0.470 

0.0813 

5.08 

1.426 

0.132 

0.082 

0.423 

0.0813 

2.54 

1.103 

0.112 

0.076 

0.377 

0,0635 

0.0 

1.000 
q- 

0.103 

0,076 

0.330 

0.0635 

-1.27 

1.008 

0.104 

0.076 

0.318 

0.0635 

-2.54 

1.038 

0.107 

0.076 

0.305 

0.0635 

-3.048 

1.055 

0.108 

0,076 

0.300 

0.0635 

-3.81 

1.090 

0,111 

0.076 

0,300 

0.0635 

-5.08 

1.169 

0,118 

0.076 

0.300 

0.0635 

-6.35 

1.259 

0.125 

0.07 6 

0.300 

0.0635 

-7.92 

1.369 

0.134 

0.076 

0.300 

0,0635 

-8.89 

1.481 

0.142 

0.076 

: 0.312 

0.0635 

-10.16 

1.590 

0,150 

0.076 

0.325 

0.0635 

-12.70 

1.831 

0.167 

0.076 

0.351 

0.0635 

-15.24 

2.085 

0.164 

0.095 

0,330 

0.0762 

-17.78 

2.362 

0.163 

0.113 

0.305 

0.0762 

-20.32 

2,628 

0.174 

0.117 

0.300 

0.0762 

-22.86 

2.812 

0.184 

0.117 

0.300 

0.0762 

-25.40 

2.924 

0.190 

0.117 

0.300 

0.0762 

-27.94 

2.962 

0.192 

0.117 

0.300 

0.0762 ' 

-30.45 

2.962 

0.192 

0.117 

0.432 

0.0762 

-32.42 

2.962 

0.192 

0.117 

0.510 

0.0762 

-34.29 

2.962 

0.192 

0.117 

0.584 

0.0762 

-35.56 

2,962 

0.192 

0.117 

0.635 

0.0762 
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TABLE 17. (Concluded) 


SLOTTED HOT-GAS WALL COOLANT CHANNEL DIMENSIONS 
(460 CHANNELS) IN ENGLISH UNITS 



Xj 

(In) 

EJ 

LAND 

(in) 

WIDTH 

(In) 

HEIGHT 

(in) 

WALL 

(in) 


9.73 

5.0 

0.1018 

0.056 

0.24 . 

0.032 


8.0 

4.052 

0.0861 

0,056 

0.24 

0.032 


7.0 

3.538 

0.0768 

0.056 

0.24 

0.032 


6.0 

3.066 

0.0724 

0.0513 

0.24 

0.032 


5.0 

2.608 

0.0676 

0.0465 

0.2217 

0.032 


4.0 

2.182 

0.0626 

0.0418 

0.2033 

0.032 


3.0 

1.782 

0.0573 

0.0371 

0.1850 

0.032 


2.0 

1.426 

0.0521 

0.0324 

0.1667 

0.032 


1.0 

1.103 

0.0442 

0.030 

0.1483 

0.025 


0.0 

1.000 

0.0407 

0.030 

0,130 

0.025 

-0.5 

1.008 

0.0410 

0.030 

0.1250 

0.025 

-1.0 

1.038_^ 

0.0421 

0.030 

0.1200 

0.025 

-1.2 

1.055 

0.0426 

0.030 

0.1180 

0.025 

-1.5 

1.090 

0.0438 

0.030 

0.1180 

0.025 

-2.0 

1.169 

0, 0464 

0.030 

0.1180 

0.025 

-2.5 

1.259 

0.0493 

0.030 

0.1180 

0,025 

-3.0 

1.369 

0.0527 

0.030 

0.1180 

0.025 

-3.5 

1.481 

0.0560 

0.030 

0.1230 

0.025 

-4.0 

1.590 

0.0591 

0.030 

0.1280 

0.025 

-5.0 

1.831 

0.0656 

0.030 

0.1380 

0.025 

-6.0 

2.085 

0.0647 

0.0373 

0.1300 

0.030 

-7.0 

2.362 

0.0641 

0.0445 

0.1200 

0.030 

-8.0 

2.628 

0.0685 

0.046 

0.1180 

0.030 

-9.0 

2.812 

0.0724 

0.046 

0.1180 

0.030 

-10.0 

2.924 

0.0748 

0.046 

0.1180 

0.030 

-11.0 

2.962 

0.0755 

0.046 

0.1180 

0.030 

-12.0 

2.962 

0.0755 

0.046 

0.1700 

0.030 

-12.76 

2.962 

0.0755 

0.046 

0.2006 

0.030 

-13.50 

2.962 

0.0755 

0.046 

0.23 

0.030 

-14.0 

2.962 

0.0755 

0. 046 

0.25 

0.030 
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xj 
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EJ 

« CHAMBER AREA RATIO 
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- CHANNEL WIDTH 
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- HOT GAS WALL THICKNESS 




APPENDIX G 


MID-CHANNEL WALL THINNING ANALYSIS ADJUSTMENT 

The analysis results used to develop the percent of mid-channel thinning vs 
duty cycles are furnished in this appendix for all configurations. Each curve 
has the analysis results plotted as a solid line. The dashed line is the 
adjustment used to eliminate the shakedown conditions that occur during the 
first one or two cycles of each 5 cycle analysis. Figure 119 depicts the 
baseline 135% hg analysis results, shows the shakedown condition that occurred 
in each 5 duty cycles, and the method of adjusting the raw analysis results to 
acount for the shakedown. Other configurations Fig. 120 through 124 and oper- 
ating condition analyzed had similar shakedown conditions. For the first 5 
duty cycles, adjustment of the residual stress and strain conditions during 
the load cycling resulted in a realistic shakedown condition. The first 
shakedown condition also included the gross distortion of the liner as it per- 
manently yields through the 0.51 mm (0.020 inch) tolerance gap and touches the 
jacket. Thereafter, at the steady-state posttest condition, the liner remains 
virtually in line-to-line contact with the jacket, so the first shakedown is 
larger than the remaining ones. The other shakedown condition, where there is 
a rapid change in thickness fo rhe first 1 or 2 cycles, followed by a much 
slower change in cycles 3, 4, and 5, occurs because of the approximations 
made. Each new 5 duty cycles start with a stress-and-strain-free distorted 
geometry rather than were used directly, a significant error in thinning would 
be obtained. The final methodology chosen to eliminate the shakedown effect 
was to assume that the distortions from the first cycles in the 5-cycle analy- 
sis were equivalent to a series of thinning cycles at the rate of change used 
for extrapolation to the restart condition. This is effectively obtained by 
extending the last 5-20-cycle extrapolation line (e.g.. A) until it intersects 
with a horizontal line drawn from the end of the shakedown condition, B. At 
this point, the B slope applies until the C slope horizontal line intersects. 
This technique was used to adjust all the analysis results in this study. 
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Figure 119. SSME MCC Life Prediction, 135% hg Baseline 




ACCUMULATED DUTY CYCLES 


Figure 120. SSME Configuration 135% hg and 100% hg Raw 
Data Thinning vs Accumulated Duty Cycles 





lO 120 140 160 180 200 

ACCUMULATED DUTY CYCLES 

Figure 121. Increased Number of Channels Configuration Raw 
Data Thinning vs Accumulated Duty Cycles 




ACCUMULATED DUTY CYCLES 


Figure 122. Contoured Wall Configuration Raw Data Thinning 
vs Accumulated Duty Cycles 
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Figure 123, Keel Rib Configuration Thinning vs Accumulated Duty Cycles 







APPENDIX H 


TYPICAL CONTOUR AND GRAPHICAL DATASET 

The enclosed figures and curves are a selected subset of the contour and 
graphical plots for the baseline SSME MCC configuration. Plots are presented 
for the first duty cycle, but they are representative of results throughout 
the analysis. The lod increment and times on the graphs are consistent with 
the duty cycle analysis presented in the Structural Modeling section of the 
report. Results are primarily presented for: <1) the steady-state operation 
condition — increment 6 — where maximum compressive strains occur, (2) cut off — 
increment 10 — when the liner is coldest and is at the maximum tension portion 
of the cycle, and (3) after cutoo off steady-state — increment 11 — wqhen the 
MCC is at 70 F, the condition used as a basis for the extrapolation. 

A graphical plot of the stress state around the model hot wall boundary is 
furnished for increments 7 and 9 during the transient phase of cutoo off. 
These depict how the stress state changes from comprehension and bending to an 
average tensile stress while the hot wall f the liner is still hot (600 to 
700 F). 
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